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I.  INTRODUCTION 


The  disposal  of  dredged  material  has  received  much  attention  from 
environmentalists.  Section  404(b)  of  the  Federal  Water  Pollution  Control  Act 
Ammendments  of  1972,  P.L.  92-500,  prohibits  discharges  of  dredged  material  to 
navigable  waters  of  the  United  States  unless  permits  are  issued  through  the 
U.S.  Army  Corps  of  Engineers.  In  1975  guidelines  on  the  issuance  of  permits 
were  published  in  the  Federal  Register.  Among  the  ecological  impacts  from 
dredged  material  disposal  to  be  examined  were  impairment  of  the  water  column 
and  the  covering  of  benthic  communities.  There  exists  a  need  for  mathematical 
models  to  predict  the  concentration  and  deposition  of  suspended  solids  and 
water  quality  constituents  resulting  from  disposal  of  dredged  material. 

The  Corps  of  Engineers  currently  dredges  portions  of  the  upper 
Mississippi  River  to  maintain  a  nine-foot  (deep)  channel  for  barge  traffic. 
The  hydraulically  dredged  material  is  discharged  into  open  water  or  onto  a 
nearby  bank.  This  water  contains  suspended  solids  and  forms  a  suspended- 
solids  plume  and  a  chemical  containment  plume  where  it  enters  the  river. 

A  multidisciplinary  consortium  called  the  Great  River  Environmental 
Action  Team,  GREAT,  was  formed  in  1975  to  determine  the  environmental  effects 
of  dredging  and  dredge  disposal  on  the  Mississippi  River.  During  the  second 
phase  of  this  study,  GREAT  II  (  from  Guttenberg,  Iowa,  to  Saverton,  Missouri), 
two  mathematical  models  were  developed  at  The  University  of  Iowa  Institute  of 
Hydraulic  Research  to  assess  water  quality  concentrations  after  dredge 
disposal.  The  three-dimensional,  steady-state  model  of  Wechsler  and  Cogley 
(the  Walden  Plume  model)  was  modified  for  the  case  of  side-bank  disposal  of 
dredged  material  in  the  GREAT  II  study  reach.  Also  a  two-dimensional,  steady- 
state  model  was  developed  by  means  of  an  analytical  solution  to  the  depth- 
averaged  equation.  It  was  a  recommendation  of  the  GREAT  II  study  that  the 
model  receive  further  refinement  and  verification. 

In  this  research,  the  Wechsler-Cogley  model  is  applied  to  three  sites  in 
the  GREAT  III  study  reach,  from  Saverton,  Missouri,  to  Cairo,  Illinois.  This 
reach  includes  different  environmental  conditions  than  those  encountered  in 
GREAT  II — i.e.,  open-water  disposal  of  hydraulically  dredged  material,  dustpan 
as  well  as  cutterhead  hydraulic  dredging  techniques,  and  conditions  existing 
below  the  lock  and  dam  and  pool  system  of  the  upper  Mississippi  River. 
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Secondly,  the  two-dimensional  analytic  model  is  developed  into  a  chemical 
contaminant  plume  model.  It  is  applied  to  a  site  within  the  GREAT  III  study 
reach.  Chemical  kinetics  and  sorption  partitioning  are  included  in  the 
model.  In  addition,  sediments  from  seven  locations  along  the  GREAT  III  study 
reach  are  screened  via  the  standard  elutriate  test.  Elutriation  kinetics  are 
delineated. 

A  primary  objective  of  this  research  was  to  refine  and  verify  the 
accuracy  of  a  mathematical  model  of  suspended-sediment  dispersion,  and  a  model 
of  desorption  and  dispersion  of  toxic  materials  from  dredged  sediments.  These 
models  were  developed  during  the  GREAT  II  study,  but  additional  sampling  of 
suspended-sediment  plumes  associated  with  dredge  disposal  is  needed  to  test 
the  accuracy  of  the  model  at  different  stream  velocities  and  for  different 
sediment  types.  The  refined  mathematical  formulas  will  provide  users  with  a 
quantitative  method  for  assessing  water-quality  impacts  under  the  Federal 
Clean  Water  Act  for  dredging.  Associated  scope  and  objectives  of  this 

research  were  to: 

1)  Measure  suspended-sediment  plumes  (length  and  width  at  mid¬ 

depth)  associated  with  open-water  disposal  channel  maintenance 
dredging  operations. 

2)  Refine  the  Wechsler-Cogley  particle-dispersion  model  as  modified 
by  The  University  of  Iowa  Institute  of  Hydraulic  Research. 

3)  Collect  water  and  sediment  samples  from  five  designated  harbor 
and  shoreline  areas  on  the  Mississippi  River,  and  analyze  them 
for  ammonia,  phosphorus,  COD,  oil  and  grease,  cadmium,  chromium, 
copper,  lead,  mercury,  nickel,  manganese,  iron,  zinc,  and  PCB's. 

4.  Perform  particle-size  analysis  and  priority-chemical -pollutant 

analysis  for  the  sediments. 

5.  Conduct  standard  elutriate  and  kinetic  tests  to  predict  the 

extent  of  release  of  chemical  pollutants  from  the  sediments  and 
to  determine  the  rate  of  desorption. 
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II.  SUSPENDED  SOLIDS  MODEL 

A.  Hechsler-Coqley  Model  Development.  The  model  with  which  this  chapter 
is  principally  concerned  is  based  on  a  turbidity  plume  computer  model  proposed 
by  Wechsler  and  Cogley  (1977).  It  was  developed  for  predicting  concentrations 
of  suspended  sediment  downstream  from  a  line  source  (e.g.,  a  hydraulic 
pipeline  dredge  discharging  into  open  water).  These  concentrations  are 
predicted  as  a  function  of  sedimentation  data  (from  jar  tests)  and  hydraulic 
parameters  such  as  eddy  diffusivity  and  current  velocity.  The  initial 
simplifying  assumptions  utilized  are  unidirectional  constant  flow,  infinite 
width,  constant  depth,  and  infinite  length. 

The  concentration  of  suspended  sediment  downstream  from  a  disposal  site 
is  determined  by  a  balance  among  four  sediment  transport  mechanisms:  removal 
by  downward  settling  and  deposition,  transport  by  vertical  diffusion  in  the 
direction  of  decreasing  concentration  gradient,  lateral  dispersion  by 
turbulent  diffusion,  and  dispersion  in  the  downstream  direction  by  both 
convection  and  turbulent  diffusion.  This  material  balance  is  given  by  the 
following  differential  equation: 

J7(^) +  !y  (/vrf(w)dw)  -  (f^  <EX  If)  ♦  Ij-  (Ey  If)  ♦  fj  (Ez  If))  =  0  (II-l) 

where  x  =  longitudinal  coordinate  (m) 

y  =  vertical  coordinate  (m) 

z  =  lateral  coordinate  (m) 

u  =  current  velocity  (m/s) 

c  =  sediment  concentration  (mg/1) 

w  =  particle  settling  velocity  (m/s) 

f(w)  =  settling  velocity  frequency  distribution  ([sediment  mass/w]  vs. 
w) 

Ex  *  diffusion  coefficient  in  the  x-dlrection  (m^/s) 

Ey  =  diffusion  coefficient  in  the  y-direction  (m^/s) 

Ez  *  diffusion  coefficient  In  the  z-direction  (m^/s) 

The  necessary  assumptions  for  this  relationship  are  steady,  uniform,  and 
fully  turbulent  flow,  as  well  as  Ficklan  diffusion  (i.e.,  mass  transport  by 
diffusion  is  proportional  to  the  concentration  gradient). 
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To  obtain  a  useful  form  of  (II-l),  further  simplifying  assumptions  must 
be  Included: 

1).  Turbulent  diffusion  In  the  x-direction  (downstream)  is  negligible 
when  compared  to  the  convection  term,  I.e., 


3x  <Ex  3  x^  K<  3 x  (uc) 


(II-2) 


For  fully  turbulent  flow,  the  velocity  profile  will  be  practically 
flat,  especially  outside  the  near-bottom  zone.  Therefore,  it  can 
be  assumed  that  the  velocity  is  constant  and  equal  to  the  mean 
velocity,  U,  and 


(II-3) 


Ey  can  be  assumed  to  be  approximately  equal  to  the  diffusion 
coefficient  for  fluid  momentum,  Em  (the  Reynolds  analogy).  This 
assumption  has  been  found  to  be  especially  true  for  fine 
particles,  which  constitute  the  principal  portion  of  the 
suspended-solids  plume  (Jobson  and  Sayre,  1970). 

A  relationship  between  Em  and  y  can  be  obtained  using  four 
equations: 

a)  The  first  equation  Is  the  definition  of  Em, 


_  c  OV 


(II-4) 


b)  The  second  equation  is  obtained  by  differentiating  the 
logarithmic  velocity  profile 


U  ■  —  In  [y/H] 


(II -5 ) 


c)  The  vertical  distribution  of  shear  Is  assumed  to  be  given  by 


xt  •  t0  [1  -  y/H] 


(II-6) 


d)  Finally,  the  shear  velocity  and  bed  shear  are  related  as 
f ol 1 ows : 

xQ  =  pu*2  ( 1 1 -7 ) 

Combining  these  three  equations,  the  following  equation  is 
obtained: 

Em  *  <u*y[l-y/H]  (1 1-8) 

in  which 

u*  *  shear  velocity  (m/s ) 

k  =  von  Karman's  constant 

H  =  channel  depth  (m) 

Tt  =  turbulent  shear  stress  (N/m2) 

p  “fluid  density  (kg/m^) 

Em  “  momentum  diffusion  coefficient  (m2/s) 
t  0  “bed  shear  stress  (N/m2) 

Thus  It  can  be  assumed  that  ^  varies  parabolically  from  zero  at 
the  surface  and  bed  to  a  maximum  at  mid-depth. 

The  bed  shear  stress  and  mean  velocity  are  empirically 
related  as  follows: 

vi'"2  <n-9> 

where  f  is  the  friction  factor  (approximately  equal  to  0.02  for 
natural,  fully  turbulent  flows). 

If  k  is  assumed  to  be  equal  to  0.4,  the  following  relation 
can  be  obtained  by  combining  (I I -7 )  through  (I I -9) : 


Ey  =  Em  =  0.02Uy(l-y/H) 


(II -10 ) 
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4) .  Based  on  experimental  studies  (Fischer,  1973),  the  lateral 

diffusion  coefficient  can  be  considered  to  be  approximately 
constant  and  given  by  the  expression 

E,  =  0.2Hu*  (II-H) 

z 

Combining  (I 1-9)  through  (11-11),  Ez  can  be  expressed  as 

E  -  2  (E  )  (H-12) 

z  y'max  ' 

where  (Ey)max  Is  the  maximum  value  of  Ey  (at  mid-depth). 

5) .  An  assumption  Is  made  that  the  sediment  can  be  divided  Into 

fractions,  and  that  each  fraction  can  be  characterized  by  a 
representative  particle  size  with  Its  corresponding  settling 

velocity.  It  also  Is  assumed  that  no  flocculation  occurs  (1 .e. ,  w 

Is  not  a  function  of  time). 

Incorporating  all  of  these  assumptions,  the  material  balance  equation, 
(I I -1 ) ,  can  be  written  as 

u  fx*  *  fy  ■  fy  <(0-<WV(l*y/l>>)  fy) 

-lr<2<EyW  If)  - 0  (»-13> 

This  equation  Is  solved  for  each  sediment  fraction.  The  total  concentration 
at  each  point  Is  obtained  by  superimposing  the  weighted  solutions  for  all 
fractions. 

Equation  (11-13)  Is  solved  numerically  by  a  finite-difference  method, 
using  the  following  boundary  conditions: 

1).  The  upstream  boundary  condition,  at  x  ■  0,  assumes  that  a  vertical 
plane  source  continuously  emits  suspended  solids  at  a  constant 
given  strength  per  unit  height,  c  ■  c0.  The  lateral  limits  of  the 
source  are  -b<z<b,  where  b  Is  small  relative  to  the  depth.  The 
source  IS,  In  general,  assumed  to  extend  from  the  water  surface  to 


the  bottom.  Variations  In  Initial  source  conditions,  such  as  a 
point  or  a  depth-varying  source,  can  also  be  modelled. 

2) .  A  second  boundary  condition  for  the  solution  of  (11-13)  is  the 

physical  limitation  of  no  sediment  flux  across  the  surface,  i.e., 
at  y  =  H 

Ey|^+Wc  =  °  (11-14) 

3) .  For  simplicity,  it  is  assumed  that  there  is  no  reentrainment  of 

sediment  reaching  the  bottom,  i.e.,  at  y  =  0 


(11-15) 


Using  these  boundary  conditions,  the  calculation  proceeds  downstream  step 
by  step.  Over  each  step.  Ax,  a  finite-difference  model  of  the  equation  is 
used  to  calculate  the  concentrations  at  the  end  of  the  step  from  those  at  the 
beginning. 

Two  aspects  of  this  process  deserve  special  mention.  First,  in  order 
for  Ax  to  be  large  enough  to  have  an  efficient  algorithm  and  to  avoid 
Instabilities  of  a  purely  numerical  origin,  an  Implicit  finite-difference 
scheme  was  selected.  For  this  implicit  system,  the  diffusion 
term,  3/3x(Ex3c/3x),  Is  expressed  at  x  +  Ax  rather  than  at  x.  This  scheme 
requires  the  solution  of  a  system  of  N  equations  at  each  step  (N=H/Ay). 

Second,  the  effect  of  lateral  diffusion  can  be  considered  separately. 
Lateral  diffusion  could  be  described  by 


3_ 

3z 


(E  — ) 
'LZ  3z' 


(11-16) 


where  C  represents  the  concentration  obtained  from  the  two-dimensional 
problem.  The  solution  of  this  equation  is  given  by 


C(x,z) 


(11-17) 
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where  v  is  a  dummy  variable  which  represents  distance  within  the  plume.  A 
change  of  variables  relates  this  solution  to  the  “error  function"  (erf). 

If  Ez  is  taken  to  be  a  constant  (largely  true,  except  near  the  bottom 
surface),  and 

*  y2/2  (II-18) 

z 

(11-17)  becomes 

_  b  , 

C(x,z)  *  1//3T  /  exp  [-y  /2]dy  (11-19) 

-b 

which  can  be  expressed  in  terms  of  the  error  function,  for  which  an  analytical 
solution  is  well  known.  The  analytical  and  numerical  results  are  computed 
separately  and  then  combined  to  obtain  an  approximation  of  the  three- 
dimensional  concentration  field. 

Modifications  to  the  Wechsler-Cogley  computer  program  are  discussed  in 
Sections  C  through  G.  Section  F  Includes  a  sample  model  run  with  its 
corresponding  computer  output.  A  final  version  of  the  modified  program  is 
listed  in  appendix  A. 

B.  Field  Trips;  Procedure  and  Results.  To  gather  data  for  refining 
and  verifying  the  computer  model  field  trips  were  made  to  Mississippi  River 
Mile  (RM)  14,  RM  254.5,  and  RM  170.  The  sampling  process  used  at  each  site 
was  as  follows: 

1).  Site  Layout.  It  was  attempted  to  establish  a  rough  grid  at  each 
site  to  systemlze  the  subsequent  sampling.  Either  buoys  or  shore 
markers,  located  from  a  common  reference  point,  were  used  to 
demarcate  the  different  sampling  transects.  It  was  estimated  that 
this  grid  would  encompass  the  sediment  plume.  The  general  layouts 
of  the  three  sites  are  shown  In  figures  I I -1  through  I 1-7 . 
Suspended-solids  data  for  each  site  appear  on  these  figures  as 
well. 
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II- 2  rm  14:  layout  and  above-ambient  suspended 
solids  concentrations  (mg/1)  -  9-ft  depth 
(plot  2) . 


Figure  II-4  RM  254.5:  layout  and  above-ambient  suspended 
solids  concentrations  (mg/1)  -  9-ft  depth. 
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Figure  II-5  RM  254.5:  layout  and  above-ambient  suspended 
solids  concentrations  (mg/1)  -  2-ft  depth. 
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Figure  U-6  RM  170:  layout  and  above-ambient  suspended 
solids  concentrations  (mg/1)  -  9-ft  depth. 
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Figure  1 1- 7  RM  170:  layout  and  above-ambient  suspended 
solids  concentrations  (mg/1)  -  2-ft  depth. 


2) .  Water  Samples  -  Suspended  Solids.  Water  samples  were  taken  at 

depths  of  2  and  9  feet.  Sediment- laden  water  was  pumped 
continuously  as  the  boat  traversed  each  transect  at  an 
approximately  uniform  velocity.  A  small  Col  e-Darwin  rotary  pump 
was  used  to  pump  discrete  samples  into  pint  jars  at  intervals  of 
approximately  thirty  seconds.  These  samples  were  analyzed  using  a 
standard  suspended-solids  analysis  (Standard  Methods,  1980)  at  the 
Water  Plant  laboratory  of  The  University  of  Iowa.  Figures  II-l 
through  I I -7  show  the  above-ambient  suspended-solids 
concentrations  for  the  three  locations  at  2- ft  and  9- ft  depths. 

3) .  Velocity  Measurements.  To  obtain  a  representative  mean  velocity 

for  each  site,  velocity  measurements  were  taken  (with  a  Universal 
Current  Meter)  at  ten  equally  spaced  intervals  along  the  depth, 
starting  at  1  ft  from  the  bottom,  plus  one  measurement  0.4  ft  from 
the  surface.  Thus,  it  was  possible  to  construct  velocity  profiles 
for  different  locations  at  each  site.  These  profiles  are  given  in 
figures  I I -8  through  11-10.  This  procedure  also  provided 

representative  mean  depths  at  each  site. 

4) .  Bed  Samples.  To  determine  a  characteristic  dredge  material  size, 

bed  samples  were  taken  with  a  Ponar  dredge  sampler.  One  sample 
was  taken  at  some  point  which  was  to  be  dredged  ("above  dredge") 
and  another  was  taken  at  the  dump  site. 

These  samples  were  analyzed  using  the  U.S.  Bureau  of 
Standards  sieve  classification  (Lambe,  1951),  either  by  the  Iowa 
Institute  of  Hydraulic  Research  (IIHR)  or  by  the  United  States 
Geological  Survey  (USGS)  in  Iowa  City,  Iowa.  The  grain-size 
distribution  curves  appear  in  figures  11-11  through  11-13. 

5) .  Water  Samples  -  Solids  Distribution.  To  determine  the  size 

distribution  of  the  solids  In  suspension  three  water  samples, 
approximately  five  gallons  apiece,  were  gathered  at  each  site;  at 
about  20  ft  (as  close  as  possible),  40  ft,  and  100  ft  from  the 
disposal  point. 
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These  samples  were  analyzed  by  the  Vsual  Accumulation  (VA) 
Tube  and  Pipette  methods  by  the  US6S.  It  should  be  noted  that  the 
settling  velocities  of  the  larger  particles  in  suspension  can 
easily  be  determined  from  the  VA  graph. 

A  brief  description  of  each  field  trip  follows: 

1).  The  first  field  trip,  to  RH  14  (near  Cairo,  Illinois),  was  made 
between  29  September  and  1  October,  1981.  The  general  location 
map  of  the  site  appears  in  figure  11-14.  The  COE  dustpan  dredge, 
"Potter,"  was  operating  using  near-shore  disposal.  The  layout  and 
suspended-solids  plumes  appear  in  figures  I I -1  through  II-3,  with 
ambient  suspended  solids  equal  to  100  mg/1. 

The  plume  at  the  9-ft  depth  is  presented  in  two  parts 
(figures  II-l  and  1 1-2 )  because  the  dredge-disposal  point  was 
moved  approximately  500  m  upstream  before  the  final  two  transects 
were  sampled.  Due  to  a  lack  of  time  and  the  difficulties  involved 
in  reestablishing  the  grid  system,  it  was  assumed  that  the  last 
two  transects  could  be  incorporated  with  the  first  four  for 
modelling  purposes. 

According  to  the  MIT  classification  (Lambe,  1951),  the 
dredged  material  at  this  site  (see  figure  1 1 -11 )  appears  to  be 
largely  medium  and  coarse  sand  with  D5q  equal  to  about  0.5  mm. 
There  was  practically  no  silt  and  clay  present  in  the  samples. 

Figure  I I -8  illustrates  two  of  the  problems  encountered  at  RM 
14:  extremely  high  velocities  and  depths.  Although  the  figure 

indicates  a  relatively  high  velocity  of  4.5  ft/sec  and  a  depth  of 
14  ft,  velocities  of  over  8  ft/sec  "gusting"  around  the  wing  dams 
were  encountered,  as  well  as  depths  of  more  than  50  ft.  These 
velocities  and  depths  made  it  extremely  difficult  to  set  up  the 
original  grid.  The  sampling  itself  was  also  difficult  because  the 
boat  was  rather  ill-equipped  to  handle  these  conditions.  The 
presence  of  a  series  of  wing  dams  contributed  to  the  general 
unsteadiness  of  the  situation.  All  of  these  conditions  taken 
together  explain,  to  a  large  extent,  the  absence  of  well- 


Figure  11-14  Site  location  map  -  RM  14. 


established  plumes  at  this  site  (see  figures  1 1 -1  through  1 1 -3). 

2) .  The  second  field  trip,  to  RM  254.5,  was  made  between  20  October 

and  25  October,  1981.  The  general  location  map  of  the  site 
appears  in  figure  11-15.  Open-water  disposal  of  the  COE 

cutterhead  dredge,  "St.  Genevieve",  was  sampled.  The  general 
layout  and  suspended-solids  plumes  appear  in  figures  I I -4  and  II- 
5,  with  ambient  suspended  solids  equal  to  156  mg/1. 

The  plume  at  the  9- ft  depth  (figure  1 1 -4)  is  fairly  well 
established,  especially  on  transects  3,  5,  and  6  with 
concentration  values  tapering  off  at  the  edges.  The  plume  appears 
to  be  on  the  order  of  200  m  wide.  The  plume  at  the  2-ft  depth 
(figure  I 1-5)  is  less  well  established. 

Figure  I 1-9  shows  representative  velocity  profiles  and  depths 
at  this  site.  The  velocities  are  about  1.5  ft/sec  and  the  depths 
are  about  15  ft.  The  sites  where  velocities  were  measured  appear 
in  figure  I I -4  as  VI  and  V2. 

D50  for  the  bed  samples  (see  figure  11-12)  lies  between  0.3 
and  0.4  mm.  Again,  there  is  very  little  silt  and  clay;  the 
dredged  material  appears  to  be  medium  and  coarse  sand.  The 
difference  between  the  two  curves  is  probably  due  to  local  sorting 
at  the  discharge  point. 

3) .  The  third  field  trip,  to  RM  170  (near  St.  Louis,  Missouri),  took 

place  between  20  November  and  22  November,  1981.  The  general 
location  map  of  the  site  is  shown  in  figure  11-16.  The  dredge- 
disposal  plume  from  the  "St.  Genevieve"  was  sampled.  Both  near¬ 
shore  and  on-shore  disposals  were  used  at  different  times.  The 
general  layout  and  suspended-solids  plumes  for  this  site  are  given 

I 

in  figures  II-6  (ambient  suspended  solids  *  70  mg/1)  and  1 1 -7 
(ambient  suspended  solids  =  84  mg/1). 

Figure  I I -6  (9-ft  depth)  shows  a  fairly  well-established 
plume  with  the  concentrations  along  most  of  the  transects  tapering 
off  at  the  edges.  The  width  of  the  plume  is  on  the  order  of  100 
m,  but  it  is  possible  that  the  whole  plume  was  not  encompassed 


Figure  11-16  Site  location  map  -  RM  170. 


within  the  grid. 

The  2-ft  depth  plume  (figure  1 1-7 )  is  considerably  less  well 
established  than  that  at  the  9- ft  depth.  The  fourth  transect 
indicates  a  drop  in  slurry  concentration. 

Figure  I I -10  indicates  representative  velocities  of  about  2.5 
ft/sec  and  depths  ranging  from  about  10  to  20  ft.  The  measurement 
sites  are  shown  as  VI,  V2,  and  V3  in  figure  I I -6. 

The  grain-size  distribution  for  this  site  (figure  11-13) 
shows  dredge-spoil  material  consisting  of  medium  and  coarse  sand 
with  a  O50  of  about  0.4  mm,  with  little  or  no  silt  and  clay 
present.  The  difference  between  the  two  curves  Is  probably  due  to 
sorting  at  the  discharge  point.  An  empirical  relationship  between 
particle  diameter  and  settling  velocity  is  presented  in  figure  II- 
17. 

One  problem  unique  to  this  site  was  the  mode  of  disposal. 
The  disposal  started  out  as  near-shore  in  water  2  or  3  ft  deep. 
As  the  sand  settled  out  of  suspension,  it  rapidly  accumulated  and 
the  disposal  became  on-shore.  This  change  naturally  affected  the 
dynamics  of  the  whole  situation  and,  therefore,  the  formation  of 
the  plume  itself. 

C.  Model  Predictions  -  GREAT  11  Data.  The  first  modelling  objective 
of  the  present  project  is  to  attempt,  through  manipulation  of  input 
parameters,  to  improve  suspended-sol ids-concentrati on  predictions  of  field 
data  from  two  GREAT  II  sites  (Schnoor  et  al.,  1980). 

Figure  11-18  shows  field  data  from  the  Rock  Island,  Illinois  site.  The 
isopleths  superimposed  on  the  data  are  those  predicted  by  the  model  as  it  was 
in  the  GREAT  II  study,  except  the  source  half-width  has  been  Increased  from 
12.5  m  to  15  m  to  improve  the  transverse  scale  on  the  computer  output.  The 
principal  problem  is  too  little  lateral  spread. 

Lateral  spreading  can  be  increased  by  increasing  either  the  source  half- 
width  or  the  lateral  dispersion  coefficient,  Ez.  Since  increasing  the  source 
width  essentially  Increases  the  width  of  the  entire  plume  proportionately.  It 
was  decided  to  concentrate  on  increasing  the  plume  spread  by  changing  the 
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Figure  11-18  Rock  Island  suspended-solids  plume  with  pre¬ 
dicted  isopleths  superimposed  (E  =0.0088  sq. 
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lateral  dispersion  coefficient. 

The  lateral  dispersion  coefficient  for  figure  11-18  was  calculated  using 
the  assumption  that  (see  Section  A) 

Ez  =  2.2  (0.005UH)  =  0.0088  m2/s  (11-20) 

A  number  of  model  simulations  were  made  systematically  increasing  the 
lateral  dispersion  coefficient  while  keeping  all  other  parameters  the  same. 
The  most  satisfactory  prediction  was  obtained  using  Ez  =  2.0  m2/s  (figure  II- 

19) .  Because  of  its  magnitude,  this  value  of  Ez  can  no  longer  be  considered 
to  represent  only  the  effect  of  turbulent  diffusion.  For  the  remainder  of 
this  report  the  lateral  dispersion  coefficient  will  be  referred  to  as  the 
coefficient  of  transverse  mixing,  Dz.  This  coefficient  includes  effects  of 
turbulent  diffusion,  transverse  dispersion  (Taylor,  1954)  associated  with 
secondary  currents,  and  any  transverse  velocity  component  imparted  by  the 
disposal  operation  itself.  This  change  from  Ez  to  Dz  will  be  discussed  in 
more  depth  in  Section  E. 

The  same  modelling  procedure  was  followed  for  data  from  the  Keithsburg, 
Illinois,  site.  Figure  11-20  shows  the  model  results  for  parameters  used  in 
the  GREAT  II  report,  except  that  the  source  width  has  been  increased  from  6  to 
12  m  (again,  to  improve  the  computer-output  scale).  This  figure  also  shows 
the  inability  of  the  model  to  adequately  predict  lateral  spread;  near-shore 
and  in-plume  suspended-solids  concentrations  are  also  poorly  predicted.  It 
should  be  noted  that  the  lateral  dispersion  coefficient  was  taken  (in  the 
GREAT  II  study)  to  be  equal  to  0.03  m2/s,  and  It  was  not  calculated  by  (II- 

20) .  The  coefficient  was  Increased  to  include  effects  of  the  sloping  channel 
bottom. 

Although  the  Keithsburg  data  were  quite  sketchy,  Dz  =  2.0  m^/s  (figure 
11-21)  again  provided  the  most  plausible  predictions  for  a  combination  of 
lateral  spread,  near-shore  concentrations,  and  to  a  certain  extent,  the  few 
available  In-plume  concentrations. 

0.  Model  Predictions  -  GREAT  III  Pata.  Having  obtained  satisfactory 
results  for  Rock  Island  and  Keithsburg  data,  the  modified  model  was  used  to 
predict  GREAT  III  data.  The  input  parameters  for  the  RM  254.5  site  were 
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Figure  11-20  Keithsburg  suspended-solids  plume  with  predicted 
isopleths  superimposed  (D  =0.03  sq.  m/s). 
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Figure  11-21  Keithsburg  suspended-solids  plume  with  predicted 
isopleths  superimposed  (D  =2.0  sq.  m/s). 


determined  as  follows: 


1) .  The  mean  velocity  of  0.5  m/s  was  determined  from  figure  II-9. 

2) .  The  depth  of  4  m  was  determined  from  figure  I I -9. 

3) .  From  personal  observation,  the  source  width  was  assumed  to  be 

equal  to  20  m. 

4) .  Above-ambient  suspended-solids  concentrations  at  RM  254.5  (see 

figures  1 1  -4  and  1 1 -5 )  were  not  substantial  enough  to  provide  a 
basis  for  determining  the  original  source  concentration.  For  this 
reason,  a  rough  solids  balance  was  used  to  calculate  the  strength 
of  the  source. 

3 

The  "St.  Genevieve"  dredges  an  average  of  1072  yd  /hr  of 
material  (personal  communication  with  Mr.  Matthew  Struckel,  U.S. 
Army  Corps  of  Engineers,  St.  Louis).  Assuming  the  dredge  spoil  to 
be  85  percent  solids,  the  mass  of  material  discharged  from  the 

O 

dredge  is  approximately  5  x  10  mg/s.  Using  the  width,  velocity, 
and  depth  calculated  above,  approximately  12,000  mg  are  discharged 
per  liter  of  receiving  water.  Less  than  one  percent  of  this 
material  actually  becomes  suspended  in  the  river  upon  disposal. 

The  next  obstacle  encountered  was  modelling  the  low 
suspended-solids  concentrations  at  the  head  of  the  plume, 
precisely  where  the  highest  concentrations  would  normally  be 
expected.  This  phenomenon  is  discussed  in  Section  E.  In  trying  a 
number  of  different  source  configurations,  the  best  results  were 
obtained  when  the  Initial  plane  source  was  assumed  to  cover  the 
bottom  tenth  of  the  depth.  Therefore,  the  initial  source  strength 
should  be  scaled  up  to  120,000  mg/1.  Since  120,000  mg/1  is  known 
to  be  too  large,  the  source  strength  was  reduced  considerably  from 
the  original  value,  because  the  bulk  of  the  material  (mainly  sand) 
settles  out  of  suspension  almost  immediately  after  discharge. 


5).  It  was  not  possible  to  accurately  determine  the  size  division  of 


suspended-dredged  sediment  from  the  field  analyses.  In  light  of 
this  situation,  the  following  procedure  was  followed: 

a)  Since  the  data  indicated  a  fairly  constant  amount  of  sediment 
in  suspension,  the  sediment  was  divided  into  two  fractions: 
fines  and  sand. 

b)  Based  on  size  analyses,  the  grain  size  of  the  sand  in 
suspension  was  taken  to  be  about  125  microns.  From  these 
analyses,  the  settling  velocity  for  125  microns  is  approximately 
0.016  m/s.  Because  the  model  is  unstable  for  a  combination  of 
this  settling  velocity  and  the  current  velocity  at  the  site  (see 
figure  11-22),  the  settling  velocity  used  was  0.012  m/s.  Since 
even  at  the  reduced  velocity  the  sand  settles  out  quite  rapidly, 
this  reduction  has  very  little  effect  on  the  final  result. 

c)  Based  on  the  pipette  analyses  the  representative  grain  size 
for  the  fines  was  assumed  to  be  0.002  mm.  The  corresponding 
settling  velocity  determined  from  figure  11-17  Is  about  3  x  10"6 
m/s. 

d)  Based  on  size  analyses,  the  material  is  assumed  to  contain 
about  50  percent  fines  and  50  percent  sand. 

6) .  The  disposal  mode  used  was  open-water  discharge. 

7) .  Since  the  plume  spreads  out  very  quickly  at  first,  and  then  more 

gradually,  it  was  decided  to  introduce  a  crude  scheme  (this 
assumption  is  discussed  more  thoroughly  In  Section  E)  for 
modelling  the  lateral  mixing  coefficient.  It  was  assumed  that  Dz 
varies  as  follows: 

Dz(x)  =  Dz*[l+b  exp(-x/100)]  (11-21 ) 

Dz*  =  eventual  "steady  state"  lateral  mixing  coefficient 
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b  =  a  constant 

In  other  words,  the  lateral  mixing  coefficient  was  assumed  to  decrease 

exponentially  until  it  reached  an  essentially  constant  value. 

These  parameters  were  modelled  for  different  combinations  of  Dz*,  b,  and 
original  source  concentration.  The  best  predictions  were  obtained  using  Dz*  = 
2.0  m^/s  and  b  =  10  for  an  original  source  strength  of  2000  mg/1.  The 

predicted  isopleths  appear  in  figures  11-23  (9-ft  depth)  and  11-24  (2-ft 

depth).  The  lateral  spread  is  predicted  well  for  both  depths.  Concentration 
predictions  for  the  9-ft  depth  are  quite  good,  while  those  at  a  depth  of  2  ft 
are  too  low. 

The  final  phase  of  the  present  study  consisted  of  verifying  the  modified 
model  using  data  gathered  at  RM  170  (see  figures  1 1 -6  and  1 1 -7).  The 

parameters  used  for  this  verification  were  as  follows: 

1) .  Mean  velocity  =  0.8  m/s  (from  figure  1 1 -10 ) . 

2) .  Mean  depth  =  4  m  (from  figure  11-10). 

3) .  Source  width  =  20  m. 

4) .  The  source  configuration  and  original  source  strength  were  assumed 

to  be  the  same  as  those  used  at  RM  254.5. 

5) .  Even  though  there  was  a  higher  percentage  of  fines  present  in 

suspension  than  at  RM  254.5,  the  distributions  were  quite 
similar.  Because  the  grain-size  distributions  at  the  two  sites 
(figures  11-12  and  11-13)  were  also  very  similar,  the  same 

sediment  division  and  settling  velocities  as  were  used  at  RM  254.5 
were  chosen.  The  smaller  amounts  of  sand  in  suspension  at  the 
head  of  the  plume  are  probably  due  to  the  fact  that  the  dredge 
spoil  often  traveled  over  land  before  entering  the  river,  and 
therefore  the  sand  had  more  opportunity  to  settle  out.  Given  the 
similarities  in  distribution  already  mentioned  and  the  fact  that 
the  fines  predominate  In  the  plume,  the  assunption  that  the 
sediment  parameters  were  also  similar  is  justified. 

6) .  Since  the  spoil  discharge  alternated  between  near-  and  on-bank, 

the  intermediate  disposal  mode  was  used. 

7) .  Since  the  depths  and  velocities  were  of  similar  magnitudes  at  both 
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Figure  11-24  ****  254.5  (2-ft  depth)  suspended- solids  plume 
with  predicted  isopleths  superimposed  (mg/1 
above  ambient) . 


sites,  b  =  10,  and  Dz*  =  2.0  nr/s  were  used  for  determining  the 
traverse  mixing  coefficient. 

The  predicted  isopleths  using  these  parameters  are  shown  in  figures  11-25 
(9-ft  depth)  and  11-26  (2-ft  depth).  The  predicted  concentration  values  for  a 
depth  of  2  ft  are  extremely  low.  This  fact  may  be  due  to  the  failure  to 

obtain  adequate  samples  for  determining  ambient  suspended-solids 

concentrations,  which  also  explains  the  difference  in  ambient  values  used  at 
the  two  depths. 

At  a  depth  of  9  ft  the  magnitudes  of  the  predicted  values  are  much  closer 

to  those  of  the  data.  In  fact,  if  the  whole  distribution  were  rotated  away 

from  the  bank,  the  predicted  values  would  agree  quite  well  with  the  data. 
This  fact  indicates  an  abrupt  increase  in  the  depth  at  about  50  m  from  the 
bank.  This  ledge  would  prevent  spreading  and  would  virtually  constitute  a 
bank  from  which  suspended  solids  below  a  certain  depth  would  "reflect".  This 
possibility  is  borne  out  by  the  inability  of  the  boat  to  approach  closer  than 
50  m  from  the  river  bank  when  sampling  at  a  depth  o^  9  ft. 

It  should  be  noted  here  that  it  was  impossible  to  model  the  data  from  RM 
14,  due  principally  to  the  fluctuating  depths,  fluctuating  velocities,  the 
presence  of  wing  dams,  and  the  moving  dredge  disposal. 

E.  Discussion.  The  purpose  of  this  section  is  to  justify  the  major 
modifications  made  to  the  model  by  explaining  why  they  were  successful  in 
improving  the  model's  predictions. 

1.  Transverse  Mixing  Coefficient.  The  amount  of  transverse  mixing  in 
an  open  channel  can  be  represented  by  the  dimensionless 
quantity  e  =  Dz/u*H,  where  u*  =  shear  velocity  and  H  =  river 
depth.  Dz  is  the  coefficient  of  transverse  mixing  which  includes 
turbulent  diffusion  as  well  as  transverse  dispersion  (Taylor  type) 
associated  with  any  secondary  currents  present  (Taylor,  1954; 
Elder,  1959;  and  Brooks,  1972). 

The  Wechsl er-Cogl ey  model  was  based  on  the  assumption 
that  e  =  0.2,  as  reported  by  Fischer  (1973).  This  value  was 
obtained  in  flume  experiments.  Yotsukura  and  Cobb  (1972)  report 
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Figure  11-25  RM  170  (9-ft  depth)  suspended- solids  plume  with 
predicted  isopleths  superimposed  (mg/1  above 
ambient) . 
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Figure  11-26 


RM  170  (2-ft  depth)  suspended-solids  plume  with 
predicted  isopleths  superimposed  (mg/1  above 
ambient) . 


values  of  6  for  natural  streams  and  irrigation  canals  varying  from 
0.22  to  0.65. 

Okoye  (1970)  showed  that  0  is  a  function  of  the  aspect 
ratio  X  =  H/W,  where  W  is  the  river  width.  Figure  11-27  shows 
this  relationship  and  indicates  values  of  9  measured  in  the  field 
of  up  to  0.72.  If  a  depth  of  15  ft  and  a  width  of  half  a  mile  are 
assumed  at  sites  sampled  for  the  present  study,  the  aspect  ratio 
would  be  equal  to  about  0.005.  A  rough  extrapolation  in  figure 
11-27  would  indicate  a  value  of  6  equal  to  nearly  2. 

Okoye  also  reported  a  variation  with  depth  of  the  transverse 
mixing  coefficient  due  to  bed  shear  and  the  nonuniform 
distribution  of  the  vertical  diffusivity  with  depth.  A  typical 
exanple  of  this  phenomenon  is  illustrated  in  figure  11-28, 
where  c  =  x/H  and  e  =  average  value  of  0.  Okoye  reported  a 
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maximum  variation  of  35  percent  fore(n)  from  its  mean  value. 
Since  the  sampling  for  the  present  study  was  done  at  what  appears 
(from  figure  11-27)  to  be  a  depth  of  above-average  e,  a  further 
increase  of  6  by  20  percent  is  reasonable. 

Another  assumption  made  for  the  Wechsler-Cogley  model  is  a 
friction  factor,  f,  equal  to  0.02.  This  value  does  not  take  into 
consideration  the  effect  of  bed- form  drag  on  the  friction 
factor.  Using  the  Alam-Kennedy  friction  factor  method  (Vanoni, 
1977),  a  rough  estimate  of  f  =  0.1  is  obtained  for  RM  254.5.  A 
five-fold  increase  in  f  would  translate  into  an  increase  in  shear 
velocity  by  a  factor  of  2.2. 

Taking  all  of  these  factors  into  consideration,  a  transverse 
mixing  coefficient  of  between  0.6  and  0.8  n^/s  could  be  expected 
at  RM  254.5  and  RM  170.  Though  these  values  are  considerably 
lower  than  the  2.0  m^/s  used  in  obtaining  the  best  model 
predictions,  the  relative  magnitude  is  quite  reasonable. 

The  effect  of  secondary  currents  (caused  by  bottom  shear  and 
channel  irregularities)  on  transverse  mixing  is  significant. 
Yotsukura  and  Sayre  (1976)  report  values  of  9  up  to  10  for  sharp 
bends  on  the  Missouri  River.  Fischer  (1968)  reports  that  higher 
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ire  11-27  Variation  of  the  dimensionless  mixing  coefficient  with  aspect  ratio 
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values  of  e  are  also  found  near  the  banks  of  rivers.  Therefore, 
even  though  the  field  values  of  figures  11-27  already  include 
effects  of  secondary  currents  present  in  all  natural  rivers,  the 
fact  that  the  dredge  spoil  is  discharged  near  the  bank  might 
increase  0  even  further. 

Any  channel  irregularities  present  would  also  increase  e. 
Irregularities  caused  by  the  dredge  operation  itself  (in  adding 
large  amounts  of  material  to  the  river  bed)  might  be  enough  to 
cause  substantial  cross-currents.  Therefore,  when  these  two 
additional  factors  are  taken  into  account  a  lateral  mixing 

p 

coefficient  of  2.0  nr/s  would  not  be  unexpected. 

2.  Decreasing  Transverse  Mixing  Coefficient.  Although  it  would  be 
impossible  to  rigourously  justify  the  decreasing-coefficient 
scheme  used,  two  factors  do  help  to  justify  it  qualitatively. 

The  method  of  disposal  (ejection  of  the  material  against  a 
deflector  plate)  tends  to  give  the  dredged  material  a  transverse 
velocity  component  which  is  not  accounted  for  in  the  usual  mixing 
coefficient.  The  transverse  velocity  component  is  even  more 
substantial  for  bank  disposal,  where  the  discharge  runs  off  the 
bank  into  the  river. 

Various  investigators  (Brooks,  1972;  Bradsma  and  Divoky, 
1976;  and  Prych,  1970)  indicate  that  density  differences  enhance 
lateral  spreading.  Since  the  effluent  is  more  dense  than  the 
receiving  water,  it  will  sink,  inducing  secondary  currents.  This 
sinking,  of  course,  is  especially  true  near  the  discharge  point, 
where  the  density  difference  is  especially  pronounced.  Both  of 
these  factors  decrease  as  mixing  progresses.  Therefore,  the  model 
scheme  for  calculating  the  transverse  mixing  coefficient  is 
qualitatively  consistent. 

3.  Submerged  Source.  As  stated  in  section  D  the  source  was  assumed 
to  be  at  the  bottom  of  the  river  to  account  for  the  low  suspended- 
solids  concentrations  at  the  head  of  the  plume,  where  the  heaviest 
concentrations  would  normally  be  expected.  The  effluent  falls 
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very  quickly  to  the  bottom  of  the  river  and  particles  which  have 
not  settled  out  of  suspension  are  distributed  over  the  full  depth 
as  the  turbulent  diffusivity  of  the  river  takes  over. 

Bradsma  and  Divoky  (1976)  refer  to  the  first  part  of  this 
phenomenon  as  "convective  descent",  during  which  the  material 
tends  to  fall  as  a  cloud  (or  a  series  of  minute  clouds)  under  the 
influence  of  gravity.  The  slurry  acts  as  a  relatively  dense 
liquid,  and  as  such  is  affected  by  negative  buoyancy. 

This  convective  descent  is  further  enhanced  by  the  method  of 
disposal.  On  both  the  "Potter"  and  the  "St.  Genevieve",  the  spoil 
slurry,  upon  emerging  from  the  pipe,  hits  a  deflector  plate.  So 
for  much  of  the  material,  the  momentum  is  directed  downward. 
Material  which  is  deflected  horizontally  also  attains  downward 
momentum  as  it  falls  into  the  river  because  the  discharge  point  is 
ten  or  fifteen  feet  above  the  water  surface. 

F.  Cowputer  Program  Modifications  and  Sanple  Model  Run.  The  principal 
modification  of  the  original  Wechsler-Cogley  model  included  as  part  of  the 
GREAT  II  study  (Schnoor  et  al.,  1980}  was  to  convert  open-water  disposal  (one 
of  the  Wechsler-Cogley  assumptions)  to  bank  disposal.  A  "reflection 
principle"  was  introduced  to  simulate  the  effect  of  the  river  bank.  This 
principle  is  illustrated  in  figures  11-29  and  11-30. 

Figure  11-29  shows  a  typical  computer  output  of  concentration  values 
predicted  at  a  given  depth,  using  open-water  disposal.  The  edge  of  the 
original  source  is  considered  to  represent  the  bank  which  becomes  a  reflection 
axis.  All  concentration  values  above  this  "folding"  axis  are  added  to  the 
corresponding  values  below  it.  Concentrations  along  the  folding  axis  are 
doubled,  except  for  the  value  at  the  source  which  remains  the  same.  These 
concentrations  are  doubled  because  the  open-water  values  of  the  lateral 
spreading  coefficient  along  the  original  source  edge  are  averaged  with  zero, 
thus  halving  them.  The  folded  distribution  corresponding  to  figure  11-29  is 
shown  in  figure  11-30. 


Although  this  principle  was  used  in  the  GREAT  II  study.  It  was  not 
Incorporated  Into  the  computer  program.  The  first  modification  Included  as 
part  of  the  present  study,  therefore,  was  to  Include  the  folding  operation  as 
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Figure  11-30  Folded  concentration  distribution  corresponding 
to  figure  7-1  (transverse  distance  (m)  vs.  down¬ 
stream  distance  (m) ) . 


« 


51 


part  of  the  program. 

The  second  modification  made  for  the  present  study  was  to  make  available 
the  option  of  different  types  of  disposal.  The  program  now  accommodates  three 
disposal  modes: 


1 


1) .  The  open-water  disposal  mode  essentially  uses  the  original  version 

of  the  model . 

2) .  The  bank-disposal  mode  makes  use  of  the  reflection  principle,  as 

described  above. 

3) .  The  near-shore  or  intermediate  disposal  mode  also  makes  use  of  the 

reflection  principle.  The  only  difference  from  bank  disposal  is 
that  the  program  operator  is  asked  for  a  folding  parameter,  NBANK, 
which  is  the  number  of  units  of  DELZ  (lateral  distance  unit)  which 
should  be  put  between  the  bank  and  the  nearer  limit  of  the  source 
plane. 


I  The  program  provides  two  sediment-division  options:  fines/sand  and 

I  silt/clay/sand.  These  divisions  affect  output  headings  and  computer  prompts, 

h[  and  can  be  overridden  simply  by  entering  (on  prompt)  a  number  of  sediment 

fractions  (NSEDF)  larger  than  3. 

The  fourth  modification  included  in  the  final  version  of  the  program  is 
tc  provide  the  following  options  for  source  configuration: 

■ 

1).  The  first  option  is  a  rectangular  source  of  constant  concentration 
reaching  from  the  water  surface  to  the  river  bottom. 

?).  The  second  option  is  a  rectangular  source  covering  any  given 
1  portion  of  the  river  depth. 

?•).  The  third  option  is  a  triangular  source  reaching  from  the  water 
surface  to  the  river  bottom.  The  original  concentration  varies 
from  zero  at  the  water  surface  to  a  maximum  at  the  river  bottom. 

Finally,  the  transverse  mixing  coefficient  can  be  assumed  to  be  constant 
(EXPCO  »  0)  or  to  decrease  exponentially  to  an  essentially  constant  value. 
This  second  option  was  discussed  previously. 

>■« 

The  computer  program  is  written  in  FORTRAN  and  appears  in  appendix  A.  As 


listed,  it  is  for  use  on  an  interactive  terminal  and  a  PRIME  750  computer. 


Figure  11-31  shows  an  example  of  the  prompts  and  entered  parameters  which 
constitute  a  model  run;  RM  254.5  parameters  and  data  are  used  (see  Section 
D).  Figures  11-32  through  11-40  show  the  output  resulting  from  the  example  of 
figure  11-31.  The  vertical  concentration  distributions  for  fines  and  sand  are 
shown  in  figures  11-32  and  11-33;  the  total  distribution  is  shown  in  figure 
11-34.  These  three  tables  are  basicaly  the  results  of  the  numerical  solution 
discussed  in  Section  A. 

The  results  of  the  analytical  solution,  the  lateral  spreading 
coefficients,  are  shown  in  figure  11-35.  It  should  be  noted  that  the  values 
appearing  in  figure  11-35  correspnd  to  the  decaying  mixing  coefficient 
discussed  in  Section  D. 

The  final  concentration  predictions  are  given  in  figures  11-36  through 
11-40.  Essentially,  these  tables  constitute  horizontal  "slices"  of  the 
suspended-sediment  plume  taken  at  five  different  depths,  from  the  water 
surface  (0.00  m)  to  a  plane  0.80  m  above  the  bottom.  These  concentration 
values  are  obtained  by  successively  multiplying  the  original  concentration 
(2000  mg/1,  in  this  case)  by  the  corresponding  coefficients  from  figures  11-34 
and  11-35.  It  should  be  remembered  that  the  coefficients  from  figure  11-34 
are  percentages;  they  should  be  divided  by  100  before  using  them  to  obtain  the 
final  concentrations. 

The  graphical  representation  of  figure  11-37  appears  in  figure  11-24,  and 
that  of  figure  11-39  in  figure  11-23.  It  should  be  noted  that  the  depths  of  2 
ft  and  9  ft  are  rounded  off  to  0.80  and  2.40  m  respectively. 

6.  Model  Sensitivity  Analysis.  The  objectives  of  this  section  are  (1) 
to  give  an  idea  of  how  the  model  behaves  by  performing  a  sensitivity  analysis, 
and  (2)  to  provide  a  method  of  quickly  and  easily  estimating  the  spread  and 
concentrations  of  a  dredge-disposal  plume. 

1).  Instability.  In  its  present  form,  the  model  is  numerically 
unstable  for  a  combination  of  low  river  velocities  and  high 
settling  velocities.  Since  only  high  settling  velocities  are 
affected,  this  problem  is  of  no  consequence  for  finer  material, 
which  constitutes  most  of  the  above-ambient  suspended  load.  The 
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Figure  11-31  Sample  interactive  session 
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Figure  11-32  Vertical  distribution  of  sand  (depth  (m)  vs.  downstream  distance  (m) ) . 
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Figure  11-34  Total  vertical  distribution  (depth  (m)  vs.  downstream  distance  (m) ) . 
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Figure  11-37  Suspended -solids  distribution  (mg/1)  in  the  horizontal  plane  at  a  depth 
of  0.80  m  (transverse  distance  (m)  vs.  downstream  distance  (m) ) . 
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Tire  11-40  Suspended-solids  distribution  (mg/1)  in  the  horizontal  plane  at 
of  3.20  m  (transverse  distance  (m)  vs.  downstream  distance  (m) ) 
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area  of  instability  is  depicted  in  figure  11-22;  instability  for  a 
source  covering  the  bottom  tenth  of  the  depth  Is  somewhat  more 
extensive  than  that  for  the  full-source  configuration. 

2) .  Vertical  and  Longitudinal  Concentration  Distribution.  The 

numerical  part  of  the  program  solution  (see  section  A,  as  well  as 
figures  11-32  and  11-33)  can  be  represented  by  the  following 

equation: 

P(y,x)  =  P*  (y)  expt-mX*]  (11-22) 

in  which 

P(y,x)  *  percentage  of  the  original  concentration  of  a  certain  grain 
size  still  in  suspension 
P'(y)  =  a  boundary  condition  percentage 

X*  =  dimensionless  distance  =  x/100  m 

m  =  a  constant. 

The  slope  m,  which  is  a  function  of  mean  depth  and  U/W  (mean  river 
velocity/particle  settling  velocity),  can  be  determined  from 
figure  11-41.  P ' (y )  is  calculated  using  figure  11-42.  Although 
(11-22)  does  not  strictly  hold  for  the  first  50  or  100  m,  the 

values  calculated  in  this  range  using  this  method  are  generally 

within  5  percent  of  the  model -generated  values. 

Figures  11-41  and  11-42  provide  a  good  method  of  estimating 
how  much  of  a  certain  size  of  sediment  remains  in  suspension 
downstream  from  the  dredge  disposal.  Figures  11-41  and  11-42  also 
give  a  good  estimate  of  concentrations  for  the  case  of  a 
fractional  source,  especially  after  the  first  200  m.  Care  should 
be  taken  in  this  case  to  scale  down  the  percentages  calculated 
because  a  less-extensive  source  would  be  used. 

3) .  Lateral  Plume  Spread.  As  calculated  by  the  model,  the  lateral 

spread  of  the  dredge-disposal  plume  is  a  function  only  of  the 
transverse  mixing  coefficient.  The  amount  of  lateral  spread  can 


LOGARITHMIC  SLOPE 


Figure  11-41  Determination  of  logarithmic  slope,  m 
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be  represented  by  the  following  expression: 

s*(x)  =  a(ozyx^ 


(11-23) 


in  which 

S*  =  dimensionless  spread  =  S/500  m. 
a(Dz)  =  spread  coefficient  (a  constant) 

X*  =  dimensionless  downstream  distance  =  x/500  m 
S  =  spreading  distance 

Equation  (11-23)  can  be  used  to  estimate  how  far  the  effect  of  the 
dredge-disposal  operation  will  extend.  The  spread  coefficient, 
a(Dz),  is  determined  from  figure  11-43.  S  is  the  distance  which 
the  plume  will  spread  (as  defined  by  a  lateral  spreading 
coefficient  =  0.001)  beyond  the  original  source  width  on  one  side 
of  the  plume.  Therefore,  if  open-water  disposal  is  used,  S  should 
be  doubled. 

4).  As  mentioned  in  Sections  0  and  F,  the  model  provides  the  option  of 
a  decreasing  rate  of  lateral  spread,  which  can  be  represented  as 
f ol 1 ows : 

Dz(x)  =  Dz*(l+b  expC-X*])  (11-24) 

in  which 

X*  =  dimensionless  downstream  distance  =  x/100  m 

Dz*  =  "steady  state"  value  of  the  transverse  mixing  coefficient 

b  =  a  constant 

Figure  11-44  shows  the  variation  of  the  factor  (1+b  exp[-X*])  with 
longitudinal  distance.  As  can  be  seen  from  the  figure,  this 
particular  spreading  scheme  has  little  effect  past  500  m. 

H.  Su—ary.  The  Wechsler-Cogley  (1977)  model  code  has  been  modified 
to  include  three  source  options  and  two  options  on  the  lateral  (transverse) 
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mixing  coefficient.  These  options  were  consistent  with  observations  of 
dredging  disposal  operations  in  the  GREAT  II  and  III  reaches  of  the  upper 
Mississippi  River.  The  model  was  applied  to  two  sites  from  GREAT  II  and  three 
sites  in  GREAT  III.  Summary  plots  of  how  to  estimate  model  coefficients  and 
suspended-solids  concentrations  were  given  in  figures  11-27,  11-28,  11-41,  II- 
42,  11-43,  and  equation  11-22. 

Suspended-solids  concentrations  were  quite  small  in  the  GREAT  III  reach 
which  prevented  further  refinement  and  verification  of  the  model.  Better 
monitoring  of  background  concentrations  in  the  field  would  have  been 
helpful.  Furthermore  the  assumption  of  steady  state  did  not  always  hold  true 
in  the  field  due  to  the  variable  dredging  operation  and  river  dynamics, 
lateral  spreading  of  the  plume  was  greater  than  expected  due  to  lateral 
velocities  of  the  discharge  and  channel  morphology. 
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III.  STANDARD  ELUTRIATE  TEST 


A  number  of  hazardous  chemicals  are  often  associated  with  suspended 
and  sediment  particles,  and  as  a  result  the  sediments  are  generally  considered 
to  be  sinks  for  a  variety  of  chemical  pollutants.  These  chemical  pollutants 
may  be  released  when  the  sediments  are  resuspended  by  dredging  and  dredge- 
asipOi,al  operations.  The  standard  elutriate  test  attempts  to  simulate  the 
open-water  disposal  of  dredge  material  in  the  laboratory.  The  purpose  of  this 
test  is  to  predict  the  extent  of  release  of  chemical  pollutants  from  the 
sediments  to  be  dredged. 

A,  Description  of  Sanpling  Sites  and  Sampling  Procedures.  Five 
sites  were  sampled  along  the  GREAT  III  reach  of  the  Missippi  River  (figures 
III-l,  1 1 1-2,  1 1 1-3,  and  1 1 1-4).  The  first  attempt  to  collect  water  and 

sediment  samples  for  the  elutriate  test  was  made  on  14  August  1981.  Sampling 
was  done  from  an  18-foot  flatbottom  Jon  boat  with  a  55  horsepower  Mercury 
motor.  The  boat  was  equipped  with  a  winch  capable  of  handling  over  100 

pounds.  The  boat  and  motor  were  supplied  by  the  Iowa  Institute  of  Hydraulic 

Research  of  the  University  of  Iowa.  In  the  first  sampling  attempt,  samples 
were  collected  from  three  of  five  recommended  sites  at  river  mile  (RM)  170, 
172  and  177.  Sediment  samples  were  collected  from  three  different  locations 
at  each  site  and  labeled  as  -1,  -2  and  -3  from  the  northern  most  location. 
Sediment  samples  were  collected  with  a  ponar  dredge  in  11  to  17  feet  of  water, 
at  100  to  300  feet  from  the  shoreline.  The  sediment  was  immediately  mixed  in 
a  plastic  tub,  and  a  homogeneous  portion  was  stored  in  glass  jars.  Water 
samples  were  collected  in  5-gallon  Nalgene  containers  from  the  surface  at  the 
northern  n.uit  location  of  each  site.  All  containers  had  previously  been  acid 
washed  and  rinsed  with  deionized-distilled  water.  The  samples  were  cooled 
durinc,  transit  and  stored  in  the  laboratory  at  4°C. 

Site  KM  170  is  located  near  the  National  Lead  Industries  on  the 

Illinois  shore.  Sample  RM  170-3  was  collected  at  the  southern-most  wing  dam 
which  was  located  by  shore  debris.  Two  other  dams  were  not  apparent. 
Sediment  samples  were  taken  at  approximately  700-foot  intervals  upstream. 

Site  RM  172  is  located  below  the  Des  Peres  River  on  the  Missouri 
shore.  The  first  sediment  sample  at  this  site  was  collected  approximately  100 


Figure  1 1 1-2.  Sampling  Site  at  River  Mile  177  near  St.  Louis,  Mo 


Figure  III-3.  Sampling  Site  at  River  Mile  218  near  Grafton,  IL. 
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feet  out  from  the  shoreline  of  the  southern  bank  of  the  Des  Peres  River.  Two 
other  samples  were  taken  at  approximately  600  foot  intervals  downstream  from 
the  barge  loading  areas. 

Site  RM  177  is  located  below  the  Monsanto  terminal  on  the  Illinois 
shore.  The  first  sediment  sample  at  this  site  was  collected  at  RM  178  on  the 
eastern  bank  (Illinois  shore).  Two  other  samples  were  taken  at  approximately 
600  foot  intervals  downstream  from  the  barge  loading  areas. 

A  second  sampling  attempt  was  made  on  26  August  1981.  Sediment  and 
water  samples  were  collected  from  the  remaining  two  of  five  designated 
sampling  sites.  The  sites  were  sampled  at  RM  218  on  the  Illinois  shore, 
corresponding  to  the  Illinois  River  miles  0  to  1,  and  at  RM  230  near  the 
Missouri  shore.  As  in  the  previous  sampling  event,  sediment  samples  were 
collected  form  three  different  locations  at  each  site,  and  water  samples  for 
use  in  the  elutriate  test  were  collected  from  three  different  locations  at 
each  site,  and  water  samples  for  use  in  the  elutriate  test  were  collected  from 
the  surface  at  the  northern  most  location  of  each  site.  All  sediment  samples 
were  uniformly  mixed  before  distribution  to  the  glass  jars.  Sediment  samples 
for  the  priority  pollutants  analysis  were  taken  at  the  northern  most  location 
(-1). 

Site  RM  218  is  located  near  Venetian  Harbor,  approximately  Illinois 
river  mile  1.  The  first  sample  (RM  218-1)  actually  was  collected  from  the 
Illinois  River  150  feet  from  shore.  The  water  depth  was  only  4  to  5  feet. 
The  water  was  slow  moving  and  ap##le  0.7,  approximately  450  feet  offshore  from 
Grafton,  Illinois.  The  water  depth  was  about  30  feet.  The  sediment  was 
composed  of  fine  sand  and  silt.  Sample  RM  218-3  was  obtained  at  150  feet  from 
the  Illinois  shore  in  22  feet  of  water  depth,  downstream  from  the  Riverway 
Shipping  company  commercial  dock  at  the  confluence  of  the  Mississippi  and 
Illinois  Rivers  on  the  Illinois  shore.  The  dissolved  oxygen  level  was  8.0 
ppm. 

Site  RM  230  is  located  near  Kirmill  Wood  Boat  club  on  the  Missouri 
shore.  All  three  locations  at  river  mile  230  were  sampled  at  125  feet  from 
the  Missouri  shore  In  23  to  25  feet  of  water.  The  three  samples  were 
collected  at  approximately  100  foot  intervals  near  boat  ramps.  All  sediment 
samples  were  composed  of  coarse  sand. 
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An  additional  sampling  atempt  was  made  on  November  20,  1981.  Water 
and  sediment  samples  were  taken  at  RM  170,  several  hundred  feet  upstream  from 
a  series  of  wing  dams  and  a  barge  docking  area  on  the  Illinois  shore.  Sample 
RM  170*- 1  was  taken  from  a  dredge  cut  area,  100  feet  upstream  from  the  dredge 
in  the  open  channel.  Samples  designated  by  *  represent  repetitions  of  samples 
at  the  same  site.  Sample  RM  170*-2  was  collected  from  a  disposal  area,  50 
feet  downstream  from  the  dredge  in  the  open  channel.  These  samples  were 
collected  approximately  100  feet  west  of  the  dredge  and  well  out  of  the  plume 
area. 


B,  Procedure  of  the  Elutriate  Preparation.  The  procedure  used  to 
prepare  elutriates  is  outlined  in  the  Federal  Register  (1975)  and  was  detailed 
by  the  Environmental  Effects  Laboratory,  WES  (1976).  The  following  is  a 
summary  of  the  elutriate  test  procedure  used  during  the  study. 

Step  1:  One  liter  of  wet  sediment  was  placed  in  a  5-gallon  Nalgene 
container.  Four  liters  of  unfiltered  dredge  site  water  were  added  to  the 
container  to  give  a  final  sediment  to  water  volumetric  ratio  of  1:4.  This 
mixture  was  then  referred  to  as  the  elutriate  slurry. 

Step  2:  A  large  coarse  bubble  stone  diffuser  was  placed  at  the 
bottom  of  the  container.  Compressed  air  was  passed  through  a  deionized  water 
trap  before  going  to  the  diffuser.  The  air  flow  rate  was  sufficiently  high  to 
assure  good  mixing  and  oxic  conditions. 

Step  3:  Every  ten  minutes  the  tanks  were  mechanically  mixed  by  hand 
(vigorous  swirling)  to  assure  complete  mixing. 

Step  4:  After  30  minutes  the  air  flow  was  shut  off  and  the 
elutriate  slurry  was  allowed  to  settle  for  one  hour. 

Step  5:  After  settling,  the  supernatant  was  decanted  and 
centrifuged  to  reduce  the  suspended-solids  concentration  and  to  reduce  the 
time  needed  for  subsequent  filtering. 

Step  6:  After  5  minutes  of  centrifugation  the  sample  centrate  was 
pressure-filtered  through  an  acid-rinsed  0.45  micron  filter  using  an  all  glass 
filter  apparatus. 

Step  7:  the  filtrate  is  the  "standard  elutriate".  It  was  stored  in 
a  plastic  container  at  4°C  in  the  dark.  All  analyses  were  completed  within 
two  weeks  of  sample  collection. 
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C.  Analytical  Procedure. 

1.  Heavy  Metals  Analyses 

a.  Digestion  Procedure.  For  the  determination  of  total  metals 
in  raw  river  water  the  sample  was  acidified  with  1:1  HNO3  to 
a  pH  of  2  at  the  time  of  collection.  The  sample  was  not 
filtered  before  processing.  A  representative  aliquot  of  100 
ml  of  the  well -mixed  sample  was  transferred  to  a  250  ml 
beaker  and  3  ml  of  concentrated  HNO3  was  added.  The  beaker 
was  covered  with  a  watch  glass  and  returned  to  the  hot 
plate.  The  temperature  of  the  hot  plate  was  increased  so 
that  a  gentle  reflux  action  occurred.  Heating  conti nued, 
with  additional  acid  added  as  necessary,  until  the  digestion 
was  complete  (generally  indicated  by  a  light  colored 
residue).  Sufficient  1:1  HCL  was  added  and  the  beaker  warmed 
again  to  dissolve  the  residue.  Beaker  walls  and  watch  glass 
were  washed  down  with  distilled  water  and  the  sample  filtered 
to  remove  silicates  and  other  insoluble  material  that  could 
clog  the  atomizer.  The  volume  was  adjusted  to  approximately 
50  ml.  Filtered  and  nonfiltered  standard  elutriate  samples 
were  digested  in  the  same  manner.  Sediment  samples  were 
dried  overnight  at  105°C  and  weighed.  The  dried  sample  was 
ground  with  a  mortar  and  pestle,  and  a  representative  aliquot 
of  1  gram  was  mixed  with  approximately  50  ml  of  distilled 
water.  This  sediment  slurry  was  digested  in  the  same  manner 
described  above. 

b.  Extraction  Procedure.  Digested  samples  were  extracted  prior 
to  analysis.  The  pH  of  the  samples  was  adjusted  to  pH  2.5 
using  30  percent  NH40H  and  2.5  percent  HC1.  The  samples  were 
transferred  to  a  100  ml  volumetric  flask,  and  2.5  ml  fresh 
APDC  (ammonium  pyrrolidine  dithiocarbonate)  solution  was 
added.  A  10.0  ml  MIBK  (methyl  isobutyl  ketone)  was  added, 
and  the  samples  were  shaken  vigorously  for  one  minute.  The 
layers  were  allowed  to  separate,  and  deionized  distilled 
water  was  added  until  the  MIBK  layer  was  completely  in  the 


neck  of  the  flask.  Sample  blanks,  standard  blanks,  and 

prepared  standards  were  subjected  to  identical  extraction 
procedures.  The  MIBK  layer  was  sampled  for  analysis  via 
atomic  absorption  spectroscopy. 

c.  Atomic  Absorption.  Extracted  samples  were  analyzed  on  a 
Jarrell-Ash  Atomic  Absorption  Unit  model  No.  280  using  a 
Fisher  Micro  Thermal  Atomizer  (Model  2).  A  20  microliter 

sample  was  transferred  to  the  tantulum  ribbon  of  the  MTA. 
Maximization  of  dry,  ash,  and  atomize  cycles  was 

accomplished.  Following  wavelength  maximization  from 
suggested  values,  20  microliter  samples  were  analyzed  until 
reproducibility  of  0.1  unit  was  achieved  on  at  least  two  runs 
of  the  sample.  Standards  and  blanks  also  were  analyzed  in 

this  same  manner. 

Traditional  Pollutants  Analyses.  The  analytical  procedures  for 
ammonia,  phosphorus,  chemical  oxygen  demand,  and  oil  and  grease 
followed  "Standard  Methods  for  the  Examination  of  Water  and 
Wastewater"  (1980). 

Priority  Pollutants  Analyses.  The  concentrations  of  toxic 
chemicals  known  as  priority  pollutants  were  determined  in  accor¬ 
dance  with  Federal  Register  3  December  1979. 

Other  Measured  Parameters.  Dissolved  oxygen  measurements  of  raw 
river  water  and  sediment-water  interface  water  were  performed  in 
the  field  with  a  YSI  Oxygen  Meter  Model  No.  57.  Oxidation-reduc¬ 
tion  potential  measurements  were  accompl ished  on  raw  river  water 
and  sediment  samples  in  the  field  using  an  Orion  Research  Model 
399A/F  analog  pH  meter.  In  the  laboratory,  pH  was  determined 
using  an  Orion  research  grade  pH  electrode.  Redox  potential 
measurements  were  accomplished  using  a  combination  platinum 
thimble  type  electrode.  Redox  potential  and  pH  measurements  of 
the  elutriate  samples  were  made  using  the  same  analog  pH  meter. 
Total  suspended  solids  and  volatile  suspended  solids  concentra¬ 
tions  were  determined  in  accordance  with  Standard  Methods  (1980) 
at  the  water  plant  laboratory  of  the  University  of  Iowa.  Other 


chemical  analyses  were  performed  at  the  University  Hygienic 
Laboratory  of  the  University  of  Iowa. 

D.  Results  and  Discussion.  Chemical  analyses  of  site  water, 

sediments  and  elutriates  were  performed  for  five  sites  and  fifteen  samples  in 
the  GREAT  III  study  reach.  The  concentrations  of  traditional  pollutants, 
heavy  metals,  and  PCBs  (polychlorinated  biphenyls)  in  unfiltered  site  water 
and  sediments  are  given  in  table  I II— 1 .  Concentrations  of  pollutants  in 
filtered  elutriates  as  compared  with  unfiltered  site  water,  are  presented  in 
table  1 1 1 -2.  The  heavy  metal  parameters  included  cadmium  (Cd),  chromium  (Cr), 

copper  (Cu),  lead  (Pb),  mercury  (Hg),  nickel  (Ni),  manganese  (Mn ) ,  iron  (Fe), 

and  zinc  (Zn).  Associated  statistical  analyses  were  performed.  Sediment 
samples  were  also  tested  for  the  129  priority  pollutants  (U.S.EPA,  1979). 

1.  Comparison  of  Elutriate  with  Site  Water.  Results  of  a 

statistical  analysis  on  the  water  quality  data  are  given  in  table 
III -2 .  T  is  t-statistic  for  the  null  hypothesis:  u^u^O,  Uj 
and  U2  are  the  mean  concentrations  of  a  pollutant  in  the 
unfiltered  site  water  and  the  filtered  elutriate,  respectively. 
Manganese,  ammonia,  and  COD  were  consistently  desorbed  (at  the 
0.05  significance  level  in  a  paired  t-test).  Zinc,  iron  and 
phosphate  showed  strong  tendencies  for  adsorption.  The  results 
are  in  good  agreement  with  those  of  the  GREAT  II  study  (Schnoor 
et  al.,  1979) 

According  to  chemical  equilibrium  theory,  manganese 
+2 

exists  predominantly  as  Mn  in  the  reduced  environment  of  the 

+2 

sediments.  Intersticial  waters  (rich  in  Mn  )  of  the  sediment 

are  mixed  with  river  water  in  the  elutriate  test,  and 

consequently  the  filtered  elutriate  water  contains  increased 

+2 

dissolved  manganese.  Precipitation  of  Mn  to  Mn0£  or  MN00H  is 
not  favored  thermodynamically,  or  else  kinetics  limit  the  nuclea- 
tion,  crystal -growth,  and  sedimentation  processes.  Ammonia  and 
COD  are  rich  in  sediment  intersticial  waters  and  tend  to  increase 
in  the  standard  elutriate  test.  Iron  exists  in  reduced  sediments 
as  Fe+2,  but  is  rapidly  precipitated  and  adsorbed  as  Fe+3,  when 
the  electrode  potential  increases  to  >  +  0.2mv  in  oxygenated 


Table  III-l.  Concentration  of  Heavy  Metals,  Traditional  Pollutants  and  PCBs 
In  Unfiltered  River  Water  (W)  and  Sediment  (S) 


S  =  Sediment 

N.O.  3  Not  detected 

—  »  Chemical  analysis  not  performed 


Table  III —1 .  (Continued) 
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Table  111*1.  (Continued) 


Restapled  on  20  November  1981 


Concentration  of  Pollutants  In  Filtered  Elutriates  and  Unflltered  Site  Water 


— chaaical  analysis  not  porforawd 

Data  on  Samples  170*-1  and  170*-2  Mart  not  usad  for  tha  statistical  analysis 


elutriate  water.  Zinc  and  phosphate  either  adsorb  to  particulate  matter  or 
precipitate  or  both.  Phosphate  Is  known  to  adsorb  to  silts  and  clays,  but  It 
could  also  precipitate  with  Iron  (III)  or  form  Ca- hydroxyapatite. 

The  ammonia  desorbed  during  the  elutriate  test  occasionally  exceeded 
the  water  quality  criterion  of  0.02  mg/1  un-ionlzed  ammonia.  The  heavy  metals 
concentrations  at  the  study  sites  were  below  the  criteria  for  freshwater 
aquatic  life  (U.S.EPA,  1976).  These  low  heavy  metal  concentrations  are  most 
likely  due  to  association  with  suspended  solids,  promoted  by  the  relatively 
high  pH  of  the  river  water. 

2.  Correlations  between  Traditional  Pollutants.  The  statistical 
analysis  used  here  was  a  simple  linear  regression  which  gives  a 
slope.  Intercept,  and  correlation  coefficient.  A  high 

correlation  coefficient  In  association  with  a  significant  slope 
was  considered  to  Indicate  the  possibility  of  association  between 
two  parameters. 

a.  Site  Water.  The  concentration  of  total  suspended  solids 
(TSS) ,  volatile  suspended  solids  (VSS,  ammonia-nitrogen  (NH3- 
N),  ortho-phosphorus  (PO4-P),  oil  and  grease,  chemical  oxygen 
demand  (COD),  and  redox  potential  are  given  In  table  III-3, 
and  the  associated  statistical  analysis  Is  In  table  III-4. 
Unflltered  COD  showed  an  Increase  with  Increasing  TSS  (figure 
III-5),  and  with  Increasing  VSS  (figure  II 1-6) .  The  high 
correlation  coefficients  and  significant  slopes  Indicate  a 
strong  tendency  for  unflltered  COD  to  Increase  as  TSS  or  VSS 
Increases.  Statistically  significant  correlations  were  not 
found  for  other  parameters. 

b.  Sediment.  Bed  sediment  consists  primarily  of  sand,  silt, 
clay,  and  organic  materials.  Generally  the  "dirty"  sediment 
samples  are  darker  In  color  and  contained  finer-grained 
sediment.  These  sediments  are  high  In  organic  matter  and  are 
more  likely  to  adsorb  pollutants  and  to  create  a  reduced 
(lower  redox  potential)  environment. 


Traditional  Pollutant  Concentrations 


COD  vs.  VSS 

MISSISSIPPI  RIVER  WATER 
r  «  0.93 
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The  concentration  of  ammonia,  oil  and  grease,  COO, 
redox  potential,  and  moisture  content  (percent)  of  the 
sediment  samples  are  given  in  table  I I 1-5 ,  and  the  associated 
statistical  analysis  in  table  I I 1-6.  COD  in  the  bulk 
sediment  showed  an  Increasing  tendency  with  increases  in 
moisture  content  of  the  sediment  (figure  1 1 1 -7 ) ,  and  a 
decreasing  tendency  with  Increases  in  redox  potential  (figure 
III-8).  The  higher  the  moisture  content  generally  the  more 
"organic"  the  sediment.  The  high  correlation  coefficients 
and  significant  slope  indicate  strong  association  among  COD, 
redox  potential,  and  moisture  (and  organic)  content. 

Note  the  that  organic  sediments  were  obtained  at 
sites  RM  170-2,  172-1,  172-3,  218-1,  218-2,  and  218-3.  These 
sediments  were  characterized  by  high  moisture  content,  COD, 
oil  and  grease,  ammonia,  and  phosphorus.  Sample  RM  170-2  was 
collected  behind  the  wing  dam  where  water  was  moving 
slowly.  The  site  RM  172  was  located  below  the  Des  Peres 
River  in  an  area  of  stagnant  water.  Site  RM  177-3  was 
located  below  the  Monsanto  Terminal  on  the  Illinois  shore,  in 
the  barge  loading  area.  A  large  number  of  rocks  and  twigs 
were  sampled  with  the  sediment.  Site  RM  218  was  located  near 
Venetian  Harbor.  The  water  was  slow  moving.  Sample  RM  218-3 
was  collected  downstream  from  the  Riverway  Shipping  Company 
commercial  dock  at  the  confluence  of  the  Mississippi  and 
Illinois  Rivers  on  the  Illinois  shore.  The  sediment  was 
composed  of  fine  sand  and  silt. 

The  organic  sediments  sampled  from  various 
locations  consisted  of  large  portion  of  silt,  clay-size 
particles,  and  colloidal  organic  particles.  The  silt  and 
clay-size  particles  are  transported  mainly  In  suspension,  and 
suspended  colloidal  particles  behave  much  like  the  ambient 
fluid.  Once  these  particles  reach  slow  moving  water,  the 
lower  current  velocities  allow  these  particles  to  settle. 
These  particles  are  eventually  deposited  on  the  bed  where 
they  behave  as  bed-material  particles  until  they  are  re- 
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Table  II 1-6.  Correlations  between  Sediment 

Pollutants, 

Moisture  Content, 

and 

Redox  Potential 

Correlation 

Parameters 

Intercept 

Slope 

Coefficient 

NH3-N  (mg/kg)  vs. 

“  0.772 

aaunnni 

0.66 

Redox  (mV) 

57.511 

-0.232 

0.56 

Moisture  (%) 

75.405 

5.428 

0.71 

Oil  &  Grease  (mg/kg) 

48.473 

0.072 

0.17 

COO  (mg/kg)  vs. 

Redox  (mV) 

16697.3 

-82.744 

0.94 

Moisture  {%) 

-20745.6 

1513.16 

0.93 

Oil  &  Grease(mg/kg) 

12756.6 

38.613 

0.43 

Redox  (mV)  vs. 

Moisture  (%) 

411.987 

-16.567 

0.90 

Oil  & 

Grease(mg/kg)  vs. 

Redox  (mV) 

62.985 

-0.340 

0.34 

Data  from  Sample 

170*-1  and  170*-2  were  not 

used  for  the 

statistical  analysis 

COD  vs.  moisture:  content 

MISSISSIPPI  RIVER  SEDIMENT 
r  =  0.93 


MOISTURE  [%] 


Figure  III-7.  Sediment  Chemical  Oxygen  Demand  vs.  Moisture  Content 


entrained  in  the  flow.  The  sedimentation  of  the  suspended 
colloidal  particles  involves  the  process  of  coagulation  which 
is  strongly  a  function  of  pH,  temperature,  and  other 
physical-chemical  factors.  Since  the  sedimentation 
characteristics  are  constantly  changing,  it  is  difficult  to 
predict  the  fate  of  these  particles.  However,  it  is  clear 
that  the  organic  sediment  will  be  found  in  slow-moving  or 
stagnant  water  areas. 

c.  Elutriate.  The  concentration  of  ammonia  and  COD  of  the 
elutriate  samples  and  the  statistical  analyses  are  given  in 
tables  1 1 1 -7  and  III -8 ,  respectively.  No  significant 
correlation  was  found  between  ammonia,  phosphorus,  and  COD. 

3.  Priority  Pollutants  in  Sediments.  The  EPA  has  been  directed  to 
regulate  discharges  of  129  priority  toxic  compounds  into  the 
nation's  waterways  to  protect  both  aquatic  and  human  life 
(Federal  Register,  Dec.,  1979).  Sediment  characterization  by  a 
"gas  chromatograph-mass  spectrometer"  scan  was  performed  for  the 
organic  parameters  at  the  University  Hygenic  Laboratory.  Table 
II 1-9  shows  the  concentration  of  the  organic  priority  pollutants 
at  the  five  sampling  sites.  All  pollutants  were  below  the 
detectable  level,  except  diethyl  phthalate  and  di-n-butyl 
phthalate  which  were  found  to  be  25.3  ug/1  and  11.1  yg/1,  respec¬ 
tively,  at  site  172.  Detection  limits  are  the  limits  given  for 
each  compound  as  less  than  (<). 

E,  Summary.  A  summary  of  all  desorbed  pollutants  at  each  site  is 
presented  in  table  III-10.  Manganese,  ammonia,  and  COD  were  consistantly 
desorbed  while  zinc,  iron,  and  phosphate  showed  strong  tendencies  for 
adsorption.  The  "dirty"sediments  were  found  in  slow-moving  or  stagnant  water 
areas,  and  were  characterized  by  high  organic  and  moisture  contents,  high  COD, 
oil  and  grease,  ammonia,  and  phosphorus,  and  low  redox  potentials. 


Table  I I I -8 .  Correlations  between  Traditional  Pollutants  In  the  Elutriate 


Intercept 

Slope 

Correlation 

Coefficient 

vs. 

NH3-N  (mg/1) 

vs. 

PO4-P  (mg/1) 

0.042 

8.343 

0.51 

P04-P  (mg/1) 

vs. 

COD  (mg/1) 

53.590 

7.850 

0.62 

Data  from  Sample  170*-1  and  170*-2  were  not  used  for  the  statistical  analysl 


Table  III-9.  Concentrations  of  the  Organic  Priority  Pollutants  In  Sediments 


dl benzo (a ,h )anthracene 

1 .2- dlchl orobenzene. . . , 

1.3- dl chlorobenzene..., 


(Continued) 


Table  111-10.  Desorbed  Chemicals  In  the  Standard  Elutriate  Test 
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IV.  DESORPTION  KINETICS 

The  kinetics  study  attempts  to  explain  the  mechanisms  of  sorption  during 
the  process  of  elutriatlon.  While  the  standard  elutriate  test  is  a  physical 
model  to  simulate  the  effect  of  dredge  disposal  on  water  quality  in  the 
laboratory,  the  kinetic  model  presented  in  this  report  is  a  mathematical  model 
to  describe  the  sorption  rates  for  various  chemicals. 


A.  Description  of 


>11  ng  Sites  and 


>11na  Procedures.  Based  on  the 


results  obtained  from  previous  elutriate  tests,  it  was  decided  that  the 
chemical  pollutants  of  interest  were  ammonia,  chemical  oxygen  demand, 
chromium,  and  nickel.  Sampling  sites  selected  were  as  follows: 


170-2  Silty  sediment  behind  old  wing  dam  on  Illinois  shore. 

172-1  Silty  sediment  at  the  mouth  of  the  Des  Peres  River  on  the  Missouri 
shore. 

172-2  Silt  and  sand  mixture  approximately  150  meters  downstream  from 
172-1. 


The  first  sampling  trip  was  made  on  September  22,  1981.  The  water  level 
was  approximately  10  feet  lower  than  when  the  site  was  originally  located  on 
August  14,  1981.  It  took  several  attempts  before  the  pocket  of  organic 
sediment  could  be  located  at  RM  170-2.  The  pocket  was  located  in 
approximately  4  feet  of  water,  8  feet  from  the  shore.  Site  RM  172-1  was 
easily  located  at  the  mouth  of  the  Des  Peres  River.  The  sample  obtained  was 
similar  In  nature  to  that  previously  collected.  Site  RM  172-2  was  beneath  the 
barge  loading  area  where  a  fine  sand  sample,  similar  In  nature  to  that 
previously  collected,  was  again  found. 

The  sediment  was  collected  with  a  ponar  dredge,  mixed  In  a  plastic  tub, 
and  stored  In  glass  jars.  Water  samples  for  use  In  the  kinetic  study  were 
collected  from  the  surface  In  5-gallon  Nalgene  containers.  The  redox 

potential,  temperature,  pH,  and  dissolved  oxygen  were  measured  in  the  field. 

An  additional  sampling  trip  was  made  on  November  20,  1981.  This  trip  Is 
designated  by  *.  Water  and  sediment  samples  were  taken  at  RM  170,  several 
hundred  meters  upstream  from  a  series  of  wing  dams  and  a  barge  docking  area  on 
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the  Illinois  shore.  Sediment  sample  (RM  170*-2)  was  taken  from  the  dredge  cut 
area,  100  feet  upstream  from  the  dredge  in  the  open  channel.  Water  sample  (RM 
170*-3)  was  collected  200  feet  downstream  from  the  dredge.  These  samples  were 
taken  at  approximately  100  feet  west  of  the  dredge  and  well  out  of  the  plume 
area. 


B.  Experimental  Procedure  of  the  Kinetics  Study  The  method  of  the 
kinetic  study  was  basically  the  same  as  that  of  the  elutriate  test.  A  seven 
liter  graduated  jar  was  equipped  with  a  mechanical  stirrer,  air  stone,  and  pH 
and  dissolved  oxygen  probes.  One  liter  of  homogeneous  sediment  was  measured 
into  four  liters  of  site  river  water  by  displacement.  When  the  addition  was 
complete,  a  500  ml  sample  of  the  well-mixed  slurry  was  immediately  extracted 
from  the  surface  (the  time  zero  sample}.  Additional  samples  were  taken  at  the 
following  intervals:  0.1,  0.2,  0.3,  0.5,  1,  2,  6,  and  12  hours.  The  samples 
were  immediately  centrifuged,  filtered  through  a  glass  fiber  filter,  and 
finally  filtered  through  a  0.45  micron  membrane  filter.  The  sediment  was 
retrieved  from  the  centrifuge  tubes,  diluted  to  500  ml  with  site  water,  and 
returned  to  the  original  jar.  An  elutriate  blank  was  prepared.  The  blank  was 
a  sample  of  site  water  run  through  the  filtration  process  and  preserved  with 
nitric  or  sulfuric  acid  in  the  same  fashion  that  the  samples  were  preserved. 
The  samples  were  analyzed  for  ammonia,  COD,  pH,  redox  potential,  chromium,  and 
nickel . 

C.  Development  of  the  Kinetic  Model  The  kinetic  tests  were  performed  in 
a  jar  reactor.  The  system  was  assumed  to  be  a  completely  mixed  batch  reactor 
and  to  be  homogeneous  throughout.  In  this  reactor,  an  initial  limited  number 
of  factors  may  be  Involved  in  the  reactions.  The  relevant  reactions  which 
remove  or  transform  pollutants  Include  adsorption,  desorption,  hydrolysis, 
oxidation,  biodegradation,  photolysis,  and  volatilization.  The  equation 
governing  dissolved  pollutant  concentration  in  the  reactor  may  be  written  as 

dC/dt  -  -  KXCM  +  K2rM  -  (Kp  +  Hh  +  KQ  +  Kb  ♦  Ky)C  (IV-1) 


where 


C  *  dissolved  pollutant  concentration  at  time  t,  ug/l» 
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r  =  adsorbed  pollutant  concentration,  mass  of  pollutant  adsorbed  per 
mass  of  solid,  ug/kg, 
t  *  time,  hours, 

M  =  suspended  solids  concentration,  kg/1, 

adsorption  rate  constant,  (yg/1)’1  (hour)"1, 

K2=  desorption  rate  constant,  hour"1, 

K„=  photolysis  rate  constant,  hour"1, 

hydrolysis  rate  constant,  hour"1, 

K0*  oxidation  rate  constant,  hour"1, 

Ky=  volatilization  rate  constant,  hour"1, 

Kb*  biodegradation  rate  constant,  hour"1. 

It  was  assumed  that  the  rates  of  photolysis,  hydrolysis,  oxidation, 
volatilization,  and  biodegradation  are  either  first  order  or  pseudo-first 
order. 

The  equation  for  pollutant  concentration  adsorbed  by  the  sediment  may  be 
expressed  as 

dr/dt  -  KXC  -  K?r  (IV-2) 

where  r  =  adsorbed  pollutant  concentration,  yg/kg  . 

It  was  assumed  that  no  reactions  occur  In  the  particulate  phase.  Since 
the  sorptive  equilibrium  Is  generally  rapid  compared  to  the  other  reaction 
mechanisms,  It  may  be  assumed  that  the  pollutant  concentration  (at  small 
times)  Is  only  governed  by  adsorption  and  desorption  rates.  If  these 
assumptions  hold  true,  (IV-1)  and  (IV-2)  reduce  to  the  following  adsorption- 
desorption  kinetics: 

dC/dt  =  -  KjCM  +  K2rM  (IV-3) 

dr/dt  =  KjC  -  K2r  (IV-4) 

The  temporal  conditions  for  the  above  equations  are  given  by 

C  *  C*  at  t  3  t* 


(IV— 5) 


Y 


C  =  Ce  at  t  =  00 


(IV— 6 ) 


where 


t  =  initial  mixing  time,  hours, 

C*  =>  Initial  concentration  of  dissolved  pollutant  at  t  ,  ug/1 , 

Ce  =  concentration  of  dissolved  pollutant  at  equilibrium,  pg/l. 

The  analytical  solution  of  (IV-3)  and  ( I V -4 )  under  the  foregoing  conditions, 
is  given  by 


C  =  Ce  +  (C*-Ce)  exp  -  [ (KjH  +  K2)(t-t*)] 


(IV— 7 ) 


P.  Evaluation  of  the  Model  Parameters.  The  parameters  necessary  to 
evaluate  the  model  are:  the  initial  concentration  of  dissolved  pollutant 
(C  ),  the  initial  mixing  time  (t  ),  dissolved  pollutant  concentration  at 
equilibrium  (Ce),  the  suspended- sol  ids  concentration  (M),  the  adsorption- rate 
constant  (Kj),  and  the  desorption- rate  constant  (K2). 

An  attempt  was  made  to  correlate  initial  filtered  ammonia  concentration 
with  sediment  ammonia  and/or  redox  potential.  The  initial  concentrations  of 
the  filtered  ammonia  and  filtered  COD  appear  to  be  correlated  with  the  initial 
redox  potential  (redox  potential  at  t  =  0)  of  the  sample.  Figure  IV-1 
presents  a  plot  of  the  initial  ammonia  concentration  versus  the  corresponding 
initial  redox  potential.  Figure  IV-2  shows  a  similar  plot  of  the  initial 
filtered  COD  versus  the  initial  redox  potential.  Only  three  experiments  were 
performed,  so  the  relationships  in  figures  IV-1  through  IV-4  are  preliminary 
and  somewhat  speculative. 

The  initial  concentrations  of  ammonia  and  COD  also  appear  to  be  a 
function  of  their  concentrations  in  the  bulk  sediments,  which  is  not  always 
the  case  for  other  pollutants.  Plots  of  the  initial  ammonia  concentrations 
versus  sediment  ammonia  concentrations  are  presented  in  figure  IV-3.  Similar 
plots  for  COD  are  shown  in  figure  IV-4.  Although  C*  is  seen  to  be  a  complex 
function  of  the  redox  potential  of  the  solution,  the  expression  of  C*  is 
assumed  to  be  linear  and  of  the  form 
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Initial  COD  (t  =0.1  hrs)  versus  Initial  Redox  Potential  (t  =  0) 


Ammonia  Concentration 


Figure  IV-4  Initial  COD  (t  =0.1  hrs)  versus  Sediment  COD 


Ill 


C*  =  aZ  +  b  (IV-8) 

where  C*  Is  the  dependent  variable  (i.e..  Initial  concentration  of  pollutant), 
Z  is  an  independent  variable,  and  a  and  b  are  constants.  An  attempt  was  made 
to  correlate  initial  ammonia  concentration  with  sediment  ammonia  and/or  redox 
potential.  Several  chemical  parameters  were  tried  for  z  in  (IV-8)  including 
redox  potential.  The  initial  concentration  of  ammonia  (C*)  possessed  a  high 
degree  of  linear  correlation  (r  >  0.99)  with  (NH3)S,  (R')"*,  (R’)"2,  (NH3)S 
(R')"*,  and  (NH3)S  (R’)"2.  A  similar  attempt  was  made  for  COD,  where  R'  =  101 
+  R,  R  =  redox  potential  (mV),  and  (NH3)S  =  sediment  ammonia  concentration.  A 
null  value  of  101  was  added  to  the  original  data  of  redox  potential  to  avoid 
negative  values  of  redox  potential.  (R ' )  had  the  greatest  linear  correlation 
(r  =  0.99)  with  the  initial  COD.  The  final  expressions  for  C*  selected  are: 

C*  =  9045.  (R')"2  +  755.  for  ammonia,  and 

C*  =  0.263  (R’)2  +  45803.  for  COD. 

Note  that  these  relations  possess  log  power  correlations,  and  the  use  of  these 
equations  may  be  limited  in  the  redox  potential  range  from  -100  to  150  mV,  and 
to  the  sediments  and  experimental  conditions  used  in  this  study. 

Approximations  of  the  values  of  Kj,  K2,  and  Ce  were  initially  performed 
using  the  nonlinear  regression  method  (SAS,  1979).  The  calibration  of  the 
model  was  then  accomplished  by  running  a  simulation  and  comparing  the  result 
to  the  corresponding  experimental  data.  Determined  values  of  Kj  and  K2  are 
20.0  (u  g/1)"1  (hour)”*  and  2.0  (hour)-*  for  ammonia,  and  20.0  (ug/1)"* 
(hour)"*  and  5.0  (hour)"*  for  COD,  respectively.  A  summary  of  the  values  for 
the  kinetic  model  parameters  is  shown  in  table  IV-1. 

E.  Results  and  Discussion.  Tables  IV-2,  3,  and  4  summarize  the  results 
of  the  chemical  analyses  of  each  experiment.  Since  the  sample  collected  on 
November  20,  1981,  showed  little  sorption  of  the  chemical  pollutants  during 
elutrlation  (table  IV-5),  the  data  were  not  used  for  the  model  simulation. 

Figures  IV-5  through  IV-7  are  plots  of  ammonia  concentrations  versus  time 
which  show  both  the  experimental  and  computed  results  for  the  three  samples. 


Table  IV-2.  Results  of  the 


NH3 (mg/1 )  COD (mg/1 


sediment  * 

“ira - 

raw  water 

.24 

18 

elutriate  blank 

.24 

10 

Time  (hours 

0 

.82 

48 

0.1 

.75 

43 

0.2 

.82 

28 

0.3 

.85 

29 

0.5 

.88 

42 

1 

.94 

35 

2 

.92 

37 

4 

.88 

34 

6 

31 

12 

.92 

33 

*mg/kg  dry  weight 
suspended  solids  3000  mg/1 


Study  (Site  170-2) 


Redox 

Potential 


Table  IV-3.  Results  of  the  Kinetic  Study  (Site  172-1) 


NH3(mg/l) 


sediment 
raw  water 
elutriate  blank 


WIi 


time  (hours) 
0 


mg/kg  dry  we 
suspended  so 


1 


COO (mg/1) 


Redox 

Potential 

(mV) 


1  18 

+180 

8.02 

!  12 
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Table  IV-4.  Results  of  the  Kinetic  Study  (Site  172-2) 


NH3(rag/l) 

COD (mg/1) 

Redox 

Potential 

(mV) 

pH 

sediment 

16 

+175 

raw  water 

.21 

29 

+210 

8.25 

elutriate  blank 

.10 

29 

time  (hours) 


0 

0.70 

26 

+150 

7.98 

0.1 

0.76 

29 

+190 

0.2 

0.70 

24 

8.21 

0.3 

0.73 

27 

+205 

0.5 

0.68 

26 

+210 

8.37 

1 

0.65 

22 

+210 

8.36 

2 

0.64 

22 

+210 

4 

0.66 

21 

+210 

8.64 

6 

-- 

+210 

8.48 

12 

0.56 

23 

+210 

8.47 

*mg/kg  dry  weight 
suspended  solids  400  mg/1 


IV-5.  Results  of  the  Kinetic  Study  (Site  170  ) 


SITE  170-2 


Figure  IY-5  Dissolved' Ammonia  Concentration  versus  Time  (Site  170-2") 


Figure  IV-6  Dissolved  Ammonia  Concentration  versus  Time  (Site  172-1) 
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170-2,  172-1,  and  172-2,  respectively.  Plots  of  COD  are  shown  in  figures  IV-8 
through  IV-10.  As  can  be  seen,  the  model  closely  approximates  the 
experimental  data  obtained  in  the  kinetic  experiments.  An  immediate  increase 
in  ammonia  and  COD  occurred  within  the  first  six  minutes  of  mixing  time.  The 
concentrations  of  the  pollutants  at  the  time  of  six  minutes  was  defined  as  the 
initial  concentration  of  pollutants.  The  concentrations  of  pollutants  were 
not  constant,  but  generally  decreased  as  the  experiment  progressed  (an 
exception  was  the  ammonia  of  sample  170-2).  The  reduction  of  the 
concentrations  of  the  pollutants  was  presumably  due  to  re-adsorption  which  may 
be  explained  by  the  increasing  redox  potential  with  time.  The  change  in  redox 
potential  can  be  expressed  simply  as  a  function  of  time. 

R  =  R0  +  Ru  (1*°  "  exP  ("M)) 


where 

R  =  redox  potential  at  time  t,  mV, 

R0  =  redox  potential  at  t  =  0,  mV, 

Ru  =  redox  potential  at  equilibrium,  mV, 

kr  =  rate  constant,  hour-1. 

For  example,  in  the  study  using  sample  RM  172-1,  R0  was  found  to  be  -100  mV. 
Ru  and  Kr  were  determined  using  a  nonlinear  regression  method  (SAS,  1979)  and 
found  to  be  326.6  mV  and  6.18  hour-*,  respectively.  The  plots  of  the 
variation  of  redox  potential  with  time  obtained  from  the  experiment  are 
presented,  along  with  the  computed  values,  in  figure  IV-11.  As  can  be  seen 
from  the  results,  a  rapid  increase  in  redox  potential  occurred  within  30 
minutes,  and  thus  the  results  are  in  accord  with  the  adsorption  and  desorption 
trends  of  the  pollutants.  It  should  be  noted  that  according  to  the  kinetic 
models  about  86  percent  of  the  concentration  change  was  complete  within  30 
minutes  and  99.9  percent  within  90  minutes,  while  about  55  percent  of  the 
ammonia  concentration  change  is  complete  within  1  minutes  and  94  percent 
within  90  minutes. 

F.  Suiwary.  The  kinetic  studies  were  very  helpful  in  understanding  the 
standard  elutriate  test.  8asically  the  following  stages  occur: 


SITE  170-2 


Figure  IV-8  Dissolved  COD  versus  Time  (Site  170-2) 


SITE  172-1 


Figure  IV-9  Dissolved  COD  versus  Time  (Site  172-1) 


SITE  172-2 


Figure  IV-10  Dissolved  COD  versus  Time  (Site  172-2) 


Redox  Potential  with  Time  (Site  172-1) 
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INITIAL  CONDITION 

1) .  Reduced  chemicals  such  as  dissolved  manganese  (Mn+2)  and  ammonium 

ion  (NH4+)  and  organics  (measured  partly  as  COD)  are  rich  in  the 
intersticial  waters  of  the  sediment  sample. 

INITIAL  MIXING 

2) .  Within  the  first  few  minutes  of  mixing  (six  minutes  according  to 

the  mathematical  model)  the  intersticial  waters  are  diluted  by  the 
river  water  in  the  elutriate  test.  Further  desorption  from  the 
mixed  sediments  occurs  during  this  stage.  The  redox  potential  of 
the  elutriate  mixture  increases  from  negative  to  positive. 

ADSORPTION/DESORPTION 

3) .  If  the  kinetics  of  desorption  and  mixing  are  fast  enough,  the 

eluted  chemicals  may  "overshoot"  their  final  equilibrium  value  and 
subsequently  readsorb  during  the  period  from  six  minutes  to  30-90 
minutes  in  the  elutriate  test.  After  90  minutes,  sorptive 
equilibrium  is  essentially  complete. 


EQUILIBRIUM 

4).  Sorptive  equilibrium  is  established  after  two  hours.  Any  further 
changes  would  likely  be  caused  by  other  mechanisms.  In  the 
Mississippi  River,  gas  transfer  or  biodegradation  could  ultimately 
reduce  the  concentration  in  the  water. 

Mathematical  modelling  of  the  desorption  kinetics  aided  in  the 
interpretation  of  the  results.  However,  it  is  still  not  possible  to  predict  ± 
priori  the  initial  elutriate  concentration,  C*,  or  the  final  elutriate 
concentration,  Ce.  While  the  elutriate  test  does  not  exactly  reflect  the  many 
environmental  variables  of  the  dredging  operation  and  disposal  (e.g.,  mixing 
intensity,  sediment:  water  ratio,  temperature,  redox  potential),  it  is  a 
useful  laboratory  procedure  to  screen  for  potential  problem  areas.  We  know 
that  organic  sediments  will  have  the  following  characteristics:  a  high 
moisture  content  (>25%),  a  low  percent  total  solids,  a  high  volatile  solids 
content,  a  fine  median  grain  size,  and  a  dark  color.  These  bed  sediments 
generally  occur  in  areas  of  slowly  moving  water,  and  the  redox  potential  is 


perhaps  the  best  Indicator  of  the  NH3,  Mn,  and  COD  that  will  adsorb  in  the 
standard  elutriate  test. 
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V.  CHEMICAL  PLUME  MODEL 

Channel  maintenance  dredging  operations  in  conjunction  with  the 
large-scale  disposal  of  dredged  material  can  promote  a  release  of  associated 
contaminants.  A  goal  of  this  chapter  is  to  develop  a  chemical  plume  model  and 
to  examine  the  extent  and  intensity  of  the  plume  of  dissolved  and  particle- 
associated  chemical  contaminants. 


A,  Description  of  Sawpllng  Site  and  Sampling  Procedures.  A 

sampling  trip  was  made  on  20  November  1981.  The  dredge  disposal  site  was 
located  at  RM  170,  several  hundred  meters  upstream  from  a  series  of  wing  dams 
and  a  barge  docking  area  on  the  Illinois  shore  of  the  Mississippi  River.  The 


sampling  effort  involved  gathering  sediment-laden  water  samples  for  analysis 
of  suspended  solids  and  the  chemical  pollutants  of  interest.  The  water 
samples  were  also  used  for  the  fall  velocity  analysis,  and  bed-sediment 
samples  for  grain-size  distribution  analysis.  These  data  all  serve  as  input 
parameters  for  the  proposed  computer  model. 


The  site  layout  and  sampling  procedures  were  described  in  chapter 
II.  Analytical  procedures  were  described  in  chapter  III. 


B,  Model  Development.  The  concentration  of  pollutant  in  the  water 
column  may  be  expressed  by  the  following  differential  equation. 


3Ct  3C,  .  dCT 

3t“  3  _U1  3x7  +  3x7  (Et  ^  *  R 


(V-l) 
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where  Cy  refers  to  the  total  concentration  (mass  per  unit  volume)  of  pollutant 
present  in  both  solid  and  liquid  phases;  u^  refers  to  the  fluid  velocity  in  a 
rectangular  coordinate  system,  refers  to  the  eddy  diffusion  coefficient 
in  the  i'th  direction;  and  R  refers  to  the  sum  of  the  relevent  reaction  rates 
of  the  pollutant.  The  reactions  may  include  photolysis,  hydrolysis, 
oxidation,  volatilization,  and  biodegradation.  The  rates  of  the  above  reac¬ 
tions  were  assumed  to  be  either  first  order  or  pseudo-first  order.  Consider¬ 
ing  the  above  reaction  rates,  (V-l)  can  be  written  as: 
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where  K  refers  to  the  sum  of  the  reaction  rates  of  the  above  reactions,  and  x, 
y,  and  z  refer  to  the  longitudinal,  vertical,  and  lateral  directions,  respec¬ 
tively.  The  initial  assumptions  are: 

1) .  The  chemical  pollutant  concentration,  Cy,  can  be  described  by  a 

vertically  averaged  pollutant  concentration, 

1  D 

CT  =  0  fQ  Cidy 

where  D  refers  to  the  depth  of  the  water  column  and  y  is  vertical 
distance  in  the  Cartesian  coordinate  system,  x,  y,  z; 

2) .  The  eddy  diffusivities  in  the  x  and  z  directions,  Ex  and  Ez,  are 

independent  of  depth; 

3) .  The  fluid  velocity  in  the  x  and  z  directions,  ux  and  uz,  are  depth 

independent;  and 

4) .  The  pollutant  transport  rates  due  to  vertical  diffusion  and  verti¬ 

cal  fluid  velocity  currents  are  much  slower  than  the  rates  of 
sorption. 


The  resulting  equation  is: 
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The  total  concentration  of  chemical  pollutant  in  the  bulk  water,  Cy, 
Is  the  sum  of  the  concentrations  in  the  particulate  and  dissolved  phases. 

CT=  C  +  Cp  (V-5) 

where  C  refers  to  the  dissolved  pollutant  concentration,  and  Cp  refers  to  the 
adsorbed  pollutant  concentration.  Cp  is  related  to  the  concentration  of 
suspended  solids,  M,  and  is  defined  as 

Cp  =  rM  (V-6) 

where  r  refers  to  the  amount  of  pollutant  adsorbed  per  unit  mass  of  dry  sedi¬ 
ment.  At  a  low  pollutant  level  in  the  dissolved  phase,  linearity  of  the 
Langmuir  isotherm  may  be  assumed,  and  C  may  be  given  by 


C  * 


(V-7) 


where  Kp  refers  to  the  linear  adsorption  partitioning  coefficient.  Adsorption 
and  desorption  rates  are  orders  of  magnitude  greater  than  the  rates  of  the 
decay  and  transport  processes  considered,  and  instantaneous  sorptive  equilib¬ 
rium  is  achieved  between  the  two  phases  very  rapidly  in  comparison  with  sedi¬ 
mentation,  volatilization,  photolysis,  hydrolysis,  oxidation,  and  biological 
degradation.  Therefore,  C  and  Cp  may  be  replaced  by  their  equivalents  in 
terms  of  CT.  Incorporating  (V-6)  and  (V-7)  into  (V-5)  and  expressing  C  and  Cp 
in  terms  of  Cy  under  the  conditions  of  sorptive  equilibrium  in  which  the 
concentration  of  suspended  solids,  M,  is  known. 


C  •  r  +  XpM 

K  MCt 
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(V— 8 ) 
(V-9) 


It  is  surmised  that  the  reactions  of  photolysis,  volatilization,  oxidation, 
and  biodegradation  occur  in  the  dissolved  phase,  and  that  the  removal  of 
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adsorbed  chemical  pollutants  is  controlled  by  sedimentation  of  suspended 
solids.  Thus,  KCj  can  be  described  by 


KCT  =  KsCp  +  kC 


(V-10) 


where  K$  is  the  chemical  pollutant  removal  rate  due  to  sedimentation  of  sus¬ 
pended  solids,  and  k  refers  to  the  sum  of  all  decay- react  ion  rates  of 
dissolved  chemicals.  Substitution  of  (V-8)  and  (V-9)  into  (V-10)  gives 


KCt  =  ( 
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where  Mj  refers  to  the  concentration  of  the  j-size  fractions  of  suspended 
solids;  K$j  and  Kp j  are  the  constants  corresponding  to  Mj.  The  concentrations 
of  the  suspended  sediment  may  fall  into  three  size  fractions;  Mj,  M2,  and  M3 
which  refer  to  sand,  silt,  and  clay,  respectively.  Substituting  (V-ll)  for 
KCj  into  (V-4)  gives  the  differential  equation: 


-u 


3Ct  3Ct  a 
"uz  Tz~  +  i7^E 


x  3  x 


3C-| 
z  Tz 


■) 


,KslKolMl  +  Ks2Kd2M2  +  KS3K03M3  +  k 
" 1 - 1  '+  k  +  k  7m;  '+  K 


pl'l  >2' '2  p3  3 


(V-12) 


In  a  river,  the  following  additional  assumptions  can  be  made  to  further 
simplify  the  above  equation  (Sayre,  1973). 
ux  »  constant, 
uz  =  0., 

D  =  constant, 

Ez  s  constant,  3C 
Ex  «  Ez,  so  |7(Ex  yt)  =  0 


Incorporating  these  assumptions  and  assuming  steady  state,  OC^/Dt  =  0,  Equa¬ 
tion  (V-12)  can  be  written  in  the  form 


3Cj 
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Ez  ^CT 


(V-13) 


The  solution  to  (V-13)  at  a  given  constant  Mj  can  be  written  (Sayre,  1979)  as: 


1 
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Cj(z»x)  = 

,  kc1knlMl  +  K-K  +  KJC  -M,  +  k 

cT<->-  -< si  %T  * plv-  — 1 1 

where  C'(z,x)  Is  the  solution  of  the  diffusion  equation: 

3Cy  Ez  32Ct 
Tx~  ~  u“  „  2 


(V-14) 


(V-15) 


For  the  case  of  a  continuous  point  source  of  flow,  Q0,  and  chemical  concen¬ 
tration,  C0,  the  solution  of  (V-15)  is  (Sayre,  1973): 
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This  equation  has  the  form  of  a  normal  probability  function  with  variance, 
2  _  x 
°z  -  2Ez  V 

Substituting  this  relationship  into  (V-16)  gives: 


Q  C  .  2 

CT  (z»x)  =  TTd2  - exp 

'  V  /3T  oz  2oz 


(V-17) 


Equation  (V-17)  describes  the  plume  resulting  from  a  continuous  point 
source.  The  mode  of  dredge  disposal  which  will  be  discussed  later,  can  be 
better  described  by  a  line  souce  of  width,  b,  perpendicular  to  the 
streamline.  Equation  (V-17)  can  be  modified  to  ('escribe  a  line  source  by  the 
method  used  by  Sayre  (1973,  1979).  The  resulting  expression  is: 


4(z,x)  -  /  C'(I*,0)  — 1 


— - —  exp  [-  )  3  dz ' 
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(V- 18) 


where  B  is  the  width  of  the  river,  and  z’  Is  a  dummy  variable  describing  any 
point  within  the  source  width.  The  initial  conditions  for  this  line  source 


are: 


132 


Q  C 

C '  (z  *  ,0 )  =-2-2  0  <  z'  <  b 

wb 
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(V— 19 ) 


where  Qb  =  uxDb  is  the  portion  of  the  river  flow  passing  through  the  source 
width,  b.  Incorporating  these  into  (V-18)  and  substituting  the  standard 
normalized  variable  s  =  z  -  z' /o  gives 


C'(z,x) 
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which  is  in  the  form  of  the  cumulative  normal  distribution  function.  The 
solution  of  (V-20)  is 


C|(*.x)  =  (F(|_,  .  F(ilfe,) 

yb  z  °z 

where  the  value  of  F(*)  can  be  obtained  from  a  cumulative  normal 
table  or  computed  by  the  following  equations  (Kennedy  and  Gentle, 


(V-21 ) 

distribution 

1980). 
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The  suspended-solids  plume  described  by  (V-21)  includes  no  effects 
from  the  side  banks  of  the  river.  It  is  assumed  that  the  channel  banks  act  as 
reflecting  barriers.  Including  the  effects  of  reflection  from  the  near  side 
bank,  the  equation  becomes; 

CT<Z>*>  *  -  F(i-j-H))  (V-22) 

This  equation  Is  not  applicable  If  the  suspended  solids  plume  disperses  In  the 
lateral  direction  enough  to  reflect  from  the  far  shoreline.  Substituting 
(V-22)  into  (V-14)  gives  the  final  solution. 
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Equation  (V- 23)  can  be  solved  simultaneously  with  the  suspended- sol ids-plume 
equation  proposed  by  Schnoor  et  al.  (1980). 

C.  Results  and  Discussions.  The  plumes  of  dissolved  and  particu¬ 
late-associated  chemical  pollutant  are  described  by  the  simultaneous  computa¬ 
tion  of  two  fields  at  a  grid  point  of  the  system  under  consideration:  a  field 
where  the  suspended-solids  concentration  is  described  by  a  processes  of  con¬ 
vection,  dispersion,  and  settling,  and  a  field  where  the  concentration  of  a 
chemical  pollutant  is  described  by  processes  of  adsorption-desorption  equi¬ 
librium  and  pollutant  reactions.  The  process  of  the  first  field  is  translated 
into  a  change  in  the  distribution  of  suspended  solids,  which  is  related  to  the 
chemical  concentrations  in  the  second  field. 

1.  Evaluation  of  the  Model  Parameters.  The  behavior  of  the  river  is 
the  result  of  interaction  between  the  flow,  the  quantity  and 
characters  of  the  sediment  discharge,  the  local  topography,  and 
local  soil  and  geological  conditions  (Shen,1979).  Because  of  the 
extremely  complicated  behavior  of  the  river,  techniques  for  model¬ 
ing  invariably  involve  some  averaging  processes.  Some  of  the 
parameter  values  are  directly  determined  in  the  field 
measurements,  while  others  are  estimated  by  calibrations  of  the 
model.  The  model  simulate  a  chemical  plume  developed  during  the 
dredging  operation  at  site  RM  170  on  20  Nov.  1981. 

a).  Mean  Depth  of  the  Plume,  D.  The  plumes  resulting  from  dred¬ 
ging  operations  are  gererally  developed  in  a  near-shore  area 
where  the  depth  varies  considerably.  Since  the  model  assumes 
a  constant  depth,  depths  were  measured  at  several  locations 
in  the  area  of  the  plume.  An  average  depth  of  4  meters  W3S 
used  in  the  model. 


b) .  Mean  Stream  Velocity,  ux.  The  stream  velocity  varies  not 

only  with  location  but  also  with  depth.  Therefore,  the 
velocity  profiles  in  the  immediate  vicinity  of  the  plume  are 
very  useful  to  determine  the  mean  velocity.  A  representative 
velocity  of  0.08  m/sec  was  used  for  the  model  simulation. 

c) .  Dispersion  Coefficient,  Ez.  The  dispersive  term  includes 

molecular  diffusion,  turbulent  diffusion,  and  mixing.  The 
dispersion  coefficient  lumps  together  all  these  effects.  A 
velocity  of  greater  than  3.0  ft/sec  at  the  disposal  site  at 
the  surface  created  a  great  deal  of  turbulence  and  large 
eddies.  These  factors  are  contained  in  the  overall  disper¬ 
sion  coefficient.  Allowing  for  all  these  effects,  a  value 
Ez  =  25  m2/sec  was  chosen. 

d) .  Width  of  the  Plume  Source,  b.  The  width  of  the  plume  source 

is  a  function  of  the  discharge,  velocity,  flow  direction,  and 
momentum  of  the  dredge  effluent.  The  mode  of  discharge  of 
dredged  material  is  so  designed  that  the  spoil  slurry,  car¬ 
ried  through  the  discharge  pipe,  hits  a  deflector  plate, 
splashes,  and  falls  through  several  meters  before  it  hits  the 
water  surface.  The  plume  resulting  from  this  mode  of 
disposal  operation  would  be  expected  to  have  a  wider  plume 
source  than  a  plume  resulting  from  side  bank  disposal. 
Considering  the  high-velocity  discharge  in  the  lateral  direc¬ 
tion,  a  plume  width  of  130  meters  was  chosen. 

e) .  Background  Suspended  Solids  and  Chemical  Concentrations.  The 

ambient  suspended-solids  concentration  at  the  site  was  84 
mg/1.  The  ambient  concentrations  of  chemicals  were:  0.03 
mg/1  for  ammonia-nitrogen,  31  mg/1  for  COD,  and  0.16  mg/1  for 
manganese.  Variations  between  ambient  site-water  concentra¬ 
tions  were  +  0.002  mg/1  for  NH^-N,  +4  mg/1  for  COD,  and 
+  0.04  mg/1  for  P0^-  P.  Since  the  phosphorus  concentration 
was  measured  as  filtrable  phosphorus,  a  total  phosphorus 
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concentration  of  0.6  mg/1  was  assumed  based  on  the 
partitioning  coefficients  of  13.0  for  sand,  130.0  for  silt, 
and  1300.0  (ug/kg) (y g/1  )"*  for  clay.  The  concentrations  of 
zinc,  copper,  '■hromium,  and  lead  were  reported  as  less  than 
(<)  0.001  mg/1  and  the  concentration  of  nickel  as  <  0.1 
mg/1.  An  assumed  concentration  of  0.005  mg/1  was  used  for 
all  these  trace  metals. 

The  model  involves  three  size  fractions;  sand, 
silt,  and  clay.  The  proportion  of  each  fraction  was  deter¬ 
mined  from  the  size-analysis  results.  The  composition  of  the 
background  suspended  solids  was  approximately  0  percent  sand, 
50  percent  silt,  and  50  percent  clay.  It  was  assumed  that  no 
background  suspended  solids  are  subject  to  settling  within 
the  plume,  a  relatively  short  period. 

f).  Initial  Suspended  Solids  and  Chemical  Concentrations  above 
Ambient.  The  initial  concentration  above  ambient  can  be 
calculated  by 

C.  =  Q  (C  -  C  )/(Q  +  Q) 
in  yov  o  a;/'yo  y/ 

where  Q0  =  flow  discharge, 

Q  =  receiving  water  flow  rate, 

CQ  =  discharge  concentration  of  suspended  solids  or  chemical 
pollutant, 

Ca  *  ambient  concentration  of  suspended  solids  or  chemical 
pollutant. 

Because  of  the  high  velocity  and  turbulent  river  conditions, 
it  is  difficult  to  correlate  the  elutriate  test  results  with 
the  Initial  chemical  concentration,  C0.  Therefore,  it  was 
decided  to  estimate  a  value  for  each  chemical  pollutant 
through  the  calibrations  of  the  model.  The  initial 
suspended-solids  concentration  above  ambient  was  estimated  to 
be  5  mg/1.  The  composition  of  the  suspended  solids  was 
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approximately  70  percent  sand,  15  percent  silt,  and  15 
percent  clay.  The  proportion  of  each  fraction  was  determined 
from  the  size  analysis  of  the  bed  sediments.  The  initial 
concentrations  estimated  for  each  chemical  contaminant 
were:  0.15  mg/1  for  ammonia-nitrogen,  0.1  mg/1  for  total 

phosphorus,  3.0  mg/1  for  COD,  0.005  mg/1  for  manganese,  and 
0.001  mg/1  for  zinc,  copper,  chromium,  lead,  and  nickel. 

g) .  Settling  Velocity  of  the  Suspended  Solids,  W 

The  terminal  settling  velocity  of  a  suspended 
particle  is  determined  by  its  specific  gravity  and  size.  The 
particle  size  can  be  estimated  by  a  fall  velocity  analysis. 
Knowing  the  particle  size  and  assuming  a  specific  gravity,  a 
modification  of  Stoke' s  law  can  be  used  to  determine  the 
terminal  velocity.  The  settling  velocities  of  the  sand, 
silt,  and  clay  were  respectively  determined  as  0.02,  0.003, 
and  0.000001  m/sec.  The  settling  rate  constant,  Ks,  was 
calculated  as  W/D,  where  W  is  the  settling  velocity  of  the 
suspended  solids  and  0  is  the  mean  depth  of  the  plume. 

h) .  Partitioning  Coefficient,  Kp 

It  was  assumed  that  when  a  dissolved  chemical 
pollutant  is  released  in  water  which  contains  particulates, 
the  pollutant  rapidly  reaches  equilibrium.  Partitioning 
coefficients  were  estimated  based  on  the  results  presented  in 
the  GREAT  II  report  (Schnoor  et  al.,  1980).  The  estimated 
values  for  the  chemicals  of  interest  are:  13.0  for  sand, 
130.0  for  silt,  and  1300.0  (u  g/kg )  (u  g/1 )  ” 1  for  clay 
particles. 

i) .  total  Decay-Reaction  Rate  Constant,  k 

The  reactions  by  which  chemicals  are  removed  or 
transformed  are  hydrolysis,  biodegradation,  oxidation, 
photolysis,  and  volatilization.  These  reactions  were  summed 
to  yield  the  total  chemical  decay  rates  in  the  liquid 
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phase.  All  the  reactions  were  assumed  to  be  first  order  in 
chemical  concentrations.  Since  the  previous  kinetics  study 
showed  that  the  rates  of  these  reactions  are  insignificant  as 
compared  to  the  rates  of  sorption,  diffusion,  and  convective 
transport,  a  null  value  of  10“^/s  was  used  for  k. 

2.  Calibrations  of  the  Chemical  Plume  Model.  As  discussed  in  Chapter 
II,  grain-size  distribution  of  the  bed  sediment  sample  showed  that 
there  was  very  little  silt  and  clay  present.  Since  the  dredge 
spoil  was  primarily  sand,  which  has  the  least  sorption  capacity, 
it  settled  out  quickly  at  the  disposal  site.  Only  a  small  portion 
of  chemicals  were  removed  by  the  settling  of  sand.  The  quantities 
of  silt  and  clay,  which  have  larger  adsorption  capacities,  were 
small  and  were  not  subject  to  rapid  settling.  As  a  result  the 
decay  of  the  chemical  pollutants  associated  with  particle  sedimen¬ 
tation  was  not  very  substantial,  and  a  large  portion  remained  in 
the  dissolved  and  suspended  phases.  The  dissolved  and  small- 
particulate  associated  chemicals  behaved  like  a  fluid,  and  thus 
dispersed  and  were  transported  a  considerable  distance  prior  to 
settling. 

Calibration  of  the  proposed  chemical  plume  model  was 
performed  using  the  data  obtained  from  site  RM  170*.  Based  on  the 
model,  three  parameters  were  investigated.  An  increase  of  the 
dispersion  coefficient,  Ez,  was  investigated  first.  Increasing  Ez 
caused  a  wider  plume,  but  it  would  decrease  chemical  concentration 

O 

rapidly  downstream.  Using  Ez  =  25  nr/sec,  the  partitioning 
coefficients,  Kpj,  Kp2»  and  Kp3,  were  investigated  next.  The 
partitioning  coefficients  were  not  a  sensitive  parameter  in  the 
model  because  sedimentation  of  silt  and  clay  particles  was  not 
substantial.  Finally,  increasing  the  width  of  the  line  source,  b, 
was  investigated.  By  increasing  b  up  to  130  m,  reasonable  results 
were  obtained.  Since  a  factor  of  the  initial  mixing  force  to 
develop  the  initial  plume  width  and  depth  is  not  incorporated  in 
the  model,  Increasing  b  may  be  a  fair  adjustment  to  the  model. 

The  field  data  for  suspended  solids,  ammonia, 
phosphorus,  COD,  manganese,  chromium,  lead,  copper,  zinc,  and 


nickel  are  presented  in  table  V-l.  The  simulated  chemical  plumes 
for  ammonia  nitrogen  are  shown  in  table  V-2  and  for  chemical 
parameters  are  shown  in  appendix  C.  Figures  V-l,  V-2,  and  V-3 
show,  respectively,  the  simulated  plumes  for  ammonia,  phosphorus, 
and  COD  superimposed  on  the  field  data.  Manganese,  zinc,  copper, 
chromium,  lead,  and  nickel  did  not  develop  plumes  and  were  not 

plotted.  From  figures  V-l  through  V-3  and  table  V-l,  it  is  seen 
that  lower  chemical  concentrations  were  detected  from  the  first 
transect  than  those  from  other  transects.  Reasons  for  the  low 

concentrations  immediately  downstream  from  the  discharge  point 
were  explained  for  the  suspended-solids  plume  in  chapter  II. 

Since  chemicals  behave  similarly  to  suspended  solids,  the  same 
discussion  could  be  applied  to  the  chemical  plume. 

As  can  be  seen  from  the  computer  simulation  results,  the 
model  concentrations  correspond  reasonably  well  with  the  chemical 
analyses.  However,  only  in  the  case  of  ammonia  nitrogen,  was 

there  a  clearly  discernible  input  to  the  water  column  from  dredge 
disposal.  This  disposal  operation  did  not  result  in  the  formation 
of  significant  chemical  plumes,  except  in  the  case  of  ammonia. 

With  the  given  river  conditions,  the  extent  and 
intensity  of  the  chemical  plume  was  preceeded  by  the  computation 
of  the  suspended-solids  plume.  In  this  way,  the  transport  and 
fate  of  dissolved  and  particulate-associated  chemicals  in  the 
river  were  satisfactorily  described,  even  though  the  model 
coefficients  might  need  further  confirmation. 

Dilution  Model.  The  dlultion  model  proposed  by  Schnoor  et  al . 
(1980)  was  used  to  examine  the  dilution  effect  of  the  receiving 
water,  and  to  predict  the  concentration  of  dissolved  chemicals  in 
the  plume.  This  approach  provided  insight  Into  the  dilution 
mechanisms  of  the  river.  The  major  assumption  involved  in  this 
method  is  that:  chemical  pollutants  previously  adsorbed  to 

sediment  particles  desorb  and  dissolve  in  the  discharge  pipe 
during  the  dredging  operation.  Once  these  dissolved  substances 
return  to  the  river,  they  are  subject  to  neither  settling  nor 


Table  V-l.  Field  chemical  plume  data  for  suspended  solids,  ammonia,  phosphorus,  COD,  Manganese,  Zinc,  Copper, 
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Table  V-2.  Simulated  Ammonia  Plume  -  CPLUME 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(9) 

(?) 

(10) 

(11) 

(12) 

(13) 

DISTANCE  WIDTH 

SAND 

SILT 

CLAY 

TOTAL 

TSS 

POLL. 

POLL. 

POLL. 

DISSOLVED 

TOTAL 

POLL. 

ABOVE 

ON 

Oh 

ON 

ABOVE 

X 

2 

SS. 

AHBI. 

SAND 

SILT 

CLAY 

POLL. 

POLL. 

AHBI. 

(M) 

(a) 

(HG/L) 

(HG/L) 

(HG/L) 

(HG/L) 

(HG/L) 

(HG/L) 

(HG/L) 

(HG/L) 

(HG/L) 

(MG/L) 

(HG/L) 

20 

0.0 

3.0? 

41.75 

41.75 

86.59 

4.59 

0.00 

0.01 

0.13 

0.04 

0.18 

0.15 

20 

7.1 

3.09 

41.75 

41.75 

86.5? 

4.5? 

0.03 

0.01 

0.13 

0.04 

0.18 

0.15 

20 

14.1 

3.09 

41.75 

41.73 

86.58 

4.58 

0.00 

0.01 

0.13 

0.04 

0.18 

0.15 

20 

21.2 

3.09 

41.75 

41.75 

86. 58 

4.58 

0.00 

0.01 

0.13 

0.04 

0.18 

0.15 

20 

28.3 

3.08 

41.75 

41.73 

86.58 

4.58 

0.00 

0.01 

0.13 

0.04 

0.18 

0.15 

20 

35. 4 

3.08 

41.75 

41.75 

86.57 

4.57 

0.00 

0.01 

0.13 

0.04 

0.18 

0.15 

20 

38.9 

3.07 

41.74 

41.75 

86.06 

4.56 

0.00 

0.01 

0.13 

0.04 

0.18 

0.15 

20 

42.4 

3.07 

41.74 

41  15 

86.55 

4.56 

0.00 

0.01 

0.13 

0.04 

0.18 

0.15 

20 

44.0 

3.06 

41.74 

41.74 

86.55 

4.55 

0.00 

0.01 

0.13 

0.04 

0.18 

0.15 

20 

49.5 

3.05 

41.74 

41.74 

86.54 

4.54 

0.00 

0.01 

0.13 

0.04 

0.18 

0.15 

20 

54.6 

3.03 

41.73 

41.74 

86.50 

4.50 

0.00 

0.01 

0.13 

0.04 

0.18 

0.15 

20 

63.6 

3.00 

41.73 

41.73 

86.45 

4.45 

0.00 

0.01 

0.12 

0.04 

0.18 

0.15 

20 

70.7 

2.94 

41.71 

41.71 

o6.37 

4.37 

0.00 

0.01 

0.12 

0.04 

0.17 

0.14 

20 

89.4 

2.72 

41.66 

41.66 

86.04 

4.04 

0.00 

0.01 

0.11 

0.04 

0.16 

0.13 

20 

106.1 

2.32 

41.56 

41.56 

85.44 

3.44 

0.00 

0.01 

0.10 

0.03 

0.14 

O.tl 

20 

123.7 

1.76 

41.43 

<1.43 

84.62 

2.62 

0.00 

0.01 

0.08 

0.03 

0.12 

0.09 

20 

141.4 

1.15 

41.28 

41.28 

83.71 

1.71 

0.00 

0.01 

0.06 

0.02 

0.0V 

0.06 

20 

176.8 

0.29 

41.07 

41.07 

82.43 

0.43 

0.00 

0.00 

0.03 

0.01 

0.04 

0.01 

50 

0.0 

2.51 

41.73 

41,73 

85,98 

3.98 

0.00 

0.01 

0.13 

0.04 

0.18 

0.15 

50 

11.2 

2.50 

41.73 

41.73 

85,97 

3.97 

0.00 

0,01 

0.12 

0.04 

0.18 

0.15 

50 

22.4 

2.48 

41.72 

41,73 

85.93 

3,93 

0.00 

0.01 

0.12 

0.04 

0.18 

0.15 

50 

33.5 

2.45 

41.71 

41.72 

85.83 

3.83 

0.00 

0.01 

0.12 

0.04 

0.17 

0.14 

SO 

44.7 

2.40 

41.70 

41.70 

85.80 

3.80 

0.00 

C.01 

0,12 

0.04 

0,17 

0.14 

50 

55.9 

2.32 

41.68 

41.68 

85,68 

6,63 

0.00 

0.01 

0.12 

0.04 

0,17 

0.14 

50 

61.5 

2.28 

41.66 

41.67 

85.61 

3,61 

0.00 

0.01 

0.12 

0.04 

0.16 

0.13 

SO 

67.1 

2.23 

41.65 

41.65 

65.53 

3.53 

0.00 

0.01 

0.11 

0.04 

0.16 

0.13 

50 

72.7 

2.17 

41.63 

41.64 

85.44 

3.44 

0.00 

0.01 

0.11 

0.03 

0.16 

0.13 

SO 

78.3 

2.11 

41.61 

41.62 

85,34 

3.34 

0.00 

0.C1 

0.11 

0.03 

0.15 

0.12 

50 

89.4 

1.96 

41.57 

41.57 

35.11 

3.11 

0.00 

0.01 

0.10 

0.03 

0.14 

0.11 

50 

100.6 

1.79 

41,52 

41.53 

64.84 

2.84 

0.00 

0.01 

0.10 

0.03 

0.14 

0.11 

50 

111.8 

1.61 

41.47 

41.47 

84.55 

2.55 

0,00 

0.01 

0.09 

0.03 

0.12 

0.09 

50 

139.8 

1.10 

41.32 

41.32 

83.75 

1.75 

0.00 

0.01 

0.07 

0.02 

0.09 

0.06 

50 

167.7 

0.64 

41,19 

41.19 

83.01 

1,01 

0,00 

0.00 

0.05 

0.02 

0.07 

0.04 

50 

195.7 

0.31 

41.0? 

41.0? 

62.49 

0.49 

0.00 

0.00 

0.03 

0.01 

0.05 

0.02 

50 

223.6 

0.12 

41.04 

41.04 

82,1? 

0.1? 

0.00 

0.00 

0.03 

0.01 

0.04 

0.01 

50 

279.5 

0.01 

41.00 

41.00 

52.02 

0.02 

0.00 

o,:-> 

0.02 

.01 

0.03 

0.00 
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Table  V- 2.  (Continued) 


too 

0.0 

1.69 

41.47 

41.67 

85.03 

3.03 

0.09 

0.01 

0.12 

0.04 

9.16 

0.13 

200 

15.8 

1.67 

41.66 

41.67 

85.01 

3.01 

0.00 

0.01 

0.12 

0.04 

0.16 

0.13 

100 

31.4 

1.64 

41.65 

41 .65 

84.94 

2.94 

0.00 

0.01 

0.11 

0.04 

0.16 

0,13 

100 

47.4 

1.57 

41.62 

41.63 

84.83 

2.83 

0.00 

0.01 

0.11 

0.03 

0. 16 

0.13 

100 

43.2 

1.49 

11.59 

41.60 

84.67 

2.67 

0.00 

0.01 

0.11 

0.03 

9.15 

0.12 

100 

79.1 

1.38 

41.55 

41.55 

84.48 

2.48 

0.03 

0.01 

0.10 

0,03 

0.14 

0.11 

100 

87.0 

1.32 

41.52 

41.53 

84.37 

2.37 

0.00 

0.01 

0.10 

,.03 

0.14 

0.11 

ICO 

94.9 

1.25 

41.50 

41.50 

84.25 

2.25 

0.00 

0.01 

0.09 

0.03 

0.13 

0.10 

100 

102.8 

1.19 

41.47 

41.47 

84.13 

2.13 

0.00 

0.01 

0.09 

0,03 

0.12 

0.99 

100 

110.7 

1.12 

41.44 

41.45 

84.00 

2.03 

3.03 

0.01 

0.03 

0.03 

0.12 

0.09 

100 

126.5 

0.97 

41.38 

41.39 

83.74 

1.74 

0.00 

0.C1 

0.98 

3,02 

0.11 

0.08 

100 

142.3 

0.82 

41.33 

41.33 

83.47 

1*47 

0.00 

0.01 

0.07 

0.02 

0,10 

0,07 

100 

158.1 

0.68 

41.27 

41.27 

03.22 

1.22 

0.00 

0.01 

0.06 

0.02 

0.08 

O.Ci 

100 

197.4 

0.37 

41.15 

41.15 

82.66 

0.66 

0.00 

0.00 

0.04 

3.01 

0.06 

0.03 

100 

237.2 

0.16 

41.07 

41.07 

82.29 

0.29 

0.00 

0.00 

0.03 

0.01 

0.04 

0.91 

100 

274.7 

0.06 

41.02 

41.02 

62.11 

0.11 

0.00 

0.00 

0.02 

3.01 

0.03 

3.00 

100 

316.2 

0.02 

41.01 

41.01 

82.03 

0.03 

o.cc 

0.00 

0.02 

3, Cl 

0.03 

C.CO 

100 

395.3 

0,00 

41.00 

41.00 

82.00 

0.00 

0.00 

0,03 

0.02 

3.01 

0.03 

0.00 

250 

0.0 

0.51 

41.51 

41. S3 

83.56 

1.56 

0.00 

0,01 

0.10 

0.03 

0.14 

0.11 

250 

25.0 

0.51 

41.51 

41.52 

83.53 

1.53 

0.03 

0.01 

0.09 

0.03 

0.13 

0.10 

250 

50.0 

0,49 

41.49 

41.50 

83.47 

1.47 

0.00 

0.01 

0.09 

0.03 

0.13 

0.10 

250 

75.0 

0.45 

41.45 

41.46 

83.37 

1.37 

0.00 

0,01 

0.09 

0.03 

0.12 

0.09 

250 

100.0 

0.41 

41.41 

41.42 

83.25 

1.25 

0.00 

0.01 

0.08 

0.03 

0.11 

0.08 

250 

125.0 

0.36 

41.36 

41.37 

83.10 

1.10 

0.00 

0.01 

0.07 

0,02 

0.10 

0.07 

250 

137.5 

0.34 

41.34 

41.34 

83.02 

1.02 

0.00 

0.01 

0.07 

0.02 

0.10 

0.07 

250 

150.0 

0.31 

41.31 

41.32 

82.94 

0.94 

0.00 

0.01 

0.07 

0.02 

0.09 

0.06 

250 

162.5 

0.28 

41.28 

41.29 

82.86 

0.86 

0.00 

0.01 

0.06 

0.02 

0.09 

0.06 

250 

175.0 

0.26 

41.26 

41.26 

82.78 

0.78 

0.00 

0.01 

0.06 

0.02 

0.08 

5.05 

250 

2C0.0 

0.21 

41.21 

41.21 

82.63 

0.63 

0,00 

0.01 

0.05 

0.02 

0.07 

0.04 

250 

225.0 

0.16 

41.16 

41.17 

32.49 

0.49 

0.00 

0,00 

0.04 

0.01 

0.06 

0.03 

250 

250.0 

0.12 

41.12 

41.13 

82.37 

0.37 

0.00 

0,00 

0.04 

0.01 

0.06 

0.03 

250 

312.5 

0,05 

41.05 

41.05 

82.16 

0.16 

0.00 

0.00 

0.03 

0.01 

0.04 

0.01 

250 

375.0 

0.02 

41.02 

41,02 

82.06 

0.06 

0.00 

0.00 

0.02 

0.01 

0.03 

0.00 

250 

437.5 

0.01 

41.01 

41.01 

32.02 

0.02 

0.09 

0,00 

0.02 

0.01 

0.03 

0,03 

250 

500.0 

0.C9 

41.00 

41.00 

82,00 

0.00 

o.oc 

0,00 

0.02 

0.01 

0.03 

0.00 

250 

625.0 

0.00 

41.00 

41.00 

82.03 

0.00 

0.00 

0,00 

0.02 

C.01 

0.03 

0.00 

LONGITUDINAL  DISTANCE,  X,  METERS 


145 


DISSOLVED  PHOSPHORUS-P  PLUME  [mg/IJ 

MISSISSIPPI  RIVER  -  RM  170 
AMBIENT  DISS.  PHOSPHORUS-P  =  0.06-O.H  mg/1 


.14  .14  .14  .14  .14  .1« 


.1*  .11  .1*  .11  .1*  -1* 


.It  .11  .IS  .It  .It 


.It  .It  .11  -It  •'* 


-soo-J  o.t4  mg/l 


0.13  mg/l 


8  DISPOSAL  POINT 
X  SHORE 


LATERAL  DISTANCE,  Z,  METERS 


Figure  V-2.  Simulated  Dissolved  Phosphorus  Plume  Superimposed  on  Field  Data 
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decay  in  the  near  field,  but  they  are  affected  only  by  dilution 
with  the  receiving  water. 

The  dilution  ratio  is  defined  as  the  number  of  volumes 
of  river  water  added  to  one  volume  of  water  at  the  head  of  the 
plume  which  yields  the  concentration  in  the  plume.  The  dilution 
ratio  can  be  calculated  by 

DR  *  Cb/C'  (V-24 ) 

where  Cb  is  the  concentration  at  the  head  of  the  plume  given  by 


Cb  =  9oCo/Qb 


(V-25) 


where  is  the  river  flow  rate;,  C0  is  the  dissolved  pollutant 
concentration  discharged  to  the  river;  C*  is  the  concentration  of 
the  dissolved  pollutant  in  the  plume,  given  by  (V-22).  Thus,  the 
dilution  ratio  (DR)  can  be  rewritten  as 


0R  ‘  F<Ct>  -  F(S±) 

°z  az 

The  concentration  of  the  dissolved  substance, 
calculated  by 


(V- 26) 

C',  is  then 


C'  =  Cb/DR  (V-27 ) 

Table  V-3  shows  the  results  of  calculation  of  the 
dilution  factors  and  predicted  chemical  concentrations.  It  can  be 
seen  from  the  tables  that  the  concentration  of  a  dissolved 
substance  calculated  by  the  Dilution  Model  Is  almost  identical  to 
the  concentrations  calculated  by  the  Chemical  Plume  Model.  This 
agreement  indicates  that  dilution  is  the  prime  mechanism  operating 
to  control  the  extent  and  Intensity  of  the  chemical  plume,  and 
that  sedimentation  played  only  a  small  part  at  site  RM  170*. 

D.  Siwary.  Table  V-3  comprises  a  summary  of  the  chemical  plume 
modeling  effort  and  Its  relationship  to  the  standard  elutriate  test.  The 
following  observations  can  be  made  from  table  V-3. 


Table  V-3.  Simulated  Ammonia  Plume  -  DILUTION 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

X 

Z 

F<«) 

F(-t) 

DILUTION 

RATIO 

POLL.  CCNC 

(H) 

(ID 

(MG/L) 

20. 

0,0 

0.9999 

0.0001 

1.0 

0.15 

20. 

7.1 

0.9999 

0.0003 

1.0 

0.15 

20. 

14.1 

1.0000 

0.0005 

1.0 

0.15 

20. 

21.2 

1.0000 

0.0010 

1.0 

0.15 

20. 

28.3 

1.0000 

0.0020 

1.0 

0.15 

20. 

35.4 

1.0000 

0.0017 

1.0 

0.15 

20. 

38.9 

1.0000 

0.0050 

1.0 

0.15 

20. 

42.4 

1.0000 

0.0066 

1.0 

0.15 

20, 

46.0 

1.0000 

0.0087 

1.0 

0.15 

20, 

49.5 

1.0000 

0.0114 

1.0 

0.15 

20. 

56.6 

1.0000 

0.0189 

1.0 

0.15 

20. 

63.6 

1.0000 

0.0303 

1.0 

0.15 

20. 

70.7 

1.0000 

0.0468 

1.0 

0.14 

20. 

88.4 

1.0000 

0.1196 

1.1 

0.13 

20. 

106.1 

1.0000 

0.2492 

1.3 

0.11 

20. 

123.7 

1.0000 

0.4298 

1.8 

0.09 

20. 

141.4 

1.0000 

0.6267 

2.7 

0.06 

20. 

176.8 

1.0000 

0.9071 

10.8 

0.01 

50. 

0.0 

0.9900 

0.0100 

1.0 

0.15 

50. 

11.2 

0.9942 

0.0168 

1.0 

0.15 

50. 

22.4 

0,9968 

0.0271 

1.0 

0.15 

50. 

33.5 

0.9983 

0.0422 

1.0 

0.14 

50. 

44.7 

0.9991 

0.0636 

1.1 

0.14 

50. 

55.9 

0.9996 

0.0925 

1.1 

0.14 

50. 

61.5 

0.9997 

0.1102 

1.1 

0.13 

SO. 

67.1 

0.9993 

0.1302 

1.1 

0.13 

50. 

72.7 

0.9999 

0.1526 

1.2 

0.13 

50. 

78.3 

0.9999 

0.1774 

1.2 

0.12 

50. 

89.4 

1.0000 

0.2341 

1.3 

0.11 

50. 

100.6 

1 , 0000 

0.2996 

1.4 

0.11 

50, 

111.8 

1,0; 00 

0.3724 

1.6 

0.09 

50. 

139.8 

4  AA.'A 

i  ♦  wvvv 

0.5693 

2.3 

0.06 

50. 

167.7 

1.0000 

0.7500 

4.0 

0.04 

50. 

195.7 

1.CC00 

0.8799 

8.3 

0.02 

50. 

223.6 

1.0000 

0.9530 

21.3 

0.01 

50. 

279.5 

1.0000 

0.9963 

267.2 

0.00 

Table  V-3.  (Continued) 


100. 

0.0 

0.9500 

100. 

15.8 

0,9674 

100. 

31.6 

0.9795 

100. 

47.4 

0.9876 

100. 

63.2 

0.9927 

100. 

79.1 

0.9959 

100. 

87.0 

0.9970 

100. 

94.9 

0.9978 

100. 

102.8 

0.9984 

100. 

110,7 

0.9988 

100. 

126.5 

0.9994 

100. 

142.3 

0.9997 

100. 

158.1 

0.9999 

100. 

197.6 

1.0000 

100. 

237.2 

1.0000 

100. 

276.7 

1.0000 

100. 

316.2 

1.0000 

100. 

395.3 

1.0000 

0.0500 

1.1 

0.0743 

1.1 

0.1067 

1.1 

0.1482 

1.2 

0.1992 

1.3 

0.2597 

1.4 

0.2931 

1.4 

0.3284 

1.5 

0.3653 

1.6 

0.4035 

1.7 

0.4823 

1.9 

0.5618 

2.3 

0.6389 

2.8 

0.803? 

5.1 

0.9124 

11,4 

0.9682 

31.5 

0.9908 

108.2 

0.9996 

2525.9 

0.13 

0.13 

0.13 

0.13 

0.12 

0.11 

0.11 

0.10 

0.0? 

0.09 

0.08 

0.07 

0.05 

0.03 

0.01 

0.00 

0.00 

0.00 


250. 

0.0 

0.8508 

250. 

25.0 

0.8925 

250. 

50.0 

0.9251 

250. 

75.0 

0.9495 

250. 

100.0 

0.9671 

250. 

125.0 

0.9793 

250. 

137.5 

0.9838 

250. 

150.0 

0.9875 

250. 

162.5 

0.9904 

250. 

175.0 

0.9927 

250. 

200,0 

0,9959 

250. 

225.0 

0.9977 

250. 

250.0 

0.9933 

250. 

312.5 

0,9998 

250, 

375.0 

1.0000 

250. 

437.5 

1.0000 

250. 

500.0 

1.0000 

250. 

625.0 

1.0000 

0.1492 

1.4 

0.2005 

1.4 

0.2611 

1.5 

0.3300 

1.6 

0.4052 

1.8 

0.4840 

2.0 

0.523? 

2.2 

0.5636 

2.4 

0.6026 

2.6 

0.6406 

2.8 

0.7123 

3.5 

0.7764 

4,5 

0,8315 

6.0 

0.927? 

13.9 

0,9750 

40.0 

0,9931 

144,0 

0,9985 

650.2 

1,0000 

26715.3 

0.11 

0.10 

0.10 

0.09 

0.08 

0.07 

0.07 

0.00 

0.06 

0.05 

0.04 

0.03 

0.03 

0.01 

0.00 

0.C0 

0.00 

0.00 
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1)  Sediments  change  with  season.  As  shown  in  table  V- 4,  the  14 
August  1981,  elutriates  showed  considerably  higher  release  of 
ammonia  and  COD  then  the  20  November  1981,  samples.  This  is 
likely  due  to  a  lower  redox  potential  of  the  sediments  collected 
during  the  summer  months.  Summer-time  dredging  during  low  flow 
would  represent  worst-case  conditions  in  terms  of  NH^-N  and  COD 
desorption. 

2)  Ammonia  -  nitrogen  was  the  only  chemical  to  clearly  form  a  plume 
in  the  field  sampling  of  20  November.  The  standard  elutriate  test 
"under-predicted"  the  magnitude  of  the  NH^  plume  source  concentra¬ 
tion. 

3)  There  was  considerable  variability  among  sediments.  Results  on  a 
single  sediment  sample  were  quite  reproducible,  but  different 
sediment  samples  from  a  given  vicinity  can  give  significantly 
different  results. 

Source  concentrations  of  the  plume  are  difficult  to  predict 
a  priori.  Also  the  source  width  is  difficult  to  estimate  without  field 
work.  Once  the  source  width  and  strength  are  known,  the  model  developed  in 
this  chapter  is  sufficiently  detailed  to  simulate  the  contaminant  plume. 

The  relationship  between  the  3-D  suspended-solids  model  of  chapter 
II  and  the  2-D  chemical  plume  model  is  readily  apparent.  Both  models,  are 
based  on  steady-flow  onditions  and  continuity  of  mass.  Because  the  3-D 
suspended-solids  model  entered  the  source  input  below  the  surface,  the  mixing 
parameters  in  the  two  model  simulations  were  not  identical.  Greater  lateral 
mixing  and  source  width  were  required  in  the  chemical  plume  model.  It  does 
not  represent  explicitly  the  lateral  mixing  due  to  density  differences  of  the 
plunging  and  subsequently  rising  plume. 
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Table  V-4.  A  Comparison  of  Site  and  elutriate  concentrations  at  RM  170  on 
two  different  dates  and  with  the  model  source  concentrations 


NH,-N 

mg/l 

PO  -P 
mg/l 

COD 

mg/l 

Sediment 

Character 

Site  Water 
170-1 

8-14-81 

0.15 

0.25 

30 

-- 

Elutriate 

170-1 

8-14-81 

0.46 

0.03 

74 

Sand 

Elutriate 

170-2 

8-14-81 

0.90 

0.06 

86 

Organic 

Elutriate 

170-3 

8-14-81 

0.06 

0.03 

73 

Sand 

Site  Water 
170-1* 
11-20-81 

0.03 

0.10 

31 

— 

Elutriate 

170-1* 

11-20-81 

0.08 

0.06 

39 

Sand 

Site  Water 
170-2* 
11-20-81 

0.01 

0.14 

35 

— 

Elutriate 

170-2* 

11-20-81 

0.01 

0.11 

42 

Sand 

Model  Source 
Concentrations 
for  11-20 

0.18 

0.2 

34 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  Conclusions 


1)  Maintenance  dredging  and  disposal  of  main  channel  sediments,  as 
practiced  by  the  Corps  of  Engineers,  was  not  a  significant  water 
quality  problem  in  the  GREAT  III  study. 

2)  A  3-D  suspended-solids  model  and  a  2-D  chemical  contaminant  plume 
model  have  been  applied  to  assess  water  quality  in  the  Great  III 
reach  of  the  Mississippi  River.  Refinements  were  made  to  the 
models.  A  variable-depth  source  option  was  included  in  the 
suspended-solids  model  to  reflect  conditions  observed  in  the  GREAT 
III  study  area.  Sorption  partitioning  and  chemical  kinetics  of 
contaminant  transformations  were  Incorporated  into  the  chemical 
plume  model.  Model  reaction  and  transport  formulations  were 
sufficient  to  describe  the  predominant  processes  at  work  in  the 
Mississippi  River. 

3)  Limitations  of  the  models  Include: 

a)  The  steady-state  assumption.  Due  to  river  dynamics  and  time- 
variable  discharges  of  dredged  material,  the  assumption  of 
steady  state  was  not  always  applicable.  Dredge  disposal  was 
halted  periodically  to  allow  barge  traffic  to  pass.  Also 
solids  concentrations  changed  abruptly  as  the  cutterhead 
moved. 

b)  The  open-field  assumption.  Wing  dams  and  variable  velocities 
and  depths  caused  spatial  structure  in  the  field  data  which  is 
not  accounted  for  by  the  models. 

c)  The  source  configuration.  It  is  not  presently  possible  to 
predict  a_  priori  the  dredge  disposal  discharge 
characteristics,  i.e.,  the  source  width,  the  source  depth,  and 


1 


the  source  strength  (concentration).  These  parameters  can  be 
estimated  from  elutriate  tests  and  historical  field  studies  to 
within  an  order-of-magnitude,  but  it  is  necessary  to  perform 
studies  at  the  site  to  make  more  accurate  simulations. 

d)  The  lateral  (transverse)  mixing  coefficient.  Lateral  mixing 
in  the  GREAT  III  field  study  was  considerably  greater  than 
lateral  dispersion  coefficients  reported  in  the  literature. 
Increased  lateral  mixing  was  caused  by  the  lateral  velocity  of 
the  discharge  jet,  the  density-driven  sinking  and  rise  of  the 
plume,  and  by  the  irregular  bed  geometry.  An  apparent  lateral 
mixing  coefficient  of  2.0  m2/sec  was  determined  by  simulations 
of  the  suspended-solids  field  data. 

Model  verification  was  complicated  by  the  low  suspended  solids  and 
contaminant  concentrations  which  were  measured  in  this  study.  In- 
si  tu  concentrations  were  several  times  lower  than  in  the  GREAT  II 
study  reach,  from  Guttenberg,  Iowa,  to  Saverton,  Missouri.  Lower 
concentrations  were  caused  by  greater  depth  and  flow-per-unit- 
width,  as  well  as  coarser  and  "cleaner"  sand  than  in  the  GREAT  II 
reach. 

Water-quality  criteria  were  not  exceeded  in  the  field  study  and 
were  only  rarely  exceeded  in  the  standard  elutriate  test. 
Ammonia-nitrogen,  chemical  oxygen  demand,  and  manganese  were  the 
most  frequently  desorbed  chemicals.  Dissolved  iron,  phosphate, 
and  zinc  were  the  most  frequently  adsorbed.  Heavy  metals  were  not 
significantly  desorbed. 

Sediments  were  relatively  free  of  contaminants  with  the  possible 
exception  at  RM  172  below  the  Des  Peres  River  at  St.  Louis, 
Missouri.  (Sediment  sample  172-3  contained  ~  0.4  mg/kg  (dry 
basis)  of  PCBs.  Sample  172-1  contained  25.3  y g/1  of  the 
plasticizer  diethyl  phthalate  and  11.1  u 9/1  of  the  plasticizer 
dilrityl  phthalate.) 


7)  The  standard  elutriate  test  was  useful  as  a  screening  test.  It 

likely  provides  an  order-of-magnitude  estimate  of  dredge-disposal 
discharge  concentrations.  It  varies  from  the  discharge 

principally  due  to  differences  in  redox  potential,  sediment-to- 

water  ratio,  and  mixing  time/intensity. 

8)  The  kinetics  of  the  elutriate  test  could  be  described  by  an 

initial  condition  and  three  stages.  The  initial  condition  depends 
on  the  redox  of  the  sediment  sample  and  the  related  parameters: 
volatile  solids  content,  percent  moisture,  and  nature  of  the 
sediment  (grain  size,  color,  organic  content).  The  first  stage  is 
described  by  an  Initial  mixing  and  desorption  which  takes  place 

over  a  few  minutes.  The  second  stage  involves 

adsorption/desorption,  ion  exchange,  and  nucleation/precipitation 
processes  which  takes  from  30-90  minutes.  The  final  stage  is 
sorptive  equilibria.  Biotransformations,  volatilization,  and  some 
chemical  oxidations  may  occur  over  a  longer  time  scale.  A 

mathematical  model  was  developed  to  account  for  the  kinetics  of 
the  standard  elutriate  test. 

B.  Reco— endatlons 

1)  Further  verification  of  the  3-D  suspended-solids  model  and  the  2-D 
chemical  plume  model  is  not  necessary  for  the  reasons  stated  in 
conclusions  1,3,  and  4.  The  models  can  be  used  to  estimate  in- 
situ  concentrations  within  an  order-of-magnitude  (or  better  if 
field  measurements  are  taken  near  the  dredge-disposal  discharge). 

2)  However,  If  further  verification  studies  of  the  model  were  to  be 
undertaken,  a  better  knowledge  of  ambient  variations  in  suspended 
solids  and  chemical  concentrations  should  be  obtained,  in  both  the 
vertical  and  lateral  dimensions. 

3)  If  care  is  taken  In  the  collection  of  sediments  and  preparation, 
the  standard  ellutriate  test  can  be  a  useful  screening  tool  to 
determine  potential  problem  cases,  especially  for  chemical  oxygen 


demand  and  ammonia  which  are  the  oxygen  demanding  pollutants.  To 
be  more  representative  of  disposal  conditions,  it  is  recommended 
that  the  ratio  of  sediment-to-water  be  changed  from  1:4  to  1:9  in 
the  elutriate  test. 

4)  The  chemical  plume  model  can  be  used  to  roughly  estimate 
contaminant  concentrations  before  dredging  by  the  following 
methodology.  See  also  the  GREAT  II  report,  "Suspended  Sediment 
Modelling  of  Dredge  Disposal  Effluent  in  the  GREAT  II  Study 
Reach",  by  Schnoor  et  al  (1980). 

a)  Measure  or  estimate  the  stream  velocity.  Range  0.4-1. 5  m/sec 
during  late  summer. 

b)  Measure  the  size  distribution  of  the  sediment  to  be  dredged 
and  the  weight  fraction  of  the  sand,  silt,  and  clay  particle 
classification.  This  can  be  performed  by  visual  accumulation 
(VA)  tube  and  micropipette  analyses.  Sand  fall  velocities 
normally  range  from  0.007  -  0.015  m/sec.  For  medium  grade 
sand  sediment  (Dgg  =  370  y),  a  direct  runoff  composition  of 
45%  sand,  25%  silt,  and  30%  clay  is  suggested.  If  the  dredge 
discharge  water  is  expected  to  run  across  an  island  before 
entering  the  river,  a  composition  of  2%  sand,  67%  silt,  and 
30%  clay  Is  suggested.  See  fall  velocities  in  chapter  II. 

c)  Estimate  the  source  concentration  of  suspended  solids  at  the 
point  where  the  discharge  water  reenters  the  river.  The  range 
is  from  20  -  200  mg/1  suspended  solids  with  a  suggested  rule 
of  thumb:  TSS  (mg/1)  =  50./u(m/sec). 

d)  Estimate  the  source  concentration  of  chemical  contaminant  from 
the  standard  elutriate  test  or  from  the  relationships  found  in 
chapter  III.  Determine  the  partition  coefficient  between 
sand,  silt,  clay,  and  dissolved  fractions  by  wet  sieving  and 
measurement  of  the  particulate  and  dissolved  phases. 


157 


e)  Estimate  the  width  f  the  discharge  based  on  the  method  of 
disposal  (i.e.,  side  bank,  near  shore,  or  open  water).  A 
suggested  average  width  Is  20  m  with  a  range  of  10-120  m. 

f)  Assume  a  lateral  (transverse)  dispersion  coefficient  of  2.0 
m2/sec.  Range:  0.1  -  20  n^/sec. 

g)  From  the  literature  determine  If  there  are  any  rapid 
transformations  (e.g.,  biodegradation,  volatilization, 
photodegradation)  which  could  occur  In  the  near  field.  Most 
chemicals  do  not  react  sufficiently  fast  to  require  the 
specification  of  k,  the  first-order  degradation  rate  constant. 

h)  Use  the  chemical  plume  model  as  described  in  Chapter  5. 


APPENDIX  A 

SUSPENDED-SOLIDS  COMPUTER  PROGRAM 
USER'S  MANUAL 
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A.  INTRODUCTION 

The  suspended-solids  computer  program  is  basically  the  WALDEN  PLUME 
MODEL  (WP)  presented  by  Wechsler  and  Cogley  (1977)  with  revisions. 

The  original  model,  WP,  was  used  to  investigate  the  turbidity 
generated  by  dredging  and  open-water  disposal  operations.  The  study  was 
conducted  to  determine  the  nature  and  amount  of  turbidity  that  may  be  produced 
by  a  given  sediment  when  it  is  suspended  by  dredging  operations  and  to 
evaluate  the  effects  of  sediment  and  water  compositional  factors  which 
influence  particle  settling  rates. 

Though  WP  was  developed  for  estuarine  conditions,  Schnoor  et  al . 
(1980),  modified  it  and  used  It  to  obtain  predictions  of  GREAT  II  field  data. 

The  present  study  modified  the  Wechsler-Cogley  model  further,  taking 
the  GREAT  II  version  of  It  as  a  starting  point.  Once  relatively  satisfactory 
results  were  obtained  for  two  GREAT  II  and  two  GREAT  III  sites,  the  refined 
model  was  used  to  predict  suspended-solids  distributions  for  a  third  GREAT  III 
site.  The  suspended-solids  computer  model  is  fully  discussed  in  Chapter  II. 

The  model  revisions  include: 

1) .  the  ability  to  model  three  modes  of  disposal— open  water,  intermediate  or 

near-shore,  and  on-shore  disposal; 

2) .  the  folding  or  reflection  principle  needed  for  near-shore  and  on-shore 

disposal  has  been  included  in  the  program; 

3) .  two  sediment  division  options  are  possible—fines/sand  or  silt/clay/sand; 

4) .  three  plume-source  options  are  now  possible-constant  concentration  from 

the  surface  to  the  bottom,  constant  concentration  over  a  fraction  of  the 
depth,  or  triangular  concentration  from  the  surface  to  the  bottom; 

5) .  two  lateral  spreading  options  are  possible— constant  diffusion 

coefficient,  Ez=  Constant  or  decreasing  diffusion  coefficient, 

Ez=  f  (spreading  factor). 

B.  ASSUMPTIONS  The  assumptions  made  in  the  derivation  of  WP  which  also  apply 
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to  the  present  model  include: 

1.  unidirectional,  steady  flow  of  water, 

2.  infinite  width  compared  to  depth, 

3.  constant  depth, 

4.  infinite  length  compared  to  depth, 

5.  eddy  diffusion  in  downstream  direction  is  negligible  compared  to  other 
diffusion  and  transport  terms, 

6.  eddy  diffution  in  vertical  direction  can  be  related  to  vertical 
position  in  the  flow  by  E^=  0.02  Uy 

7.  eddy  diffusion  in  lateral  direction  in  E  ~  2.2  (E„l  „ 

y  y  max 

Any  other  assumptions  which  apply  to  the  present  model  are  listed  and 
discussed  in  chapter  II. 

C.  SUSPENDED  SOLIDS  MODEL  INPUT  VARIABLES 

All  input  data  are  entered  as  format-free  variables. 


Variable 

Definition 

ISCE 

type  (configuration)  of  sediment  source 

IUP 

upper  limit  of  suspended-solids  source 

IBT 

lower  limit  of  suspended-solids  source 

ITYPE 

disposal  type 

NBANK 

distance  between  river  bank  and  source  location 

CONOR 

original  concentration  (mg/1) 

EZ 

lateral  eddy  diffusivity  (m2/sec) 

NSEDF 

number  of  sediment  fractions 

U 

average  stream  velocity  (m/sec) 

H 

average  stream  depth  (m) 

NSTEP 

number  of  downstream  steps 

XL 

initial  source  half-width  (m) 

DELZ 

lateral  step  size  (m) 

EXPCO 

spreading  factor 

W 

particle  settling  velocity  (m/sec) 

CO 

sediment  fraction  (%) 
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0.  SUSPENDED-SOUPS  MODEL  OUTPUT 

1.  Instructions  describing  input  variables 

2.  input  values  of  concentration  distributions  and  fall  velocities 

3.  vertical  distribution  of  sediment  downstream  for  each  sediment  fraction 
(without  lateral  spreading) 

4.  the  summation  of  the  vertical  distributions  for  all  size  fractions 

5.  the  lateral  spreading  coefficients 

6.  horizontal  slices  through  the  3-dimensional  plume  at  five  preselected 
depths  showing  the  concentration  distribution  at  each  depth 

7.  folded  distributions  of  (6.)  if  reflection  is  needed  for  disposal  type 

The  interpretation  and  application  of  this  model  are  fully  explained 
in  chapter  II. 
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9  SUSPENDED  SOLIDS  COMPUTER  PROGRAM  LISTING 


C 

c 

C  PROGRAM  SOURCE:  "A  LABORATORY  STUDY  OF  THE  TURBIDITY 
C  GENERATION  POTENTIAL  OF  SEDIMENTS  TO  BE  DREDGED,''  BY 
C  B.A.  WECHSLER  &  D.R.  COGLEY,  TECH.  REPORT  D-77-14, 

C  U.S.  ARMY  ENGINEEER  WATERWAYS  EXPERIMENT  STATION, 

C  NOV.  1977 

C 

C 

C  WALDEN  PLUME  MODEL 

Q  ****************************************************** 

C  ASSUMPTIONS  INCLUDE  STEADY  STATE,  NO  MOMENTUM  EFFECTS 
C  AND  NO  RESUSPENSION  OF  MATERIAL 
C 

C  ****************************************************** 

C  THIS  PROGRAM  CONTAINS  A  SUBROUTINE  FROM  THE  IMSL 
C  LIBRARY,  A  PROPRIETARY  PACKAGE  FROM  THE  INTERNATIONAL 
C  MATHEMATICAL  &  STATISTICAL  LIBRARIES,  INC.,  HOUSTON, 

C  TEXAS.  THIS  ROUTINE  MAY  NOT  BE  REDISTRIBUTED  OR 

C  REMOVED  FROM  THIS  SOFTWARE  FOR  USE  IN  OTHER 

C  SOFTWARE  DEVELOPMENT.  THE  IMSL  ROUTINE  INCLUDED  IS 

C  "ERF." 

Q  ****************************************************** 

c 

C  THE  PRESENT  VERSION  OF  THE  PROGRAM  IS  BASED  ON  A 
C  MODIFICATION  OF  THE  ORIGINAL  WECHSLER- COGLEY  MODEL 
C  FROM  "SUSPENDED  SEDIMENT  MODELING  OF  DREDGE-DISPOSAL 

C  EFFLUENT  IN  THE  GREAT  II  STUDY  REACH,"  BY  J.L. 

C  SCHNOOR  ET  AL.,  IIHR  REPORT  NO.  228,  MARCH  1980. 

C 

q  ****************************************************** 

C 

C  AS  PART  OF  THE  PRESENT  PROJECT,  "REFINEMENT  AND 

C  VERIFICATION  OF  PREDICTIVE  MODELS  OF  SUSPENDED 

C  SEDIMENT  DISPERSION  AND  DESORPTION  OF  TOXICS  FROM 

C  DREDGED  SEDIMENTS,"  DONE  FOR  THE  CORPS  OF  ENGINEERS  - 

C  GREAT  III,  A  NUMBER  OF  MODIFICATIONS  HAVE  BEEN 

C  INCORPORATED: 

C 
C 


1.  THREE  MODES  OF  DISPOSAL: 


A.  OPEN  WATER 

B.  INTERMEDIATE  OR  NEAR-SHORE 

C.  ON-SHORE 


2.  "FOLDING”  FOR  INTERMEDIATE  AND  BANK  DISPOSAL 
HAS  BEEN  INCLUDED  AS  A  PART  OF  THE  PROGRAM. 

3.  TWO  SEDIMENT-DIVISION  OPTIONS  ARE  INCLUDED: 

A.  FINES/SAND 

B.  SILT/CLAY/SAND 

4.  THREE  PLANE-SOURCE  OPTIONS  ARE  INCLUDED: 

A.  CONSTANT  CONCENTRATION,  SURFACE  TO  BOTTOM 

B.  CONSTANT  CONCENTRATION,  FRACTION  OF  DEPTH 

C.  TRIANGULAR  CONCENTRATION,  SURFACE  (C=0. )  TO 
BOTTOM  (C=CONOR) 

5.  TWO  LATERAL  SPREADING  OPTIONS  ARE  INCLUDED: 

A.  CONSTANT  E2  (EXPCO=0) 

B.  DECREASING  E2  =  F(EXPCO) 

******************4******************************** 

REAL  IA, MA, ML, I IA 
DIMENSION  WM(21, 21 ) 

DIMENSION  FOLD1 (30,51), FOLD2 (30, 51 ) , F0LD3 (30, 51 ) 
DIMENSION  C(20,51),D(20,20), IA(20, 20) ,MA(20, 20) 
DIMENSION  RA(20, 20) ,CJ(20),IIA(20,20),B(20) 
DIMENSION  ADELX(51 ) , ADELY(20 ) , Z( 21, 51 ) , AVG(51 ) 
DIMENSION  AVGSUM( 51 ) , B8( 20 ) , IOUT( 51 ) , CSUM( 20, 51 ) 
DATA  SILT, CLAY, SAND,FINE/' SILT* , 'CLAY' , 'SAND' , 'FINE 
DATA  BD1 , 0WD1 , 0WD2 , 0WD3/' BANK' , ' OPEN' , '  WAT' , ' ER' / 
DATA  BD2 , BD3 , BD4/ ' INTE ' , ’ RMED ' , ’ I ATE ' / 

DATA  GDI , GD2 , GD3/'  DIS' , 'POSA' , ’ L'/ 

E(Y)=0 . 02*U*Y* ( 1-Y/H) 

CALL  OPENSA( ASWRIT, 'STODAT' ,6,4) 

REWIND  8 

ENTER  THE  SOURCE  TYPE: 

1  =  RECTANGULAR,  SURFACE  TO  BOTTOM 

2  =  RECTANGULAR,  FRACTION  OF  DEPTH 

3  =  TRIANGULAR,  SURFACE  (C=0)  TO  BOTTOM  (C=CONOR) 

WRITE(1, 18) 

18  FORMAT (2X, 'ENTER  SOURCE  TYPE'/2X, 

+'1=RECT  -  SURF  TO  B0T'/2X, ' 2=RECT  -  FRAC  OF  DEPTH'/ 
+2X, ' 3=TRIANG' ) 

READ ( 1 , * ) I SCE 
I F ( I SCE . EQ . 1 ) GO  TO  22 
I F ( I SCE . EQ . 2 ) GO  TO  24 


ooo  oonoonoooo  ooon  ooon 


r 
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I F ( I SCE . EQ . 3 ) GO  TO  23 
22  IUP=1 
IBT=20 
GO  TO  23 
24  WRITE ( 1, 17) 

ENTER  THE  UPPER  LIMIT  OF  THE  SUSPENDED  SOLIDS  SOURCE 
(IUP)  IN  TENTHS  (OF  DEPTH)  FROM  THE  BOTTOM  UP 

17  FORMAT(2X/ 'ENTER  THE  UPPER  LIMIT  OF  SOURCE,  IUP, ' , 

♦  '  IN  TENTHS  (OF  DEPTH)  FROM  THE  BOTTOM  UP') 

READ  ( 1 ,  *  )  I  UP 
IUP=20-2* IUP+1 


ENTER  THE  LOWER  LIMIT  OF  THE  SUSPENDED  SOLIDS  SOURCE 
( IBT)  IN  TENTHS  (OF  DEPTH)  FROM  THE  BOTTOM  UP 

WRITE(1, 19) 

19  FORMAT ( 2 X, 'ENTER  THE  LOWER  LIMIT  OF  SOURCE, IUP, ' , 

+'  IN  TENTHS  (OF  DEPTH)  FROM  THE  BOTTOM  UP') 

READ ( 1 , * ) I BT 
IBT=20-2*IBT 

ENTER  THE  TYPE  OF  DREDGE  SPOIL  DISPOSAL: 

1  =  BANK  DISPOSAL 

2  =  OPEN  WATER  DISPOSAL 

3  =  INTERMEDIATE  DISPOSAL 

NBANK  =  NUMBER  OF  UNITS  OF  DELZ  WHICH  YOU  WANT  TO 
PUT  IN  BETWEEN  THE  BANK  AND  THE  START  OF  THE  ORIGINAL 
CONCENTRATION  ("FOLDING"  IS  DONE  AT  THE  BANK) 
DELZ=LATERAL  STEP  SIZE 

23  WRITE (1,15) 

15  FORMAT (2X, 'ENTER  DISPOSAL  TYPE ' /2X, ' 1=BANK' /2X, 
+'2=OPEN  WATER '/2X, ' 3= INTERMEDIATE ' ) 

READ ( 1 , * ) I TYPE 

IF( ITYPE . EQ . 3 ) GO  TO  20 

GO  TO  21 

20  WRITE (1,16) 

16  FORMAT(2X, 'ENTER  THE  UNITS  OF  DELZ  YOU  WANT  TO  PUT', 
+'  BETWEEN  THE  BANK  ("FOLDING"  LINE)  AND  THE  START', 

+'  OF  THE  ORIGINAL  CONCENTRATION') 

RE AD (1,*) NBANK 

21  IF( ITYPE. EQ. 1)NBANK=0 

IF( ITYPE. EQ. 1) WRITE (8,30) BD1, GDI, GD2 , GD3 
30  FORMAT ( ////2 X, 6A4) 

IF( ITYPE . EQ . 2 )WRITE( 8 , 30 )0WD1 , OWD2 , OWD3 , GDI , GD2 , GD3 
IF( ITYPE . EQ . 3 ) WRITE ( 8 , 30 )BD2 , BD3 , BD4, GDI , GD2 , GD3 

INITIALIZE  ARRAYS 


J 


o  o  o  o 
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CALL  MCON  (0,20,20,0. ) 

CALL  MCON(CSUM,20, 51,0. ) 

CALL  MCON  (Z,21, 51,0. ) 

CALL  MCON  (F0LD1, 30,51,0. ) 

CALL  MCON  ( F0LD2 ,30,51,0. ) 

CALL  MCON  ( F0LD3 ,30,51,0. ) 

CALL  MCON  (AVGSUM,51, 1,0. ) 

CALL  MIDN  ( I A, 20 ) 

CALL  MIDN  ( 1 1 A, 20 ) 

CALL  MSCL  (IIA, IIA, 20, 20,-1. ) 

DO  10  1=1,20 

IF( I . GT . 1 ) I IA( I , 1-1 )=1 . 

10  CONTINUE 
C 

C  ENTER  MODEL  PARAMETERS: 

C 

C  CONOR=OR I G I NAL  CONCENTRATION  (MG/L) 

C  EZ=LATERAL  EDDY  DIFFUSIVITY  (M**2/SEC) 

C  NSEDF=NUMBER  OF  SEDIMENT  FRACTIONS 

C  U=AVERAGE  STREAM  VELOCITY  (M/SEC) 

C  H=  "  "  DEPTH  (M) 

C  NSTEP=NUMBER  OF  DOWNSTREAM  STEPS  FOR  F.D.  SOLUTION 
C  XL=INITIAL  SOURCE  HALF-WIDTH  (M) 

C  DELZ=LATERAL  STEP  SIZE  (M) 

C  EXPCO=SPREADING  FACTOR 

C 

WRITE( 1,13) 

13  FORMAT (  2X, ' ENTER  ORIGINAL  CONCENTRATION  (MG/L)') 
READ(1, *) CONOR 

NSED=1 
WRITE( 1, 14) 

READ ( 1 , * ) EZ 

14  FORMAT (2X, 'ENTER  E-SUB-Z  (LATERAL  EDDY  DIFFUSIVITY)') 
WR I TE (1,9) 

9  FORMAT (2X, 'ENTER  NSEDF' ) 

READ ( 1 , * ) NSEDF 
WRITE(1, 1000) 

1000  FORMAT ( 2X, ' ENTER  VALUES  OF  U, H, NSTEP , XL, DELZ , EXPCO' ) 
READ ( 1 , * )  U,H, NSTEP, XL, DELZ, EXPCO 
WR I TE ( 8 , 1010)U,H, XL, CONOR, EZ 
1010  FORMAT ( //2X, ' INPUT  PARAMETERS //2X, 

+' STREAM  VELOCITY  (M/SEC):  U  =',F6.2/2X, 

+ ' STREAM  DEPTH  (M) :  H  =',F5.1,/2X, 

+' SOURCE  HALF-WIDTH  (M) :  XL  =', F5 . 1 , /2X, ' ORIGINAL ' , 

+'  TOTAL  SEDIMENT  CONCENTRATION  (MG/L):  CONOR  =' , 
+F5.0,/2X, ' LATERAL  EDDY  DIFFUSIVITY  (M**2/SEC):  EZ  =', 
♦F8.6,////) 

PROGRAM  CONTROL  IS  TRANSFERRED  HERE  FOR  EACH 
SEDIMENT  FRACTION 
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11  IF(NSED  .GT.  NSEDF)GO  TO  99 

12  FORMAT (2X,A4) 

IF(NSEDF . EQ. 3 . AND . NSED .  EQ.  1 )WRITE( 1 , 12 ) SILT 
I F ( NSEDF . EQ . 3 . AND . NSED . EQ . 2 ) WRI TE ( 1 , 12 ) CLAY 
I F ( NSEDF . EQ . 3 . AND . NSED . EQ . 3 ) WRITE ( 1 , 12 ) SAND 
IF(NSEDF. EQ. 2 . AND . NSED . EQ . 1 )WRITE( 1 , 12 ) FINE 
IF (NSEDF . EQ. 2 . AND .NSED . EQ. 2 )WRITE( 1 , 12 ) SAND 
C 

C  ENTER  INDIVIDUAL  SEDIMENT-FRACTION  PARAMETERS: 

C  CO=SEDIMENT  FRACTION  (%) 

C  (SUM  OF  ALL  C0=100%) 

C  W=P ARTICLE  SETTLING  VELOCITY  (M/SEC) 

C 

WRITE(1, 101) 

101  FORMAT (2X, 'ENTER  VALUES  OF') 

WRITE (1, 102) 

102  FORMAT (2X, 'W,CO' ) 

READ(1, *)W,CO 

IF ( NSEDF. EQ. 3 . AND . NSED . EQ . 1 )WRITE( 8, 1011 ) SILT, W, CO 
IF ( NSEDF . EQ . 3 . AND . NSED . EQ . 2 ) WRITE ( 8 , 1011 ) CLAY , W , CO 
IF(NSEDF. EQ. 3 .AND .NSED. EQ . 3 )WRITE(8, 1011 ) SAND, W, CO 
IF(NSEDF.EQ.2 . AND . NSED . EQ . 1)WRITE(8, 1011 ) FINE, W, CO 
I F ( NSEDF . EQ . 2 . AND . NSED . EQ . 2 ) WR I TE ( 8 , 1011 ) SAND, W, CO 
1011  FORMAT (////2X, 'VERTICAL  CONCENTRATION  DISTRIBUTION', 
+'  OF ' , A4 , '  (%) '/2X, 'FALL  VELOCITY  (M/SEC):  W=',F8.6, 
+/2X, 'ORIGINAL  CONCENTRATION  (%) :  CO  =',F5.0//) 

N=20 

THETA=1 . 0 

IF  (NSTEP.GT. 50)  NSTEP=50 

NSTEP1=NSTEP+1 

DELX=10 . 

DELY=H/FLOAT(N) 

DELI 2=1 . / ( DELY*DELY ) 

AVGSUM ( 1 ) =AVGSUM ( 1 ) +C0 
DO  501  J=1 , NSTEP1 
ADELX ( J ) = ( J- 1 ) *DELX 
AVG ( J ) =0 . 

501  CONTINUE 

DO  500  1=1, N 

ADELY ( I )=( 1-1 ) *DELY 

IF  ( (I-1)*DELZ.LE.XL)  Z(I,1)=1. 

500  CONTINUE 
C 

C  CALCULATE  CONCENTRATION  DUE  TO  LATERAL  DISPERSION 
C 

DO  680  J=1 , NSTEP 

FOREX=SQRT (ADELX(J+1)*4.*(EZ*( 1 +EXPCO*  EXP ( - ( 

+ADELX( J+l )/100 . ) ) ) ) 

DO  680  IZ=1 , 21 
AZ=( IZ-1 ) *DELZ 
TOP= ( AZ+XL ) /FOREX 
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ET=ERF ( TOP ) 

BOT=(AZ-XL) /FOREX 
EB=ERF (BOT) 

Z( IZ, J+l )=0 . 5* (ET-EB) 

680  CONTINUE 
C 

C  START  FINITE  DIFFERENCE  SOLUTION 
C 

I F ( I SCE . EQ . 3 ) GO  TO  25 
DO  100  I=IUP,20 
C(I,l)=CO 
100  CONTINUE 
GO  TO  26 

25  TRI=0. 

DO  103  1=1,20 
C(I, l)=CO*TRI 
TRI=TRI + . 05 
103  CONTINUE 

26  AVG ( 1 ) =CO 
D(l/1)=-E(l. 5*DELY ) 

D  ( 1 , 2 ) =E (1.5  *DELY ) 

N1=N-1 

C 

C  CALCULATE  EDDY  DIFFUSIVITY  AT  VARIOUS  DEPTHS 
C 

DO  200  1=2, N1 
X1=(2*I-1 ) * . 5*DELY 
X2=( 2*1+1)*. 5*DELY 
D( I , 1-1 )=E(X1 ) 

D( J, 1)=-E(X1)-E(X2) 

D ( I , 1+1 )=E(X2 ) 

200  CONTINUE 

D(N,N)=-E( ( FLOAT (N) - . 5 ) *DELY) 

D(N, N- 1 ) =-D( N, N ) 

DO  300  1=1, N 
DC  300  J=1 , N 
D( I , J )=DELI2*D( I, J) 

C 

C  CALCULATE  COEFFICIENTS  OF  SYSTEM  OF  EQUATIONS 
C  RESULTING  FROM  FINITE  DIFFERENCE  SOLUTION 
C 

MA( I , J )=U/DELX* I A( I , J ) -THETA*D( I , J) 

RA( I , J ) =U/DELX* I A ( I , J) +( 1 . -THETA ) *D( I , J) +W* I IA( I , J) 
+/DELY 

300  CONTINUE 

DO  400  J=1 , NSTEP 
DO  405  1=1, N 
CJ( I )=C( I, J) 

CSUM( I , J ) =CSUM( I , J ) +CJ ( I ) 

405  CONTINUE 
C 
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C  FINAL  SOLUTION  BY  COMBINATION  OF  ANALYTIC  AND 
C  NUMERICAL  SOLUTIONS 

C 

CALL  MMLT(B,RA,CJ,N,N, 1) 

CALL  LINEQ(B8,B,MA,WM,20,21,  IER) 

IF  (IER.GT.O)  WRITE ( 1,4747) 

IF  (IER.GT.O)  WRITE ( 8 , 4747 ) 

4747  FORMAT ( 1  NO  SOLUTION  FOUND') 

C 

C  OUTPUT  OF  FINAL  RESULTS,  FOLDING 

C 

DO  399  K=1 , N 
C ( K , J  + 1 ) =B8 ( K ) 

AVG( J+1)=B8(K)+AVG( J+l) 

399  CONTINUE 

AVG( J+1)=AVG( J+l)  /FLOAT (N) 

AVGSUM ( J+ 1 ) =AVGSUM ( J+ 1 ) + AVG ( J+ 1 ) 

400  CONTINUE 

WRITE(8,5001) (ADELX(I), 1=1 , NSTEP1 , 5 ) 

WRITE (8, 5010) 

5010  FORMAT ( / ) 

5001  FORMAT (10X,20F6.0) 

DO  410  1=1 ,N 

CSUM ( I , NSTEP 1 ) =CSUM ( I , NSTEP 1 ) +C( I , NSTEP 1 ) 

WRITE ( 8 , 5000 )  ADELY(I), (C ( I , J) , J=1 , NSTEP1 , 5 ) 

5000  FORMAT (3X, F5 . 2 , 2X, 20F6 . 2 ) 

410  CONTINUE 

WRITE (8, 5003)  ( AVG(KK) , KK=1 , NSTEP1 , 5 ) 

5003  FORMAT (4X, 'AVG' ,3X,20F6. 2) 

NSED=NSED+1 
GO  TO  11 
99  CONTINUE 

WRITE ( 8 , 5005 ) 

5005  FORMAT  ( 1H1 , //T25 SUMMATION  OF  SUSPENDED  SEDIMENT', 

+'  (%)'//) 

WRITE (8, 5001 ) (ADELX( I ) , 1=1 , NSTEP1 , 5 ) 

WRITE(8, 5010) 

DO  412  1=1, N 

WRITE(8, 5000)  ADELY(I), (CSUM( I , J) , J=l, NSTEP 1, 5 ) 

412  CONTINUE 

WRITE(8, 5003 )  (AVGSUM(KK) , KK=1 , NSTEP1 , 5) 

WRITE(8, 6668) 

6668  FORMAT (1H1,//T2 5, 'LATERAL  SPREADING  COEFFICIENTS'//) 
WRITE( 8, 5006) ( ADELX(KK) , KK=1 , NSTEP 1, 5 ) 

5006  FORMAT ( 1 IX , 2  0F6 . 0 ) 

WRITE(8, 5010) 

DO  411  1=1,21 
AZ=( I-1)*DELZ 

WRITE (8, 5002)  AZ, ( Z ( I , J) , J=1 , NSTEP 1 , 5 ) 

5002  FORMAT ( 4X, F5 . 1 , 2X, 20F6 . 3 ) 

411  CONTINUE 
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DO  800  IY=1,20,4 
YVAL=( IY-1 ) *DELY 

IF( ITYPE . EQ. 1 . OR. ITYPE . EQ. 3 )GO  TO  27 
WRITE (8, 5555)  YVAL 

5555  FORMAT (1H1,////,T7, 'DISTRIBUTION  OF  SEDIMENT  IN', 

+'  HORIZONTAL  PLANE  AT  DEPTH' , F6 . 2 , '  M  (MG/L)',//) 
WRITE ( 8 , 5001 ) ( ADELX( KK ) , KK=1 , NSTEP1 , 5 ) 

WRITE( 8, 5010) 

27  IF1=0 
IF2=0 

DO  810  IZ=1,21 
IAZ=22- IZ 

AZ= ( -DELZ* ( I AZ- 1 ) ) 

IF1=IF1+1 

DO  801  IX=1 , NSTEP1 

IOUT( IX)=( (CSUM( IY, IX) *Z ( IAZ, IX) *1000 . +0 . 5 )/800 . ) 

+* (CONOR/125. ) 

801  CONTINUE 

IF( ITYPE. EQ.l. OR. ITYPE. EQ. 3 )GO  TO  28 
WRITE (8, 6666)  AZ, ( IOUT(KK) , KK=1 , NSTEP1 , 5 ) 

6666  FORMAT (2X, F6 . 1 , 2X, 20F6 . 0) 

28  I F ( ITYPE . EQ . 2 ) GO  TO  810 
IXLDZ=2 1- I NT { XL/DELZ ) -NBANK 
IF(IF1.GE. IXLDZ)GO  TO  860 
DO  850  IFL1=1 , NSTEP1 , 5 
FOLD1 ( IF1, IFL1 )=IOUT( IFL1 ) 

850  CONTINUE 
GO  TO  810 

860  IF2=IF2+1 

DO  852  IFL2=1 , NSTEP1 , 5 
FOLD2 ( IF2 , IFL2 )=IOUT( IFL2 ) 

852  CONTINUE 
810  CONTINUE 

DO  802  IZ=2 , 21 
AZ=DELZ*( IZ-1) 

IF2=IF2+1 

DO  803  IX=1 , NSTEP1 

IOUT( IX)  =  ( (CSUM( IY, IX)*Z(IZ, IX ) *  1000 . +0 . 5 )/800 . ) 

+* (CONOR/125. ) 

803  CONTINUE 

IF( ITYPE. EQ.l. OR. ITYPE. EQ. 3) GO  TO  29 
WRITE (8, 6666)  AZ, ( IOUT(KK) , KK=1 , NSTEP1 , 5 ) 

29  I F ( I TYPE . EQ . 2 ) GO  TO  802 
DO  851  IFL2=1 , NSTEP1 , 5 
F0LD2 ( IF2 , IFL2 )=IOUT( IFL2 ) 

851  CONTINUE 

802  CONTINUE 

I F ( I TYPE . EQ . 2 ) GO  TO  800 
DO  853  1=2 , IXLDZ 
IT=IXLDZ+1- I 
DO  853  J=1 , NSTEP1 , 5 
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F0LD3 ( I , J)=F0LD1 ( IT, J) 

853  CONTINUE 

DO  855  1=1 , NSTEP1 , 5 
FOLD3 ( 1 , I )=0 . 

855  CONTINUE 

DO  857  1=2 , NSTEP1 
FOLD2 (1, I )=2 . *FOLD2 (1,1) 

857  CONTINUE 

CALL  MADD  ( FOLD1 , FOLD2 , FOLD3 ,30,51) 

AZ=0 . 

5777  FORMAT(  1H1 ,  ////,  T7 ,  '  ''FOLDED”  DISTRIBUTION  AT'  ,  F6. 2, 

+'  M  (MG/L)'//) 

WRITE (8,5777 ) YVAL 

WRITE(8, 5001 ) ( ADELX(KK) , KK=1, NSTEP1 , 5 ) 

WRITE ( 8 , 5010 ) 

DO  854  1=1,30 

WRITE ( 8 , 6666 )AZ , ( FOLD1 ( I , K) , K=1 , NSTEP1 , 5 ) 

AZ=AZ+DELZ 

854  CONTINUE 
800  CONTINUE 

CALL  CLOA ( 4 ) 

STOP 

END 

****************************************************** 

IMSL  ROUTINE  NAME  -  MERF=ERF 

****************************************************** 


PURPOSE 


USAGE 


-  EVALUATE  THE  ERROR  FUNCTION 

-  RESULT  =  ERF(Y) 


ARGUMENTS  Y  -  INPUT  ARGUMENT  OF  THE  ERROR  FUNCTION. 

ERF  -  OUTPUT  VALUE  OF  THE  ERROR  FUNCTION. 

COPYRIGHT  -  1978  BY  IMSL,  INC.  ALL  RIGHTS  RESERVED. 

WARRANTY  -  IMSL  WARRANTS  ONLY  THAT  IMSL  TESTING  HAS 

APPLIED  TO  THIS  CODE.  NO  OTHER  WARRANTY 
EXPRESSED  OR  IMPLIED,  IS  APPLICABLE. 

****************************************************** 


REAL  FUNCTION  ERF(Y) 


REAL 

INTEGER 

DIMENSION 

REAL 


SPECIFICATIONS  FOR  ARGUMENTS 
Y 

SPECIFICATIONS  FOR  LOCAL  VARIABL 
I SW  I 

P(5),Q(3),P1(8),Q1(7),P2(5),Q2(4) 
P,Q,P1,Q1,P2,Q2, XMIN, XLARGE, SSQPI , 


o  o 
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+ 

c 

.477 

DATA 

+ 

+ 

+ 

+ 

DATA 

+ 

+ 

DATA 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

DATA 

+ 

+ 

+ 

+ 

+ 

+ 

C 

DATA 

+ 

+ 

+ 

+ 

DATA 

+ 

+ 

+ 

C 

DATA 

DATA 

C 

X  =  Y 
ISW  =  1 


RES , XSQ , XNUM, XDEN, XI , X 
COEFFICIENTS  FOR  0.0  .LE.  Y 

P( I )/- . 44422647396874/, 
P(2)/10. 731707253648/, 

P( 3 )/15. 915606197771/, 
P(4)/374. 81624081284/, 

P(5)/2 . 5612422994823E-02/ 
Q(l)/17. 903143558843/, 

Q( 2 )/124. 82892031581/, 

Q( 3 )/332 . 17224470532/ 
COEFFICIENTS  FOR  .477  .LE.  Y 
.LE.  4.0 

PI (l)/7. 2 117582508831/, 
Pl(2)/43. 162227222057/, 

PI ( 3 )/152 . 98928504694/, 
Pl(4)/339. 32081673434/, 

PI ( 5 )/451. 91895371187/, 
Pl(6)/300. 45926102016/, 

PI ( 7 )/-l . 3 6864857382 72E-07/, 

PI (8)/. 56419551747897/ 
Ql(l)/77. 000152935229/, 
Ql(2)/277. 58544474399/, 

Q1 (3 )/638 . 98026446563/, 
Ql(4)/931. 35409485061/, 

Q1 ( 5 )/790 . 95092532790/, 

Q1 ( 6 )/300 . 45926095698/, 

Ql( 7 )/12. 782727319629/ 
COEFFICIENTS  FOR  4.0  .LT.  Y 
P2 ( 1 )/- . 22695659353969/, 

P2 (2 )/-4. 947309 1062325E-02/, 

P2 ( 3 )/-2 . 9961070770354E-03/, 

P2 ( 4)/-2 . 2319245973418E-02/, 

P2 ( 5 )/-2 . 7866130860965E-01/ 

Q2 ( 1 )/l . 0516751070679/, 

Q2 ( 2 )/. 19130892610783/, 

Q2 ( 3 )/l . 0620923052847E-02/, 
Q2(4)/l. 9873320181714/ 
CONSTANTS 

XMIN/1 . 0E- 8/ , XLARGE/5 . 6875E0/ 
SSQPI/. 56418958354776/ 

FIRST  EXECUTABLE  STATEMENT 


IF  (X.GE.0.0E0)  GO  TO  5 
ISW  =  -1 
X  =  -X 

5  IF  (X. LT. . 477E0)  GO  TO  10 
IF  (X.LE.4.0E0)  GO  TO  25 
IF  ( X . LT . XLARGE )  GO  TO  35 
RES  =  1.E0 


.LT. 
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GO  TO  50 


ABS(Y)  .LT.  .477,  EVALUATE 
APPROXIMATION  FOR  ERF 
10  IF  (X.LT.XMIN)  GO  TO  20 
XSQ  =  X*X 
XNUM  =  P( 5 ) 

DO  15  1=1,4 

XNUM  =  XNUM*XSQ+P(I) 

15  CONTINUE 

XDEN  =  ( (Q( 1 ) +XSQ) *XSQ+Q(2 ) ) *XSQ+Q(3 ) 

RES  =  X* XNUM /XDEN 
GO  TO  50 

20  RES  =  X*P<4)/Q(3) 

GO  TO  50 

.477  .LE.  ABS(Y)  .LE.  4.0 
EVALUATE  APPROXIMATION  FOR  ERF 


25  XSQ  =  X*X 

XNUM  =  P1(7)*X+P1(8) 
XDEN  =  X+Ql ( 7 ) 

DO  30  1=1,6 

XNUM  =  XNUM*X+P1 ( I ) 
XDEN  =  XDEN*X+Q1(I) 
30  CONTINUE 

RES  =  XNUM /XDEN 
GO  TO  45 


4.0  .LT.  ABS(Y),  EVALUATE 
APPROXIMATION  FOR  ERF 

35  XSQ  =  X*X 

XI  =  1 . OEO/XSQ 
XNUM  =  P2(4)*XI+P2(5) 

XDEN  =  XI+Q2(4) 

DO  40  1=1,3 

XNUM  =  XNUM*XI+P2 ( I ) 

XDEN  =  XDEN*XI +Q2 ( I ) 

40  CONTINUE 

RES  =  ( SSQPI+XI *XNUM/XDEN)/X 
45  RES  =  RES*EXP( -XSQ) 

RES  =  1 . OEO-RES 
50  IF  (ISW.EQ.-l)  RES  =  -RES 
ERF  =  RES 
RETURN 
END 
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APPENDIX  B 

v  GLOSSARY  OF  SELECTED  VARIABLES 

^  SUSPENDED  SOLIDS  PROGRAM 

1.  CO  =  percentage  of  individual  size  fraction. 

2.  CONOR  =  original  source  concentration  (mg/1). 

3.  DELX  =  longitudinal  step  size  (m). 

4.  DELY  =  vertical  step  size  (m). 

5.  DELZ  =  transverse  step  size  (m) . 

6.  EXPCO  =  spreading  factor  (decreasing  EZ). 

7.  EZ  =  transverse  mixing  coefficient  (m2/21//s). 

8.  E(y)  =  function  for  determining  vertical  dispersion 

coefficient. 

9.  H  -  river  depth  (m). 

10.  NSEDF  =  number  of  sediment  fractions. 

11.  NSTEP  =  number  of  longitudinal  steps  (usually  50). 

12.  THETA  =  model  stability  factor  (set  equal  to  1.0). 


13.  U  =  river  velocity  (m/s). 


14.  W  =  particle  settling  velocity  (m/s). 


15.  X  =  downstream  coordinate. 

16.  XL  =  source  half-width  (m) . 

17.  Y  =  vertical  coordinate. 

18.  Z  =  transverse  coordinate. 


174 


APPENDIX  C 

CHEMICAL  PLUME  AND  DILUTION  COMPUTER  MODEL 
USER'S  MANUAL 


The  User's  Guide  provides  information  that  should  help  users  to  run  the 
Chemical  Plume  and  Dilution  computer  models. 


Introduction 

The  program  'CPLUME'  refers  to  Chemical  plume  model  and  'DILUTION'  to  the 
Dilution  model.  These  two  programs  were  written  to  simulate  a  chemical  plume 
resulting  from  a  dredging  operation  in  a  river  system.  They  use  a  subroutine 
MDNORD  (Normal  probability  distribution  function  of  a  double  precision 
argument)  to  calculate  the  normal  distribution  function  in  the  models.  The 
subroutine  MDNORD  is  available  from  IMSL  (International  Mathematical  and 
Statistical  Library)  at  the  Weeg  Computing  Center  (WCC)  on  the  IBM  370,  the 
Vax,  and  Prime  750  systems.  The  models  were  programmed  in  the  FORTRAN 
language  on  a  Prime  750  computer  at  The  University  of  Iowa,  Iowa  City,  Iowa. 

Algorithms 

The  computer  programs  were  set  up  in  the  following  manner. 

CPLUME 

First,  the  program  reads  in  all  the  data  from  the  data  file  needed  to 
solve  the  equation.  Next,  the  computer  calculates  the  suspended-solids 
concentration  at  a  given  point  (z,x).  Using  the  suspended-solids 
concentration,  it  calculates  the  chemical  concentrations  in  the  plume  and 
writes  out  the  result.  A  flow  chart  and  a  program  listing  follow  as  well  as 
example  input  and  output  files,  and  an  input  and  output  variable  dictionary. 


Flow  Chart  of  CPLUME 


READ 


D,UX,EY,B,BC,MCIN 
BPSD,BPST,BPCY ,PSD 
PST,PCY,WSD,WST,WCY 
BGC,CIN,KP1,KP2,KP3 
SK,YMAX,YS 


o  o  o 
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Listing  of  the  Computer  Program  "CPLUME" 


CPLUME 

C 

C  'CPLUME* 

C  THE  PROGRAM  'CPLUME'  WAS  WRITTEN  TO  CALCULATE 
C  THE  CHEMICAL  CONCENTRATION  IN  PLUME  RESULTED 
C  FROM  A  DREDGING  OPERATION 
C  THE  MODEL  SIMULATES  THE  GREAT  III  REACH  IN 
C  THE  MISSISSIPPI  RIVER  <RM  170) 

C 

C 

DIMENSION  YS(20)*Yl(20)»Y2<20)»r(20)»Pl(20)*P2(20)»DPSD(20)» 

*C ( 20 ) » CT ( 20 ) >  CT  A ( 20 ) >  BP  <  20 ) >  Li  .  20) >CP2<20) »CP3(20) > CSIK20) » 
*CST(20>  >DP3Y(20) >CCY(20) >TC(20) >DPCY(20;  >M1 (23) >n2(20) *n3(20>  > 
HXKPM1 (20) »XKPM2(20) *XnPM3(20) »XXt (20) » <KPM(20) >A1(2>/') 

DIMENSION  UKM ( 20 ) >  UKM1 <  20  >  >  UKM2 ( 20  >  >  WKM3 ( 20  >  >  E ( 20  >  > DPST ( 20 ) > 
*TCSD(20) >TCST(20) >TCCY(2Q) >TTC(20) > TM1(20) >TM2(20) »TM3(20) > 
*XKPTM1(20) »XKPTM2<20) »XKPTM3(20) >X11 (20) >X21 (20) »X31 (20) 

C 

c 

REAL*8  Xf EY*SK1 »SK7*SK3»M1 »M2»M3»KP1»KP2»KP3 
REAL*8  Pll» P7I » Y1 1 > Y2I 

REAL*8  BCSD  >BCST»BCCY>  BPSD » BE'ST > BPCY » BC > BGC » MC IN 
INTEGER  YHAX>XMAX>CF 

SET  DATA 

DATA  CF/ 1000000.0/ 

DATA  INi IOUT/11 >10/ 

C 

READ(IN>»)  D>UX>EY 
READ(IN>«)  B 
READ(IN>«)  BCfMCIN 
READ(IN>*)  BPSD >8FST> BPCY 
READ(IN>*)  PSDfPSTvPCY 
READ(IN>*)  USPf'JSr >WCY 
READdNr*)  BGC > C I N 
REAIK IN>*)  NR1>I\F'2>T.F*3 
READ(IN>*)  SK 


I 
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C  COMPUTE  BACKGROUND  CONCENTRATION'S  OF  SAND  (BCSD), 

C  SILT  (BCST)  AND  CLAY  (BCCY).  BC  IS  THE  BACKGROUND 

C  CONCENTRATION  OF  SUSPENDED  SOLIDS*  BPSB*  BF'ST *  AMD 
C  BPCY  ARE  FRACTION  OF  SAND*  SII  T  AND  CLAY*  RESPECTIVERY. 

C 

BCSD=BC*BF'SD 

BCST=BC*BPST 

8CCY=BC*BPCY 

C 

C  READ  NUMBER  OF  INPUT  DATA*  YHAX 
C 

READ! IN,*)  YMAX 
C 

C  READ  YS,  Z/SIGMA  Z 
C 

READdN,*)  (YS<I), 1=1, YMAX) 

WRITE( IOUT ,6001 )  (I*I=1*YHAX) 

6001  FORMAT (//'INPUT  YS'//20( I2*2X) ) 

WRITE! IOUT *6000)  ( YS(I )* 1  =  1  * YMAX) 

6000  F0RMAT(20F4. 1 ) 

C 

C  SET  TABLE  FORMAT  FOR  THE  FINAL  RESULT 
C 

WRITE (I OUT *207) 

207  FORMAT (IX*///) 

WRITE (IOUT, 204) 

204  FORMAT (1X,T3*'(1),,T12*'(2)'*T19*'(3)'* T27* ' (4) ' *T35* ' (5) ' * 

*T43, ' (6) '  ,T50*  '(7)',I58,,(8)'»T66,'(9)','l74,'(10)'  *T83* ' ( 11 ) '  * 
*T92, ' ( 12) ' *T100* ' ( 13) '//) 

WRITEdOUT  ,202) 

202  FORMATdX*  'DISTANCE'  *Tll  *  'WIDTH'  ,T19,  'SAND'  *T27* 

*'SILT' * T35* 'CLAY' ,T42, 'TOTAL' *T50* 'TSS' *T58* 'POLL. ' *T66, 'POLL. ' * 
*T74, 'POLL. ' *T81 , 'DISSOLVED' *T92* 'TOTAL ' *T100, 'POLL. ' ) 

WRITEdOUT, 203) 

203  FORMAT ( 1X»T50» 'ABOVE' , T59, ' ON' ,T67, 'ON' ,T75» 'ON' ,T100, 'ABOVE' ) 
WRITE (IOUT, 201) 

201  FORMATdX, T5,'X',Tt2,'Z',T16,'  ',T24,'  ',T32,'  ',T42,'SS.'* 

*T50, 'AMBI. ' , 

*T58, 'SAND' ,T66, ' S I L  f ' »T74, 'CLAY' ,T83, 'POLL. ' ,T92, 'POLL. ' , 
iTIOO, 'AMBI . ' ) 

WRITEdOUT, 205/ 

205  FORMAT! 1X,T4, ' <M) ' ,Tll,'(M)',T19, ' (MG/L) ' ,T27, ' (MG/L) ' , 

*T35, ' (MG/I.) ' ,T42 , ' (MG/L ) ' ,1 50, ' (M6/L) ‘ , T58, ' (MG/L) ' ,T66, ' (MG/L) ' , 
tT74, ' (MG/L) ' ,T02, ' (MG/L)' ,T92, ' (MG/L) ' ,T100, ' (MG/L)' 

*///) 

C 


* 
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c 

C  1.  CALCULATION  OF  SAND,  SILT  AND  CLAY  CONCENTRATIONS  ABOVE  AMBIENT 
C  AT  GIVEN  LOCATION  (Z»X).  SEE  GREAT  II  REPORT  (SCHNuOR  ET  AL.  1980) 

C 

C  PRELIMINARY  CALCULATION 
C 

210  READ(IN»*)X 

IF(X  .LT.  .0)  GO  TO  999 
C 

A=X/(D*UX) 

ESD=FXP(-WStiA) 

EST=EXP(-WST*A> 

ECY=EXP(-WCY*A> 

s=<2.*ey*x/ux)**< 

DO  400  1-1 , YMAX 
Y(I)=YS(I)*S 
Y1(I)  =  (Y(I  H-B'/S 
Y2(I)=(Y<I)-B)/S 
Y 1 I =Y 1 < I ) 

Y2I=Y2( I ) 

CALL  SUBROUTINE  * MDNORD '  FROM  THE  IMSL  LIBRARY  TO  CALCULATE 
NORMAL  DISTRIBUTION  FUNCTION 

CALL  MDNORD  (Y1I,P1I) 

CALL  MDNORD  <Y2I»P2I) 

F'l  < I  )=F'l  I 
P2( I )=P2I 
riP<I>=Pl<I)-P2(I) 

DPSD< I ) -DP( I )*ESD 
DPST(I)*DP{I)*EST 
DPCY ( I )=DF( I ) *ECY 

COMPUTE  SUSPENDED  SOLIDS  CONCENTRATION  AF-OVE  AMBIENT 

CSD(I?-DP5D(T)»=SD*nCIN 
CST  < I )-DPST ( I )*PST*MCIN 
CCY  < I )  =  DPCY ( I >*FTY*MCIN 
TC( I )=CSD< I H CST( I ) TCCYt I ) 

COMPUTE  TOTAL  SAND.  SILT  AND  CLAY  CONCENTRATION  (AMBIENT  PLUS  PLUME) 

TCSD( I )=CSD( I )tBC3D 
TCST ( I ) =CST  < I HPf.ST 
TCCY  ( I  )=CCY  ( I  HBCCY 
C 

C  COMPUTE  TOTAL  SUSPENDED  SOLIDS  CONCENTRATION  (AMBIENT  PLUS  PLUME) 

C 

TTC< I )=TC5D( I ) +TC5T  < I )+TCCY( I ) 


A120  725  REFINEMENT  AND  VERIFICATION  OF  PREDICTIVE  MODELS  OF 
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c 

C  2.  CALCULATION  OF  CHEMICAL  CONCENTRATION  AT  GIVEN  (Z.X)  POINT 

C 

C 

C  COMPUTE  SAND*  SILT  AND  CLAY  CONCENTRATION  IN  KG/L  (ABOVE  AMBIENT) 
C  CF»CONVtRSION  FACTOR  TO  CHANCE  M  FROM  MG  TO  KG 
C 

Mi<n=csnm/CF 

M2(I)=CST(I)/CF 

M3U)*CCY(I)/CF 

C 

C 

C  CONVERT  TOTAL  SAND*  SILT  AND  CLAY  CONCENTRATION  FROM  MG/L  TO  KC/L 
C 

TM1 ( I )*TCSD( I )/CF 
TM2(I)*TCST(I)/CF 
TH3(I)*TCCY(I)/CF 
C 

XKPM1 ( I ) =KP1*M1 ( I ) 

XKPM2( I )*KP2tM2( I ) 

XKPM3( I )=KP3*M3( 1 ) 

C 

XKPM(I)=1 .  +XKPM1 ( I )+XKPM2( I )+XKPM3( I ) 

C 

UKM1 ( I )»(HSD/4.  MtXKPMl ( I ) 

UKH2( I ) ■ ( WST/4 . ) *XKPM2< I ) 

MKM3(I)*(UCY/4«)4XKPN3(I) 

UKM(I) “HKM1 ( I ) +WKM2 ( I ) +WKN3 ( IHSK 
A1(I)«HKM<I)/XKPM(I) 

E(I)*EXP<* A1(I)*X/UX) 

C 

C  COMPUTE  CHEMICAL  CONCENTRATION  ABOVE  AMBIENT  (CTA) 

C 

CTA( I )*CIN$DP< I )*E ( I ) 

C 

C  COMPUTE  TOTAL  CHEMICAL  CONCEN1 RATION  (AMBIENT  4  PLUME) 

C 

CT(I)“CTA(I)+BGC 

C 

XKPTM1 ( I )*KP1$TM1 ( I ) 

XKPTM2(I)=KP2*TM2(I) 

XKPTM3( I)*KP34TM3<  I ) 

C 

XXI ( I ) «1 .  +XKPTM1 ( I ) +XKPTM2 ( I ) +XKPTM3 ( I ) 

C 

Xn(I)«XKPTHl(I)*CT(I> 

X21 ( I ) *XKPTM2 ( I ) tCT  < I ) 

X31 ( I )*XKPTM3 ( I ) *CT ( I ) 


uuu  uuuu 


COMPUTE  CPI'  CP2»  CP3'  AND  C 

CP1(D*X11(I)/XX1(I) 

CP2( I )=X21 ( I ) /XXI  ( I ) 

CP3( I )=X31 ( I )/XXl { I ) 

C< I )=CT< I >-<  CPI < I )+CP2< I >+CP3< I ) > 

400  CONTINUE 

URITE  OUT  RRESULT 

DO  300  I«1»YMAX 

URITE ( IOUT » 301 ) X » Y  < I ) » TCSD( I ) » TCST  < I ) » TCCY ( I > » TTC( I ) » TC ( I ) » 
tCPl(I)'CP2(I)'CP3(I)iC(I)'CT<I)'CTA(I) 

301  FORMAT<lX'F7.0.T8»F7.1.T14.F;.2'T24»F7.2'T32'F/.2f 
*T40»F7.2'T48»F7.2»T56»F7.2»T64*F7.2»T72'F7.2»T81'F7.2» 
tT90»F7.2»T98»F7.2) 

300  CONTINUE 

URITE ( IOUT' 109) 

109  FORMAT (IX,//) 

END  OF  COMPUTATION  AT  A  GIVEN  DOWNSTREAM  DISTANCE  X 
READ  NEU  X 


GO  TO  210 
999  STOP 


I 


Example  of  Input  and  Output  Files  -  CPUHE 


a.  Input  Data 

4.0  0.8  2S.0 
130.0 
82.0  S.O 
0.0  0.50  0.S0 
0.70  0.15  0.15 
0.02  0.0003  0.000091 
0.03  0.15 
13.0  130.0  1300.0 
0.0090O00901 
18 

0.0  0.2  0.4  0.4  0.8  1.0  1.1  1.2  1.3  1.4  1.4  1.8  2.0  2.5  3.0  3.5  4.0  5.0 
20.0 
50.0 
100.0 
250.0 
-0 

b.  Output  File 
INPUT  vs 

1  2  3  4  5  4  7  8  0  10  11  12  13  14  15  14  17  18 

0.0  0.2  0.4  0.4  0.8  1.0  1.1  1.2  1.3  1.4  1.4  1.8  2.0  2.5  3.0  3.5  4.0  S.O 


(1)  (2)  (3)  <4>  (5)  (4)  (7)  (8)  (9)  (10)  (11)  (12 

NSTMiCE  HIDTH  SMS  Sill  041  TOTAL  TSS  POLL.  POll.  POLL.  OISSOLVED  TOT 


75  84.5? 

75  84.5? 
75  84>S8 
75  84.58 

/S  84.58 
75  86.5/ 

75  84.54 

75  84.54 

71  84.55 

74  84.54 

74  84.50 

7J  84.45 

a  fcg 

54  85.44 

4$  84.42 


09  0. 

00  0. 
00  0. 
00  0. 
00  9. 
00  9. 


00  0. 
CO  0. 
00  0. 

CO  o. 


Cl  0 

13  4 

.94  0. 

Cl  0 

13  0 

.04  0. 

Cl  0 

13  0 

.04  0. 

Cl  0 

13  0 

.04  0. 

Cl  0 

13  0 

.04  0. 

Cl  0 

13  0 

.04  0. 

01  0 

13  fl 

.04  0. 

01  0 
01  0 

13  0 

13  0 

.04  0. 

.04  9. 

01  0 

01  0 

01  0 

SI  l 

13  0 

13  0 

.04  0. 

.04  0. 

.04  0. 

.04  0. 

.04  0. 

01  0 

10  0 

.03  0. 

01  0 

08  0 

.03  0. 

01  0.94  I 

1.02  0. 

00  0.93  « 

.01  0. 
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Chemical  Plume  Model  Input  Variables  and  Format 


arlable 

Format 

Units 

Definition 

B 

* 

m 

Width  of  the  plume  source ,  b 

D 

* 

m 

Mean  stream  depth  In  plume,  D 

UX 

* 

m/sec 

Mean  stream  velocity,  Ux 

EY 

* 

m^/sec 

Dispersion  coefficient,  Ez 

WSD 

* 

m/sec 

Mean  settling  velocity  for  sand 

WST 

* 

m/sec 

Mean  settling  velocity  for  silt 

WCY 

* 

m/sec 

Mean  settling  velocity  for  clay 

PSD 

* 

None 

Fraction  of  sand 

PST 

* 

None 

Fraction  of  silt 

PCY 

* 

None 

Fraction  of  clay 

BC 

* 

mg/1 

Ambient  suspended-solids  concentration 

MIN 

* 

mg/1 

Initial  suspended-solids  concentration 
above  ambient 

BSPD 

* 

None 

Fraction  of  ambient  sand 

BPST 

* 

None 

Fraction  of  ambient  silt 

BPCY 

* 

None 

Fraction  of  ambient  clay 

BGC 

* 

mg/1 

Ambient  chemical  concentration 

CIN 

* 

mg/1 

Initial  chemical  concentration  above 

ambient 

KP1 

* 

(ug/l)(wg/kg)-1  Partitioning  coefficient  for  sand 

KP2 

* 

(yg/lHug/kg)”1  Partitioning  coefficient  for  silt 

KP3 

* 

(wg/lMug/kg)"1  Partitioning  coefficient  for  clay 

SK 

* 

(sec)"1 

Total  decay  rate  of  the  chemicals 

YMAX 

* 

None 

Number  of  Input  data,  yflz 

YS 

* 

m 

Distance  coordinate  In  the  lateral 
dl recti on/o2 

X 

* 

m 

Distance  coordinate  In  the  longitudinal 

direction 


*  Free  format 


Chemical  Plume  Model  Output  Variables  and  Format 


Variable 

Format 

Unit 

Definition 

X 

F7.0 

m 

Distance  coordinate  in  the  longitudinal 
direction,  x 

Y 

F7.1 

m 

Distance  coordinate  In  the  lateral 

direction,  z 

TCSO 

F7.2 

mg/1 

Total  sand  concentration  In  bulk  water 

TCST 

F7.2 

mg/1 

Total  silt  concentration  In  bulk  water 

TCCY 

F7.2 

mg/1 

Total  clay  concentration  In  bulk  water 

TC 

F7.2 

mg/1 

Total  suspended-solids  concentration 

CPI 

F7.2 

mg/1 

Chemical  concentration  adsorbed  onto  sand 
particulates  In  bulk  water 

CP2 

F7.2 

mg/1 

Chemical  concentration  adsorbed  onto  silt 
particulates  in  bulk  water 

CP3 

F7.2 

mg/1 

Chemical  concentration  adsorbed  onto  clay 
particulates  In  bulk  water 

C 

F7.2 

mg/1 

Dissolved  chemical  concentration 

CT 

F7.2 

mg/1 

Total  chemical  concentration  In  bulk  water 

CTA 

F7.2 

mg/1 

Total  chemical  concentration  above  ambient 

Procedure 

The  following  procedure  and  an  example  show  how  to  run  the  program 
'CPLUME*.  In  this  manual  user  Inputs  are  underlined,  and  computer  prompts  and 
responses  are  capitalized  (see  example  run). 

The  PRIME  750  computer  responds  with  OK  when  login  is  accepted. 

Step  1.  Open  a  file  'DATA*  to  read  In  data  needed  for  conputation. 

Step  2.  Open  a  file  'REST'  to  store  the  result  after  computation. 

Step  3.  Run  the  program  'CPLUME'. 

Step  4.  Close  all  files. 

Step  5.  List  the  computed  result  from  the  file  'REST*. 


Chemical  Plume  Model  Output 


(1)  Distance  coordinate  In  the  longitudinal  direction. 

(2)  Distance  coordinate  In  the  lateral  direction. 

(3)  Concentration  of  sand  at  a  given  point  (z,x). 

(4)  Concentration  of  silt  at  a  given  point  (z,x). 

(5)  Concentration  of  clay  at  a  given  point  (z,x). 

(6)  Total  suspended- sol  ids  concentration  at  (z,x). 

(7)  Suspended-solids  concentration  above  ambient  at  (z»x). 

(8)  Concentration  of  chemical  pollutant  adsorbed  Into  sand 

(9)  Concentration  of  chemical  pollutant  adsorbed  Into  silt 

(10)  Concentration  of  chemical  pollutant  adsorbed  Into  clay 

(11)  Concentration  of  dissolved  chemicals. 

(12)  Total  concentration  of  chemical  pollutant. 

(13)  Total  chemical  concentration  above  ambient. 
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Example  Run 

OK*  0  PAT A  7  1 
OK*  0  REST  6  2 
OK*  SEG  ♦Cf'LUilE 

m*  stop 

OK*  C  ALL 
OK*  SLIST  REST 

INPUT  VS 

1  2  3  4  5  6  7  8  9  1G  11  12  13  14  IS  1&  17  18 
0.0  0.2  0.4  0.4  0.8  1.0  1.1  1.2  1.3  1.4  1.4  1.8  2.0  2.5  3.0  3.5  4.0  5.0 


(1) 

(2) 

(3) 

(4) 

(5) 

(4) 

(7) 

(8) 

(91 

(10) 

(11) 

(12) 

(13) 

DISTANCE  WIDTH 

SAND 

SILT 

CLAY 

TOTAL 

TSS 

POLL. 

POLL. 

POLL, 

DISSOLVED 

TOTAL 

POLL. 

• 

ABOVE 

ON 

ON 

ON 

ABOVE 

X 

Z 

SS. 

AHBI. 

SAND 

SILT 

CLAY 

POLL. 

POLL. 

AH6:, 

(N) 

(H) 

(H6/L) 

(HG/L) 

(N6/L)  (H6/L) 

(HG/L) 

(HG/L) 

(HC/L) 

(NS/L) 

(HG/L) 

(HG/L) 

(HG/L) 

20 

0.0  . 

3.09 

41.75 

41.75 

84.59 

4.59 

0.00 

0.01 

0.13 

0.04 

0.18 

0.15 

20 

7.1 

3.09 

41.75 

41.75 

84.59 

4.59 

0.00 

0.01 

0.13 

0.04 

0.18 

0.15 

20 

14.1 

3.09 

41.73 

41.75 

84.58 

4.58 

0.00 

0.01 

0.13 

0.04 

0,18 

0.15 

20 

21.2 

3.09 

41.75 

41.73 

84.1*8 

4.58 

0.00 

0.01 

0.13 

0.04 

0.18 

0.15 

20 

28.3 

3.08 

41.75 

41.75 

84.58 

4.56 

0.00 

0.01 

0.13 

0.04 

0.18 

0.15 

20 

35.4 

3.06 

41.75 

41.75 

84.57 

4.57 

0.00 

C.01 

0.13 

0.04 

0.18 

0.15 

20 

38.9 

3.07 

41.74 

41.75 

84.34 

4.54 

0.00 

0.01 

0.13 

0.04 

0.18 

0.15 

20 

42.4 

3.07 

41.74 

41.75 

84.54 

4.54 

0.00 

0.01 

0.13 

0.04 

0.18 

0.15 

20 

44.0 

3.04 

41.74 

41.74 

84.55 

4.55 

0.00 

0.01 

0,13 

0.C4 

0.18 

0.15 

20 

49.5 

3.05 

41.74 

41.74 

84.54 

4.54 

0.00 

0.01 

0.13 

0.04 

0.18 

0.15 

20 

54.4 

3.03 

41.73 

41.74 

84.50 

4.50 

0.00 

0.01 

0.13 

0.04 

o.ia 

0.15 

20 

43.4 

3.00 

41.73 

41.73 

84.45 

4.45 

0.00 

0.01 

0.12 

0.C4 

0.1S 

0.15 

20 

70.7 

2.94 

41.71 

41.71 

84.37 

4.37 

0.00 

0,01 

0.12 

0.04 

0.17 

0.14 

20 

88.4 

2.72 

41.44 

41.44 

84.04 

4.04 

0.00 

0.01 

0.11 

0.04 

0.14 

0.13 

20 

104.1 

2.32 

41.54 

41.54 

85.44 

3.44 

0.00 

0.01 

0.10 

0.03 

0.14 

0.11 

20 

123.7 

1.74 

41.43 

41.43 

84.42 

2.42 

0.00 

0.01 

0.08 

0.03 

0.12 

*•44 

0.09 

20 

141.4 

1.15 

41.28 

41.28 

83.71 

1.71 

0.00 

0.01 

0.04 

0.02 

0.09 

0.06 

20 

174.8 

0.29 

41.07 

41.07 

82.43 

0.43 

o.c-o 

0.00 

C.C'3 

0.01 

0.04 

0.01 
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DILUTION 

First,  the  program  reads  In  all  the  data  from  the  data  file  needed  to 
solve  the  equation.  Next,  the  computer  calculates  the  dilution  ratio  at  a 
given  point  (z,x).  Using  the  given  Initial  concentration  of  a  chemical 
pollutant.  It  calculates  the  dissolved  chemical  concentration  In  the  plume, 
and  writes  out  the  result. 

A  flow  chart  and  a  listing  of  the  program  follow  as  well  as  example  Input 
and  output  data  files,  and  an  Input  and  output  variable  dictionary. 


oor>  onn  r>oo  o  noo  n  no  o  n 


Listing  of  the  Computer  Program  "DILUTION" 


DILUTION 


THE  PR08RAM  'DILUTION'  WAS  WRITTEN  TO  COMPUTE 
THE  DILUTION  RATIO  ANO  A  DISSSOLVED  POLLUTANT 
CONCENTRATION  IN  A  PLUME 

THE  MODEL  SIMULATES  THE  GREAT  III  REACH  IN  THE 
MISSISSIPPI  RIVER  (RM  170) 

DIMENSION  YS(20) >Y1(20)»Y2(20) »Y!20) fP1(20)»P2(20) » 
*DP(20)fCD<20)»DV<20) 

REAL68  P1I»P2I»Y1I»Y2I»X»EY 
INTEGER  YMAX 

DATA  1N»I0UT/11»10/ 

READ  DATA 

READ<IN>«)  DtUXrEY 
READ(IN>S)  B 
READ(IN»*)CIN 

READ  NUMBER  OF  INPUT  DATA*  YMAX 

READ(IN»«)  YMAX 

READ  YS»  Z/SIOMA  Z 

READ(IN»f )  (YS(I) • 1=1 » YMAX ) 

WRITE (IOUT* 6001)  ( It  1=1  * YMAX) 

6001  FORMAT!//' INPUT  YS  '//20!I2»2X> > 

WRITE! IOUT, 6000)  !YS!I>rW»YHAX> 

6000  F0RMAT!20F4«1) 

C 

WRITE !IOUT* 207) 

207  FORMAT! IX »///) 


oooooo  o  oooo  ooooo  o  o  non 
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SET  TABLE  FORMAT  FOR  THE  FINAL  RESULT 
WRIT£(I0UT*2C4) 

204  FORMAT ( IX »T3» ' (1)' #T13» ' (2) ' »T27* ' *3) ' *T42* ' (4) '» 
tT57» '  (5) ' *T71» ' *6) '// ) 

WRITE (IOUTf 201) 

201  FORMAT* 1X*T4* 'X' *T14* 'Z' »T2A* 'F<+*) ' »T41 * 'F(-t) ' * 

*T55, 'DILUTION' *T6H,'  ') 

WRITE* I0UT>202) 

202  FORMAT* 1X*T55. 'RATIO' *T<48* 'POLL.  C0NC.') 

WRITE* ICUT*203) 

203  FORMAT ( IX *T3* ' (M) ' *T13» ' (M) ' *T70* ' (MB/L) ' »//) 


READ  DOWNSTREAM  DISTANCE*  X 


1  READ* IN»*)X 
IF(X  .LT.  .0)  CO  TO  999 
S=<2.*EY*X/UX)**(.5) 

DO  200  Ial»YMAX 
Y<I)=YS<J)*S 
Y1(I)=<Y(I)+B)/S 
Y2(I)=(Y*I)-B)/S 
Y1I=Y1*I) 

Y2I=Y2(I) 

CALL  SA6R0UTINE  MDNORD  FROM  LIBRARY  IMSL 
TO  COMPUTE  THE  NORMAL  DISTRIBUTION  FUNCTION 

CALL  MDNORD  (Y1I.P1I) 

CALL  MDNORD  *Y2I»P2I) 

P1<I)*P1I 
P2< I )=P2I 
DP(I)=P1(I)-P2(I> 

COMPUTE  THE  DILUTION  RATIO 

DV(I)=1./DP(I) 

COMPUTE  DISSOLVED  POLLUTANT  CONCENTRATION 

CD(I)*CIN/DV(I) 

C 

200  CONTINUE 
C 


C  WRITE  OUT  COMPUTED  RESULT 
C 

DO  300  1=1 *YMAX 

WRIT£<I0UT*301)X*Y(I)*P1<I)*P2<I)*DV<I)*CD(I) 

301  FORMAT (iX*F5.0*T10*F7»l*T20*F10«4*T35*F10»4*T50*F10«l 
*TA5*F10.2> 

300  CONTINUE 

WRITE  UOUT*  10?) 

109  FORMAT ( IX »//) 

C 

C  END  OF  COMPUTATIOPN  OF  DISSOLVED  POLLUTANT 
C  CONCENTRATIONS  AT  THE  CROSS  SECTION  QF  THE 
C  STREAM  AT  GIVEN  DOWNSTREAM  DISTANCE*  X. 

C 

C  READ  NEXT  DOWNSTREAM  DISTANCE*  X 
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Example  of  Input  and  Output  Files  -  DILUTION 


a.  Input  Data  File 

4.0  0.8  25.0 
130.0 
0.15 
18 

0.0  0.2  0.4  0.4  0.8  1.0  1.1  1.2  1.3  1.4  1.4  1.8  2.0  2.5  3.0  3.5  4.0  5.0 
20.0 
50.0 
100.0 
250.0 
-9 


b.  Output  File 

INPUT  YS 

1  2  3  4  5  4  7  8  9  10  11  12  13  14  15  14  17  18 

0.0  0.2  0.4  0.4  0.8  1.0  1.1  1.2  1.3  1.4  1.4  1.8  2.0  2.5  3.0  3.5  4.0  5.0 


(1) 

(2) 

(3) 

(4) 

(5) 

(4) 

X 

Z 

F<«) 

F(-J) 

DILUTION 

RATIO 

POLL.  CONC 

<H) 

(«) 

(N6/L) 

20. 

0.0 

0.9999 

0.0001 

1.0 

0.15 

20. 

7.1 

0.9999 

0.0003 

1.0 

0.15 

20. 

14.1 

1.0000 

0.0005 

1.0 

0.15 

20. 

21.2 

1.0000 

0.0010 

1.0 

0.15 

20. 

28.3 

1.0000 

0,0020 

1.0 

0.15 

20. 

35.4 

1,0000 

0,0037 

1,0 

0.15 

20. 

38.9 

1.0000 

0.0050 

1.0 

0.15 

20. 

42.4 

1.0000 

0.0044 

1.0 

0.15 

20. 

44.0 

1.0000 

0,0087 

1.0 

0.15 

20. 

49.5 

1,0000 

0.0114 

1,0 

0.15 

20. 

54.4 

1.0000 

0.0189 

1.0 

0.15 

20. 

43.4 

1.0000 

0.0303 

1.0 

0.15 

20. 

70.7 

1.0000 

0.0448 

1.0 

0.14 

20. 

88.4 

1.0000 

0,1194 

1.1 

0.13 

20. 

104.1 

1.0000 

0.2492 

1.3 

0.11 

20. 

123.7 

1.0000 

0,4293 

1.8 

0.09 

20. 

141.4 

1.0000 

0.4247 

2.7 

0.04 

20. 

174.8 

1.0000 

0.9071 

10.8 

0.01 

Dilution  Model 
Input  Variables  and  Format 


Variable 

Definition 

B 

* 

m 

Width  of  the  plume  source,  b 

D 

* 

m 

Mean  stream  depth  in  plume,  D 

UX 

* 

m/sec 

Mean  river  velocity,  ux 

EY 

* 

n^/sec 

Dispersion  coefficient,  Ez 

CIN 

* 

mg/1 

Initial  chemical  concentration  above 

ambient 

Output  Variables  and  Format 

Variable 

Format 

Unit 

Definition 

X 

F7.0 

m 

Distance  coordinate  In  the  longitudinal 

direction 

Y 

F7.1 

m 

Distance  coordinate  In  the  lateral 

direction 

PI 

F7.2 

None 

Normal  distribution  function:  F(*), 

*  > 

P2 

F7.2 

None 

Normal  distribution  function:  F(*), 

*  < 

DV 

F7.2 

None 

Dilution  ratio 

CD 

format 

F7.2 

mg/1 

Dissolved  chemical  concentration 

Dilution  Model  Output 

(1)  Distance  coordinate  In  the  longitudinal  direction 

(2)  Distance  coordinate  In  the  lateral  direction 

(3)  Normal  distribution  function:  F(*)»  *2° 

(4)  Normal  distribution  function:  F(*),  *  <  0 

(5)  Dilution  ratio 

(6)  Dissolved  chemical  concentration 


APPENDIX  D 

REFERENCED  PROCEDURES  FOR  CHEMICAL  ANALYSES 
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Copper 

Iron 


Load 

Manganese 

Mercury 

Nickel 

Zinc 

Total  Solids 

Volatile  Solids 


Referenced  Procedure 

Atonic  Absorption,  Direct  Aspiration  -  Method  213.1, 
Atonic  Absorption,  Heated  Graphite  Furnace  -  Method 
213.2-1.  Methods  for  Chenlcal  Analysis  of  Water  and 
Wastes,  £PA^M)/4-7$-<J26  (March  1574) - 

Atonic  Absorption,  Direct  Aspiration  -  Method  218.1, 
Atonic  Absorption,  Heated  Graphite  Furnace  -  Method 

218.2.  Methods  for  Chenlcal  Analysis  of  Water  and 
Wastes.  EPA-8flfl/4-79-02ti  (March  19797 

Atonic  Absorption,  Direct  Aspiration  -  Method  220.1, 
Methods  for  Chenlcal  Analysis  of  Water  and  Wastes, 
EWTSOT/l-73-020  (March  1379) 

Atonic  Absorption,  Direct  Aspiration  -  Method  236.1, 
Atonic  Absorption,  Heated  Graphite  Furnace  -  Method 

236.2.  Methods  for  Chenlcal  Analysis  of  Water  and 

Wastes.  EPA-fflQ/4-79-02fl  (Kirch  079) - 

Atonic  Absorption,  Direct  Aspiration  -  Method  239.1, 
Atonic  Absorption,  Heated  Graphite  Furnace  -  Method 

239.2.  Methods  for  Chenlcal  Analysis  of  Water  and 

wastes.  EPJUMfl/4-79-0»  (Kirch  W9) - 

Atonic  Absorption,  Direct  Aspiration  -  Method  243.1, 
Methods  for  Chenlcal  Analysis  of  Water  and  Wastes. 
EWUW0/4-79-OH)  (March T979J - 

Cold  Vapor  Method  with  Preliminary  Digestion  - 
Method  245.1.  Methods  for  Chenlcal  Analysis  of 
water  and  Wastes.' EWLgW/T-79-0ZP7*rcK  T979Y  ~ 

Atonic  Absorption,  Direct  Aspiration  -  Method  249.2, 
Methods  for  Chenlcal  Analysis  of  Water  and  Wastes, 

EWC8M74-7WKB  (March  1979] - 

Atonic  Absorption,  Olrect  Aspiration  -  Method  289.2, 
Methods  for  Chenlcal  Analysis  of  Water  and  Wastes, 
(March  T5791 - 

Gravimetric,  Dried  at  103-105*C  -  Method  160.3, 
Methods  for  Chenlcal  Analysis  of  Water  and  Wastes, 
EWEWB/4-7MgB  (March  1979) - 


Gravlnetrlc,  Ignition  at  550*C  -  Method  160.4, 


Referenced  Procedure 


Ammonia  Nitrogen 
Orthophosphate 


Gravimetric,  Separatory  Funnel  Extraction  -  Method 
413.1,  Methods  for  Chemical  Analysis  of  Hater  and 
Wastes.  ENUMOM-79-flM  (Harch  T9TS) - 

Automated,  Phenate  Method  -  Method  417F,  Standard 
Methods,  14th  Edition,  p.  363 

Automated  Ascorbic  Acid  Reduction  Method  -  Method 
424G,  Standard  Methods,  14th  Edition,  p.  422 


Chemical  Oxygen  Demand  Tltrlmetlc,  Mid-Level  -  Method  410.1,  Methods  for 

Chemical  Analysis  of  Hater  and  Hastes,  EPa-S0o/4- 
"(Harch  1979)  - 


Polychlorinated 

Biphenyls  (PCBs)  “Methods  for  Benzidine  Chlorinated  Organic  Compounds, 

Pentachlorophenol  and  Pesticides  In  Water  and 
Wastewater,"  U.S.  EPA  Environmental  Monitoring  and 
Support  Laboratory,  September  1978,  p.  43 


Priority  Pollutants 


EPA  Method  625,  Federal  Register  Vol.  44.  No.  233 
(3  December  1979),  p.  59540 


